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§ 1. INTRODUCTION 


From necessity papers published in any given Journal must conform to certain 
conventions in order to give uniformity to the style of the Journal ; MSS. submitted for 
publication have therefore to be scrutinized and amended before being placed in the 
printer’s hands. : 

The following instructions refer to the main points which have to be considered, and 
authors offering papers to any of the Physical Society’s publications are asked to conform 
to these recommendations. 

In general the author’s degrees or honours will not be printed, but the mode of address, 
e.g., Professor, Dr., Miss, will be inserted. 


§2. Copy 


Manuscripts should be typed in double spacing on paper not wider than 8 in. and not 
longer than 13 in. Only one side of the paper should be used, and a margin of 1-1} in. 
should be left. As alterations in the text cannot normally be allowed once the paper is 
set up in type, authors should aim at absolute clarity of meaning and of typing, and should 
check the typescript carefully before submission. It is advantageous where possible to 
have this read also by a colleague before it is communicated. 


(i) Title and Sectional Headings 


The title should be as concise as is compatible with lucidity ; if it is too long for 
reproduction at the head of right-hand pages, a short title may be given for this purpose. 

Papers are usually divided into sections beginning with § 1. INTRODUCTION and 
ending with §... CONCLUSION. ‘These sections are numbered serially and their titles 
printed in small capitals centrally on the page. If required, sub-section headings may be 
used ; they are printed centrally in italics and if necessary may be numbered in roman 
numerals (1), (ii), ..... Sub-sub-section headings are printed in italics at the beginning 
of the paragraph ; these may be designated (a), (b) .. . if desired. 

At the end of the paper there are usually two further headings: ACKNOWLEDGMENTS 
and REFERENCES. Please note that headings are not part of the text, which should thus 
be complete if headings were omitted. 


(ii) Abstract 


An abstract is printed at the beginning of the paper. It is not part of the paper, and, 
perhaps because its purpose is not well understood, an author’s abstract is often unsatis- 
factory. ‘The abstract is intended for two classes of readers, those directly concerned 
with the subject of the paper and those with a general interest in the subject; for the 
former the abstract should serve as an index, for the latter es asummary. It should 
therefore be designed to state briefly the general aspects of the subject on which new 
information is presented and the main conclusions reached and final numerical results 
obtained. Details of the methods employed, the definition of the field within which the 
results are true, and all the finer considerations involved in the work are matters for the 
paper itself and not for the abstract, 


2 Instructions to authors 


The abstract is placed immediately after the title, name(s) of author(s), and place of 
employment of author(s). It ‘5 thus convenient for the first page of the typescript to 
contain just these items, and for the second page to start with § 1. INTRODUCTION. This 
simplifies assessment of the length of the final printed paper. 

This “title and abstract ” page should be submitted in duplicate. 


(iii) Footnotes 


Footnotes should be used only for brief notes which cannot conveniently be fitted into 
the text, but should be avoided whenever possible. They should be inserted in the MS., 
not at the bottom of the page, but zmmediately below the line to which they refer, and must 


be separated from the rest of the text by horizontal lines across the page above and below. 


(iv) Tables 


Tabular matter should be avoided as far as possible except where there is a gain in 
conciseness by its use. Valuable space can usually be saved by arranging tables with 
more columns than rows. For instance, comparison of monthly values for, say, five years 
is better tabulated : 


Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec. 
1939 
1940 
1941 
1942 
1943 


than with dates along the top and months down the side, giving twelve rows and five 
columns. When it enables space to be saved columns should be numbered or lettered 
and explanations added at the foot of the table. Rules should not be inserted unless 
special circumstances make them necessary. 

Tables are numbered serially and referred to in the text by number, table 1 (lower 
case t) etc., as the exact position of the table in the text is governed by its position on the 
printed page. 


(v) Mathematics 


It is not necessary to give detailed derivations of mathematical expressions and 
formulae in a published paper when the work is straightforward ; it is quite sufficient to 
‘ndicate the method of treatment and the final results. When interested readers are likely 
to find the argument difficult to follow fuller treatment is called for. 

Extraordinary care should be taken with mathematical scripts ; especially for ex- 
pressions in the least complicated, legible longhand is preferable to typing, as most 
typewriters have no characters suitable for indices and subscripts, for which our printers 
use smaller type; further, the letter “1” is often used for the figure one, which may lead 
to ambiguity. For simple fractions, the solidus / should be used instead of a horizontal 
line, e.g., A/27 rather than 7 care of course being taken to insert brackets where necessary 
to avoid ambiguity. Exceptions to this rule are the proper fractions 4, 4, $, etc., which are 
available as single types and so are better left in this form, though such fractions as 9/17 
12/25 are better written with the solidus. 

Full formulae or equations should be displayed, i.e. written on a separate line, and 
equalities numbered serially (1), (2), ... onthe right-hand side. In displayed expressions 
the solidus should be used with discretion. Short expressions occurring in the text may 
be incorporated with the text, as in the expression \/2m above. 

When complicated expressions occur repeatedly in mathematical work they should 
not be written in full every time. Both the reader and the compositor will be helped if 
they are represented by a special symbol of which the meaning is defined when the 
expression occurs first. 


Instructions to authors 3 


_Care should also be taken with exponential expressions and with subscripts ; the 
principle to bear in mind is simply that matter above or below the level of the normal type 
involves special treatment which leads both to less elegant work and to greater expense : 


P subscript » is simple, P,,, but a subscript ny (P»,) is not; e” is simple but e” is not, and 


may on occasion be written exp «?. (a?+4b2)! and \/(a?-+-62) are preferable to Va2_B?. 
As other fractional indices commonly occur the use of the fractional index is recommended. 

For the logarithm of x to base 10 write ‘log x’ and for that to base e write ‘In x’. 
Do not write ‘log,) x’ or ‘log, x’. 

An interesting and explanatory document on the setting out of mathematics is the 
London Mathematical Society’s ‘‘ Notes on the Preparation of Mathematical Papers”, 
obtainable from Messrs. C. F. Hodgson & Son Ltd., 2 Newton Street, Kingsway, Strand 
W.C.-2, price 1s, ; 


(vi) References 


In the text bibliographical references are made by giving the name of the author 
and the year of publication in brackets, e.g. (Jones 1942), and details are given in 
the last section, ‘“‘ REFERENCES”. Here the references are arranged in alphabetical 
order of authors’ names and in date order for each author ; the form of reference is: 
Author, initials, year, journal, volume, page, e.g. 


KOEHLER, J. S., 1941, Phys. Rev., 60, 397. 


The abbreviations used are those on the World List of Periodicals, with the slight 
difference that capital letters instead of lower case are used for adjectives. It will pro- 
bably be most convenient for authors to give all journal names in full ; the appropriate 
abbreviations will then be made editorially. Books are referred to as follows :— 


Love, A. E. H., 1944, Elasticity (Cambridge : The University Press). 


§ 3. DRAWINGS 


Drawings should be in Indian ink on tracing cloth, tracing paper or white card, with 
lettering in soft pencil. ‘The drawings should in general be sufficiently large to allow of 
reduction in printing and the lines should therefore be bold ; the frame lines of graphs 
should be slightly finer than those of the plotted curves. Essential photographs can be 
reproduced as plates, but they should be avoided when possible because of their expense. 

If authors are unable to submit drawings in the required form, arrangements can be 
made for these to be drawn and the cost price charged to the author. 

If the drawings are made on tracing cloth or paper, it is advisable to send a set of 
prints with the MS. so that the originals may not be handled unduly. The pencil 
lettering will show on the print if it is done in soft (2B or 3B) pencil with a flat “ chisel ” 


point. 
Figures should be numbered serially in arabic numerals and referred to in the text as 
figure 1, figure 2, ... i.e. word in full with lower case f. 


§ 4. ABBREVIATIONS AND JOURNAL STYLE 


For symbols indicating position or quality, use roman type, e.g. the point A; ha, 
value of & in air, kg, value of k in glass. In MS. these are left with no comment. 

For symbols indicating quantity use italic type ; e.g. current J, voltage V, angles 6, 
6,...0,. In MS. these have single underlining, e.g. V. 

For vectors use bold face letter (r). In MS. these are to be underlined with a wavy 
line, e.g. r. Italic bold face type (R) should be used for universal constants. 

Greek letters should be explained by a pencil note in the margin, e.g. ‘ Greek alpha’. 

Units are used in abbreviated form in the singular, e.g. cm. (not cms.), mm., km., 
kc/s., Mc/s., sec., wa., ma. 

Lower case “a” is used for unit of current to distinguish from Angstrom units ; 
most other units, often derived from names, are rendered in small capitals (indicated by 
double underlining) ; Angstrém units a., volts v., degrees Kelvin °x., degrees centi- 
grade °c., kilo-x-ray units kx. A short list of abbreviations is given in an appendix. 


4 Instructions to authors 


Care should be taken to see that spelling is consistent throughout each paper. The 
preferred style is -ize in such words as recognize, polarize. 


§ 5. LENGTH OF PAPER 


Papers should be kept as short as possible. This can often be achieved by attention 
to mathematics and diagrams. The time taken for publication of a paper increases with 
its length ; speedy publication can usually be assured for papers of less than 3000 words. 


§ 6. CONCLUSION 


It is hoped that authors will cooperate by seeing that these recommendations have 
been adhered to before communicating their papers ; such assistance would go far to 
remove the difficulties experienced in editing and in setting up type, which are the major 
problems in the publishing of a paper. 


APPENDIX 


ABBREVIATIONS AND SYMBOLS 


a. Amp. 

A. Angstrom unit. 

INGO Atomic unit. 

ALC. Alternating current. 
DCs Direct current. 
E.M.F. Electromotive force. 
Chine Ks Degrees Centigrade, Fahrenheit, Kelvin. 
V. Volt. 

X.U. X-ray unit. 

kva. Kilovolt-amp. 

ev. Electron volt. 

Mev. Mega electron volt. 
h, m, e, c, k Universal constants. 


= Is approximately equal to 
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The Structure of Photo-sensitive Lead Sulphide and Lead 
Selenide Deposits and the Effect of Sensitization by Oxygen* 


By H. WILMAN 
Applied Physical Chemistry Laboratory, Imperial College, London 


MS. received 2 Fune 1947 


ABSTRACT. An exploration has been made, by electron diffraction and photocon- 
ductivity experiments, of the structure of photosensitive lead sulphide and lead selenide 
deposits prepared by chemical deposition and by sublimation in vacuo or in oxygen. In 
spite of the large variations in photoconductive properties which the deposits show as a 
result of treatments such as baking in vacuum or oxygen, the PbS and PbSe crystals were 
found to have a lattice axial dimension which was constant to 0-1 % in all samples measured, 
even in strongly oxidized deposits in which there was estimated to be 10 to 20% of the 
oxidation product. ‘The oxidation product was identified from.its diffraction pattern as 
PbO.PbSO,. The photoconductivity measurements, however, indicate that the sensitivity 
is associated with deviations from the stoichiometric compositions PbS and PbSe, and that 
oxygen treatments increase the sensitivity of PbS to wavelengths between 1 and 3 », but do 
not do so for PbSe. 


Sle NE ROD WiC LON 


URING the war lead sulphide photoconductive cells were manufactured in 
Germany by Gudden, Kaspar, Kutzscher and others, but the methods of 
production have not been published. In 1944 methods of manufacture 
of PbS and PbSe cells were developed at the Admiralty Research Laboratory, 
Teddington (for details of the vacuum-sublimed deposits see Starkiewicz, 
Sosnowski and Simpson 1946, Lee and Parker 1946), and a brief account is 
given below of the structure of the deposits, as revealed by the use of electron- 
diffraction methods. Electrons of about 60kv. voltage were used with a camera 
length of approximately 50cm. in a Finch type electron-diffraction camera 
(Finch and Wilman 1937). Early in 1946 further electron-diffraction investi- 
gations provided definite information on the amount, chemical nature and form of 
the oxidation product obtained by heating or subliming (at temperatures above 
500°c.) the PbS layers in oxygen at about 0-1 mm. Hg pressure, during their 
sensitization as photoconductive or photovoltaic cells (Starkiewicz, Sosnowski and 
Simpson 1946). 
The crystal lattice structure and parameters of the PbS and PbSe were 
practically identical for all the deposits examined ; these are stated in $6. Of the 
photosensitivity data only a brief outline of some general characteristics is given 


in $6. 


* This work was carried out in the Applied Physical Chemistry Laboratory, at Imperial College, 
London, by arrangement with Professor G. I. Finch, M,B.E., D.Sc., F.R.S. 


PROC, PHYS. SOC. LX, 2 9 


118 H. Wilman 


The vacuum-sublimed deposits were all prepared in well degassed pyrex 
pumping systems (0-5 to 1 hour degassing at 500°c.) evacuated by mercury 
diffusion pumps to about 10~-* mm. Hg, with liquid-air trap to avoid mercury 
contamination. 


§2. THE STRUCTURE OF LEAD SULPHIDE AND LEAD SELENIDE 
DEPOSITS FORMED BY CHEMICAL PRECIPITATION 


Deposits of lead sulphide with good photoconductivity prepared by Dr. F. 
Kicinski at A.R.L. were precipitated as mirror-like layers on pyrex substrates from 
a solution of lead acetate and thiourea by adding sodium hydroxide solution, and 
after washing with weak ammonium sulphide solution were dried in air at 100°C. 
Figure 1 shows an example of the electron-diffraction patterns obtained from such 
mirror-like slightly translucent deposits, which are dark olive-green by transmitted 
light. The arc positions show that the PbS crystals tended to be orientated with 
cube faces, i.e. {001}, parallel to the polished pyrex substrate. Similar patterns, 
but with less distinct arc development, were given by thicker opaque deposits 
(~1 thick) in which the PbS crystals in the surface layers had only to a very 
slight extent the preferred {001} orientation. The mean crystal size was estimated 
to be at least 150 a. and was in some cases of the order of 500. 

Translucent olive-brown PbS deposits were also prepared for comparison with 
the above, by passing H,S gas over a lead-acetate solution at room temperature, 
slightly acidified with acetic acid. These floating mirror-like films were picked 
up on nickel gauze, transferred to the surface of an H,S water solution to wash 
them, removed, and a drop of Aquadag graphite added for lattice-constant 
reference purposes; they were then dried im vacuo in the diffraction camera. 
The ring-patterns obtained at normal transmission, figure 2, became weakly arced 
on inclining the film, showing that the PbS crystals had their cubic faces, {001}, 
more or less parallel to the surface of the solution on which they were formed. The 
mean crystal diameter was about 250 a. 

Similar mirror-like, translucent reddish-brown layers of PbSe were obtained 
by passing H,Se over lead-acetate solution acidified with acetic acid, and these gave 
patterns like figure 13, where the {110} graphite ring is also present for calculation 
of lattice dimensions. ‘The PbSe crystals had a mean diameter of about 250 a. and 
had tended to grow with {001} planes parallel to the surface of the solution. Such 
PbSe layers showed no change in diffraction pattern after heating for 30 minutes 
in air at 100° to 130°c. or in steam at 130° to 150° c, 


§3. THE STRUCTURE OF LEAD SULPHIDE AND LEAD SELENIDE 
DEPOSITS PREPARED BY VACUUM SUBLIMATION 
Some of the early experimental PbS cells made by the vacuum sublimation 
method were cut open and immediately re-evacuated in the electron-diffraction 
camera, and an electron beam allowed to graze the inner cylindrical surface of the 
deposit on the electrode region of the pyrex cell wall. ‘Three types of deposit were 
examined : (a) PbS volatilized and condensed zm vacuo then superficially acted upon 
by sulphur vapour at low pressure and temperature about 500°c. (figure 3); 
(b) PbS volatilized and condensed in presence of sulphur vapour (figure 4); 
(c), as for (6), but resublimed in oxygen at 1 mm. Hg pressure. ; 
In all these the mean crystal diameter of the PbS was estimated to be between 
150 and 500, or more, and there was usually a marked tendency for the PbS 
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crystals to grow in a one-degree orientation, either {111}, {110}, {210}, or {311}; 
in some cases two, or even three, of these were present in the. same layer. For 
example, figure 4 shows that {311} planes in some crystals and {210} in others were 
parallel to the pyrex substrate. "The type of orientation devetaped presumably 
depends on the temperature of the surface and the rate of condensation. The one 
deposit of the third type examined had a larger crystal size than the others, 500. 
or more, as shown by the incompletely formed spotty rings, cf. figures 5-8. This 
was to be expected, for the temperature of sublimation of the PbS was higher 
(~580° c.) in oxygen than in vacuo (450°-500° c.), and since that of the substrate 
was also higher, there was more rapid aggregation of the PbS to form larger crystals. 
Two exceedingly faint unidentified rings were present in this pattern in addition 
to the PbS rings. 

Seven vacuum-sublimed translucent orange-brown PbSe deposits in experi- 
mental cells were examined. The PbSe used was prepared by passing H,Se into 
lead-acetate solution followed by much repeated washing of the precipitate and 
drying at about 100° to 130°c. All these deposits yielded clear ring or arc patterns, 
as in figure 14, due to Pb Se crystals either randomly disposed or in strong one- 
degree {001} orientation on the substrate. The estimated mean crystal size 
varied between about 150. and 500. or more, increasing with increasing duration 
and temperature of the vacuum bakes which the layer had undergone after the 
initial condensation. In one case where a cell had developed a noticeable leak at 
one of the glass-tungsten seals during formation of the layer, which had become 
grey in colour, an unidentified diffraction pattern was obtained (figure 16, table 4). 
This presumably corresponds to one or more of the compounds of Pb, Se andO. 


Shih oObRUCLURE OF LEAD SULPHIDE DEPOSITS SENSITIZED 
BY SUBLIMING AND BAKING IN OXYGEN 


(1) Photoconductive deposits 


Two cells A and B (made by J. Starkiewicz at A.R.L.) which had undergone the 
normal type of activation treatments by slow sublimation and further baking in 
oxygen, yielded the electron diffraction patterns, figures 5 and 6; net-plane spacings 
are given in table 1. In both cells the pumping stem had a narrow constriction 
impeding escape of SO, gas to the liquid air trap, and the final (modulated-light) 
sensitivity reached the average good high level attained with most cells of this type. 
Both PbS deposits had a matt grey-black surface and were opaque. 

The patterns of spotty arcs (figures 5 and 6) showed that the PbS crystals had 
relatively large size, probably ~1000 a. diameter, and had grown so that some had 
{110} planes parallel to the pyrex substrate, others {001}. ‘The lattice dimension 


| was a=5-917kx. relative to graphite a=2-456,kx. he additional fainter arcs in 
_ the patterns showed the presence of an oxidation product (forming about 10-20%, 
_ of the diffracting layer) which was also crystalline, with mean crystal diameter of 
the order of 500 ., and was also (especially in cell A) strongly orientated, the net- 


plane parallel to the substrate having a spacing of 1-228kx. ‘T’his oxidation product 


_ was identified from its net-plane spacings as PbO.PbSO, (see (iii) below), the 


crystal habit apparently varying according to the conditions of oxidation and causing 


- some variation in ring intensities, without much change in the lattice dimensions 


y 
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and net-plane spacings. 
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(11) Photovoltaic deposits 


Similar patterns, figures 7 and 8, were also obtained from two cells, C and D, 
which had been activated as photovoltaic or barrier-layer cells by J. Starkiewicz at 
A.R.L., i.e. after sublimation in oxygen at 0-1 mm. Hg pressure, followed by baking 
in oxygen, a direct current of 2 ma. was passed through the cell during bakes at 250° 
to 400° c. (Starkiewicz, Sosnowski and Simpson 1946). Both had about average 
(modulated-light) photoconductivity. | During oxidizing activation of cell C the 
SO, was impeded from reaching the liquid-air trap by constrictions in the pump- 
ing stem, while in the case of D no impedance was applied and the SO, formed 
was quickly trapped. The conditions of oxidation were thus appreciably different 
in the two cases; nevertheless the diffraction patterns showed (table 1) that many 
of the net-plane spacings of the oxidation product present besides the PbS, were 
essentially the same within the accuracy of measurement of the radii of the rings or 
arcs which were often very faint.- Moreover, the spacings agreed closely with most 
of those from cells A and B, and well with many of those found for PbO.PbSO , 
especially from the pattern obtained by heating PbS in air to 350°c. (see (iii) 
below and table 2). The relative intensities and circumferential positions of the 
corresponding rings in the two patterns were different, and differed from those 
of cells A and B, showing that the PbO.PbSO, crystals were orientated differ- 
ently relative to the substrate. In cell C the {001} planes of the PbO. PbSO, 
tended to be parallel to the substrate (cf. tables 1 and 2), while the PbS 
crystals were mainly in {001} but some in {110} orientation; and in cell D some 
PbO.PbSO, crystals had {111} orientation (as in cells A and B), and in others 
{102} planes were parallel to the substrate, the PbS being again strongly orientated 
partly with {001} and partly with {110} planes parallel to the substrate. 


(111) Identification of the oxidation product as PbO.PbSO, 


‘The net-plane spacings and relative diffraction intensities of the oxidation 
product from table 1 were compared with the available x-ray data for iead com- 
pounds, including that of Hanawalt, Rinn and Frevel (1938) for Pb, PbO (yellow 
rhombic form), PbO,, Pb;04, PbSO, and PbS,O, (lead thiosulphate), PbCO,, 
2PbCO,.Pb(OH), and Pb(NOs),, and with electron-diffraction data for the last 
two of these compounds previously obtained in this laboratory. 

Other recent published work on the Pb—-O and Pb-—O-S compounds was also 
compared with the data of table 1, but the following gave no agreement: PbO 
(yellow rhombic)—Halla and Pawlek (1927), Rencker and Bassiére (1936), 
Petersen (1941), Bystrém (1943); PbO (red, tetragonal)—Dickinson and Friauf 
(1924), Darbyshire (1932), Clark and Tyler (1939), Clark and Rowan (1941), 
Moore and Pauling (1941); Pb,;0,—Clark, Schieltz and Quirke (1937), Gross 
(1941, 1943), Straumanis (1942), Bystrom and Westgren (1943); Pb,O;— 
Clark, Schieltz and Quirke (1937), Baroni (1938), Gross (1941); Pb;Q,—Clark, 
Schieltz and Quirke (1937), Clark and Rowan (1941), Davidson (1941); PbO,—- 
van Arkel (1925), Darbyshire (1932); PbO,—Clark and Rowan (1941); 
5PbO.2H,O—Clark and Tyler (1939); PbSO,—Basche and Mark (1924), 
James and Wood (1925); 2PbO.PbSO,, 3PbO.PbSO,, 4PbO.PbSO,—Clark, 
Mrgudich and Schieltz (1936). 

A useful series of x-ray patterns of many of the above is given by Clarke (1940), 
and x-ray data for the mineral “lanarkite”, PbO.PbSO,, by Richmond and 
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Wolfe (1938). The latter authors give the size of the monoclinic pseudohexagonal 
unit cell as a= 13-73 a.,b=5-68a., c=7-074., and B=116° 13’. Thecell is base- 
centred, thus only AAl diffractions appear which have (h+k) even. A calculation 
of the net-plane spacings corresponding to possible diffractions gave the values in 
table 2, which includes all values above 2 a. and many below. Relative intensities 
were obtained by the experiments described in §5, below. 

In the comparison of the lists of spacings and diffraction intensities in tables 1 
and 2 the following points must be realized: (i) owing to the orientation of the 
crystals in the photocell deposits, not all the diffractions in table 2 will be visible 
above the shadow edge in the reflection patterns; (ii) the relative intensities of the 
arcs which do appear will be different from those in table 2, owing to the orientation 
and to absorption which weakens beams emerging near the shadow edge; (iii) 
planes which are parallel to the surface in most of the crystals give especially strong 
arcs whose centres lie on the plane of incidence ; (iv) there may bea slight difference 
in the lattice dimensions and the relative diffraction intensities such as that observed 
in the patterns from PbS heated in air at 350° c. and 550° c. (table 2), although the 
net-plane spacings are all standardized by comparison with the graphite 110 
spacing; (v) strong PbS diffraction arcs may obscure other weaker arcs in their 
vicinity ; (vi) smallness of the radii and of refraction effects decrease the accuracy 
of measurement of the larger net-plane spacings. 

Bearing in mind the above conditions the oxidation product in the photocell 
deposits was identified from the spacings given in table 1 and from the circum- 
ferential positions of the arcs in the patterns, figures 5 to 8, as PbO . PbSO, (or 
approximating closely to that composition and lattice structure) for the following 
reasons : 

(i) The spacings in the first four lists in table 1 agree well with themselves and 
with many of those in table 2, notwithstanding the faintness of some of the diffrac- 
tions. (ii) The lanarkite structure accounts for the large spacings of the two 
innermost diffractions of figure 8 (cell A). (iii) On the basis of the lanarkite 
structure, the net-planes which are orientated parallel to the substrate are densely 
populated planes having small indices, {111} in cells A, B and D, and {111} and 
{102} together in cell D; the second and third orders of reflection from the {111} 
planes give as the {111} spacing the value 3-68, A. in agreement with the value for 
lanarkite in table 2. (iv) The circumferential positions of the arcs on the other 
rings appeared to correspond well with those to be expected from lanarkite 
crystals in these orientations. (v) The plane spacings of the transmission patterns 
obtained after complete oxidation of PbS films in air at 350°c. and 550° c. agreed 
closely with those calculated for lanarkite (table 2). 

In order to obtain diffraction patterns from other Pb-O-S compounds 
several other materials were examined by electron diffraction, and though they did 
not correspond to that of the oxidation product discussed here, they are of some 
interest and provide another example of appreciable changes of diffraction 
intensities associated with only slight changes in lattice dimensions. Figure 11 
and table 3 (spacings calibrated by graphite) show the diffraction data obtained by 
transmission from a floating white surface layer precipitated by passing SO, gas 
over lead-acetate solution. ‘This is certainly not PbSO,, but the pattern seems to 
have some similarities with those of the 2-, 3- and 4-PbO.PbSO, oxysulphates. 
When this material was heated in air for 5 minutes at 300°, it gave a pattern (see 
table 3), in ywhich the ring intensities were appreciably changed though their 
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relative radii were nearly unaltered; after a further 45 minutes in air at 550°C. it 
gave figure 12 (table 3), which may possibly be closely related to the x-ray pattern 
from 2PbO.PbSO, though it is not identical with this. 


‘Table 3. (a) and (d) are nearly identical and are closely related to the pattern 
(b) of table 2; the spacings of (c) correspond with most of those of 
table 2(c) but the intensities are different. 


(a) ; (b) (c) as 
White pp. from Pp. of (a) heated 5 min. Spec. of (b) heated 45 min. 
Pb acetate+-SO,. at 300° c. in air : at 550° c in air 

dix) Intensity d (kx.)* Intensity d (kx.)* Intensity 
8:07 VVI — as == = 
6°75 VI 6-4 VF = = 
5-182 F 5-14 VE 4-45 F(D) 
4-261 M 4-24 ME 4-11 F(D) 
S758 MS SI07/S: MS 3-68 MF(D) 
3-414 MF 3-40 MF 3-48 F(D) 
3-218 MS 3-215 M 3-286 MF 
3-122 MS 3-110 M 2-963 Vs 
2°905 MS 2:905 S 2:849 S 
D123 S PINS) S (2:63) VF 
2°576 F Dos 7/5 VF 2-528 M 
2-469 MF 2-444 F 2-395 MF 
2:341 VF — — — == 
2-248 M 2-246 F 2-263 F 
2:136F ME 2:127t MF 2-182 ¥F 
2-058 MS 2-064 M 2:047 M 
2025 MS 2-025t M — a 
1-958 MF 1-960 F 1-954 F(D) 
1-903 MF 1:901 F 1-896 F(D) 
1-865 F 1-852 F 1-843 F(D) 
1-831 F 1-825 VF — a 
1-760 VE — = 1:754 F(D) 
1-719 MF i 7h MF 1-716 F(D) 
1-704 F — — = — 
1:659 FE 1-657 F 1-669 F(D) 
1-644 F — = — = 
1-606 MF 1-613 r 1-615 F(D) 
1:571 MF 1-566 F — a 
1:502 F 1-502 VF 1-520 F(D) 
1-457 M 1-454 MF 1-462 M(D) 
1-446 M = — 1-419 VF 
1-393 MF 1-391 VF 1-386 VE 
1-368 F — — 12352) VF 
1-339 F — — = es 
1-318 VF — — ASA MF(D) 
1-305 VF a= = aes i= 
1-283 VE — — 1-289 MF(D) 
1-254 M e258 MF 1-249 VF(D) 
1-205 VF — — 1-206 VE(D) 
1-189 VE — = 1-169 VF(D) 


* Relative to graphite, d)4)>=1:228 kx. 
+ May be entirely the graphite 100 ring. 
} May be partly but only in small part due to graphite 101 diffraction. 


§5. THE OXIDATION OF PbS AND PbSe IN AIR AT ATMOSPHERIC 
PRESSURE 


Translucent olive-green films of PbS with strong metallic reflexion were 
formed by passing H,S over lead nitrate solution. After washing by floating on 
distilled water these gave normal PbS diffraction patterns like figure 2, When 
such layers were heated in air at atmospheric pressure at 350-370° c. for 10 minutes 
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Figure 1. Chemically deposited PbS layer Figure 2. PbS layer from H.S over Pb 
(good sensitivity). - acetate solution. 7 


Figure 3. PbS cell; layer formed by Figure 4. PbS cell; Jayer formed by 
vacuum sublimation and treatment with S. vacuum sublimation in presence of 5S. 


Figure 5. PbS cell“ A”; layer sublimed Figure 6. PbS cell“ B”’ ; layer sublimed 
in Oy. in Oy. 


Figure 7. PbS cell ““C” ; formed in Figure 8. PbS cell “D” ; formed in 


? . 
oxygen, but photovoltaic type. O,, but photovoltaic type. 
To face page 124 
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Figure 9. PbO.PbSO, by heating PbS in Figure 10, PbO.PbSO, by heating PbS 
air at 350° c. for 10 minutes. (Two in air at 550°c. for 5 minutes. 
outermost rings graphite.) (Graphite present as reference.) 


Figure 11. White layer by SO, over Pb Figure 12. Deposit of figure 11 heated in 
acetate solution. air at 550° c. for 45 minutes. 


Figure 13. PbSe layer from H,Se over Pb Figure 14. PbSe cell ; layer formed by 
acetate solution. vacuum sublimation. 


Figure 15. PbSe heated in air at 350° c. Figure 16. PbSe cell after a leak had 
for 10 minutes. caused oxidation. 
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they became transparent yellowish-grey, giving a different pattern, figure 9, which 
corresponds exactly to PbO.PbSO,. ‘Table 2 shows the plane spacings (cali- 
brated by comparison with graphite) and the relative intensities of the diffraction 
rings, and these are compared there with the data for lanarkite, PbO.Pb SOR 

A fresh PbS layer, heated in air up to 550° c. for 7 minutes and kept there for 
5 minutes, became milky-transparent and gave a pattern, figure 10 and table 2, with 
spacings nearly identical with those of figure 9, but with a difference in the relative 
intensities of the corresponding rings, the random disposition of the crystals in both 
cases was checked by photographs with the specimen inclined to detect any arcing 
of the rings. This slight change in lattice dimensions may indicate a slightly 
different degree of oxidation in the two cases with only a small difference in 
crystal lattice dimensions. The change in relative intensities of the rings may 
be attributed to a difference in the external shape of the crystals, which would 
cause different absorption of the various diffracted beams passing through 
different thicknesses of crystal in different directions (BOhm and Gantner 1928, 
Fordham 1940, Brandenberger 1945). 

Similar PbSe layers formed by passing H,Se gas over lead-acetate solution 
acidified with acetic acid gave unchanged ring patterns by transmission electron 
diffraction after heating the layers at 100° to 130°c. in air for 30 minutes or in 
steam at 130°to 150° c. for 30 minutes; but heating in air at 350° c. for 10 minutes 
converted them completely to a transparent yellowish oxidation product (figure 
15), the diffraction data for which are given in table 4. Comparison with the 
pattern data from the PbSe layer oxidized by a leak in a cell during a vacuum 
bake (see §3) seems to show some similarities in spacings though differences in 
ring intensities. 


SO. SUMMARY AND DISCUSSION OF RESULTS 
A. Electron-diffraction results 


The following results were obtained by electron diffraction: 


1. All the lead sulphide and lead selenide deposits (~500-10000 4. thick) 
which were examined by electron diffraction were found to be crystalline with no 
noticeable signs of any amorphous material; diffraction methods would only 
detect proportions of the order of 1° or more of the material. 

2. All the deposits had the face-centred-cubic structure of sodium chloride 
type, previously found for PbS and PbSe by x-ray diffraction (Strukturbericht, 
Z. Kristallogr., 1, 1928). 

3. The mean crystal size and the extent and nature of the orientation present 
in these polycrystalline deposits varied widely from one specimen to another. 

4. All the PbS deposits showed an absence of appreciable variation of the 
lattice dimensions within the crystals of each deposit, and all the deposits examined 
had the same lattice dimensions within the limits of accuracy of the measurement 
(to approximately 0-1°/,), whether the deposits had been prepared by (a) chemical 
deposition from a lead acetate + thiourea + sodium hydroxide reaction on to pyrex 
substrates, (6) as a floating mirror-like layer by passing H,S gas over concentrated 
slightly acidified lead acetate or lead nitrate solution, the film being washed by 
floating on an H,S-water solution and dried in vacuum at room temperature, 
(c) sublimation of well-washed and dried precipitated PbS, in an evacuated pyrex 
cell at about 10-® mm. Hg pressure, with or without further short bakes in vacuum 


126. H. Wilman 


Table 4. 
(a) (bd) ; (c) 
PbSe PbSe heated in air PbSe cell which developed 
(= OvsM/kesen) i 1 hour at 350° c. a leak during a vacuum bake 
(figure 15) (figure 16) 
d (kx.) Intensity hkl d (kx.)* Intensity d (kx.)f Intensity 
6-16 MF == = 
5-34 F 4-67 F 
4-45 F(D) 4:33 F 
3:92 MF 4-06 F 
3-68 MF 3-69 a 
3e532 M 111 3-55 M 3-413 
max pe 3.068 MIF 
‘058 M 200 — — : i 
: 2-970 VIS(B) 2:988 VE 
2°845 F 2838 MF 
2-763 MS 7 ia S 
2°578 MF = = 
2°506 F 2-468 F 
2-404 F — a 
2:291 MF 2:309 see 
2°163 S 220 Dest M 27191 
5063" M 2-110 MF 
1-992 MF 2-009 F 
1-946 MF 1-924 F 
1-906 F — = 
1-844 MF 311 1eSSr7 MF 1:841 F 
1-766 MS 222 1-764 MF 1-764 MS 
1-702 MF 1-706 F 
1-678 MF 1-667 F 
1-634 M 1-627 F 
1-588 F 1-570 M 
Hey 28) VE 400 1-529 MF 1-540 MF 
— = (1-520) VF 
1473 M 1-465 F 
1-403 MF 331 1-424 F 1-425 F 
1-368 iS) 420 ey VF A352 MS 
1-304 MF 1-296 MS 
1-287 MF 1-269 VE 
1-248 S 422 1-254 F 1-238 M 
1-177 F on 3 a 1-172 M 
1-149 VF 
1:127 VE 
1-097 VF 
1-076 M 
1-023 M 


* Relative to graphite d,,,3=1-228 kx. 
t Not relative to graphite ; absolute values as a whole may be a few % in error. 
{ May be entirely due to the graphite 100 diffraction ring. 


at between 200° and 400°c., (d) sublimation in oxygen at 0-1 to 0-3 mm. Hg 
pressure at 500° to 580° c. 

A similar constancy (to ~0:1%) was also found for PbSe deposits prepared 
(a) as a floating mirror-like layer by passing H,Se gas (from either FeSe + dilute 
HCl, or Al,Se,+ boiled distilled water) over concentrated lead acetate solution 
acidifiedwith acetic acid, (b) bymethod (a) followed byheating in air at 130°to150°c. 
for 30 minutes, (c) by method (a) followed by heating in steam at 150°c. for 
30 minutes, (d) by sublimation on to pyrex in a cell evacuated to 10-8 mm. Hg, 


with or without further short bakes (of the order of 10-30 minutes) in vacuum at 
200° to 400° c. 


A te a 
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‘The mean lattice constants found for PbS and PbSe were a= 5-917 + 0-005 kexts 
and 6-117 + 0-005 kx. respectively*, relative to a= 2-456, kx. for graphite (Trzebia- 
towski 1937, Nelson and Riley 1945, Finch and Fordham. 1936, Finch and 
Wilman 1936, cf. also Wilman 1940), the 110 graphite diffraction ring with 
d=1-228kx. being used as the comparison standard. 

Previous X-ray measurements, mainly for the more or less impure minerals 
galena (PbS) and clausthalite (PbSe) have given the following lattice dimensions: 
for PbS, a=5-96,..,4. (Lehmann 1924), 5-914. (Kolderup 1924-5), for very 
pure galena from Pribram Bohemia, 5.93. (Ramsdell 1925), 5-954. (Gold- 
schmidt 1927-28), 5-964. (Bravo 1926), 5-91+0-03a. (Vlasak and Trousil 
1934) 5-935. (von Zeipfel 1935 a, corrected for more accurate value of A of 
Al-Ke x-rays by S.von Friesen 1935), 5-92334kx. (von Zeipfel 1935 b), 5-94 a. 
(Hanawalt, Rinn and Frevel 1938), 5-92 a. and 5-93 a. for galena containing 6°, Bi 
(Goldschmidt 1927-28), for PbSe, a=6:14a. (Ramsdell 1925), 6-162. (von 
Olshausen 1925); and 6-13;,,a. from Pb>Se melted together in N, (Gold- 
schmidt 1927-28). 

The constancy of the lattice axial dimension to within 0-1°% for PbS and PbSe 
in the various degrees of oxidation can be taken as showing that the proportion of 
excess Pb atoms or 5S defect (“‘holes’’), or of S or Se respectively replaced by O ions 
in the PbS or PbSe lattice is very small, and probably less than about 1%, even 
in the strongly oxidized deposits in which up to 10 to 20% of PbO. PbSO, has 
segregated out. For example, if the ionic radii of Pb++, S~~ and O ~ are 1-32, 
1-74 and 1-32 a. respectively, then replacement of 1°% of the S~~ by O-~ would be 
expected to result in a reduction in the lattice dimensions of the PbS by approxi- 
mately 0-:08°%. However, presence of excess Pb or of “holes”? due to missing 
S~~ ions might tend to increase the lattice dimension and counterbalance a reduc- 
tion due to presence of O- “ions. Since the outer diffraction rings are not notice- 
ably less sharp than the inner rings, however, it seems that lattice dimension 
variations larger than 0-1.°4 do not occur throughout the PbS crystals of 500 to 
1000 a. diameter, or occur only locally to a very small extent. 

5. The fact that well-defined diffraction rings were obtained from the photo- 
sensitive deposits by the “reflection”? method, shows that these deposits were 
relatively rough even when they appeared to be highly reflecting and optically flat, 
the height of the “‘ peaks and valleys”’ being at least 150.4., and mostly 500. or 
more in the vacuum-sublimed deposits. Correspondingly the electron-diffraction 
patterns show the structure of the deposit to this depth below the outer parts of the 
surface. 

6. While the electron-diffraction patterns from the chemically-prepared, and 
also the vacuum-sublimed, PbS and PbSe deposits showed no trace of any other 
material present, those from PbS sublimed or baked in oxygen at ~0-1 mm.Hg 
pressure (the SO, formed being trapped by liquid air cooling) showed besides the 
strong PbS rings or arcs a much fainter pattern due to an oxidation product, which 
was estimated to be present in the surface layers to the extent of about 10-20%. 
When the oxidation was slight this pattern was represented by only one or two very 
faint diffuse rings, but in the four typical cells of this type examined, two of photo- 
conductive and two of photovoltaic type, it consisted mainly of clearly defined rings 

* This can be converted to Angstroms by multiplying by 1-00202, the factor recently agreed on 
by the X-ray Analysis Group of the Institute of Physics and the American Society for X-ray and 


Electron Diffraction. 
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or sharp arcs, from the positions of which the oxidation product was identified as 
lead oxysulphate (lanarkite), PbO.PbSO,. Translucent PbS layers heated in 
atmospheric air at 350°c. for 10 minutes were completely converted to a pale 
yellow transparent material also identified as PbO.PbSO,; and heating at 550 C. 
in air, followed by cooling to room temperature, caused only very slight changes in 
the lattice dimensions but distinct changes in the relative intensities of the diffrac- 
tion rings, probably due to a different crystal shape having been developed, causing 
greater absorption for some diffracted beams than for others. 

Lead selenide deposits, when similarly heated or sublimed in oxygen, were also 
converted into a greyish-yellow oxidation product; this was not identified owing to 
lack of diffraction data for the lead oxyselenium compounds. No improvement in 
photosensitivity was obtained from these oxidizing treatments of PbSe. 


B. Photoconductivity results 


With regard to the photosensitivity of the PbS and PbSe cells, only the 
following results will be stated here for consideration with the diffraction results, 
and a more detailed account of PbS cells will be published shortly by Starkiewicz 
et al. 

1. The photoconductivity properties of the deposits varied widely according 
to the conditions of formation of the deposit (e.g. in vacuum-prepared cells the 
temperature, rate of condensation, thickness of deposit, pressure of residual gas in 
the apparatus, nature of substrate etc.), and the further heat treatments in vacuum 
or in oxygen. No definite relationship was observed between the magnitude or 
the spectral distribution of the photosensitivity and the crystal size or kind and 
degree of orientation of the PbS*or PbSe as shown by the diffraction electron 
patterns. Since no full and systematic investigation has yet been made, only brief 
general indications will be stated here about the photoconductive and other 
properties of the deposits. It was noticeable, however, that the sensitization of 
PbS cells by subliming the material (or baking it) in oxygen led to a large increase 
in the crystal size as shown by the incomplete formation (spottiness) of the sharp 
diffraction rings due to the small number of crystals in the path of the electron 
beam. It is well known that crystal growth is much faster in the presence of 
oxygen than in vacuum, though it is not clear whether such accelerated growth is 
accompanied by appreciable diffusion of oxygen into the crystal lattice in solid 
solution. In the present case the sensitization in oxygen was found (Lee and 
Parker 1946) to be associated with development of a peak in the spectral sensi- 
tivity curve at a wavelength near 2°7 in addition to that at 1-1, evidently indi- 
cating strong adsorption or absorption of oxygen. The relation of the rise in infra- 
red sensitivity to the observed increase in crystal size and the formation of rela- 
tively large crystals of PbO.PbSO, is not yet clear, but the state of over-oxida- 
tion evidently corresponds to segregation of a large proportion of PbO.PbSO, 
among the PbS crystals. 

The observed absence of beneficial effects of oxygen treatments on the photo- 
sensitivity of PbSe deposits is probably due to the absence of appreciable 
“reducing” action, 1.e. liberation of lead, and to too rapid segregation of the oxida- 
tion product. Segregation is accelerated in this case because selenium is not 
readily removable as SeQ,, which is a solid at temperatures up to several hundred 
degrees Centigrade, whereas SO, is gaseous and is removed from the layer. 
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2. Asa brief general indication of the electrical properties of the deposits it may 
be stated here that the resistances of the PbS and PbSe deposits between graphite or 
platinum electrodes 1 cm. long and 1 mm. apart were in the range 10 to 1000 kQ 
and were approximately doubled when the cells were cooled in liquid air. The 
PbS cells, when illuminated by a broad, roughly collimated beam of unmodulated 
white light from a 36-watt lamp, showed a decrease of resistance of up to about 
half the dark resistance and a ratio of signal to noise of unity was observed 
with the cells exposed to 800 cycles/sec. modulated radiation from a 3 mm. 
aperture in front of a black body at 200°c., at a distance of the order of 70 cm. 
Further details of the sensitivity are given by Lee and Parker (1946). The 
PbSe cells showed only negligible photoconductivity (<2°% change in resist- 
ance) with continuous illumination, and with the modulated radiation had 
mainly a sensitivity of the order of 1-2°/, of that of good PbS cells. 

3. In spite of the observed degree of constancy of the lattice dimensions 
calculated from the electron-diffraction patterns ($4) obtained from cells whose 
properties differed widely, the electrical properties of both the PbS and PbSe 
deposits, measured at room temperature, varied within a wide range when 
previously well-outgassed cells were baked (at between 200° and 500° c.) while 
evacuated to about 10-6mm. Hg pressure. In the case of PbSe, for example, the 
translucent reddish-orange deposits made by vacuum sublimation in presence of 
excess selenium, had low resistance (~10 000 Q), very small “noise” and also very 
low signal and signal/noise ratios, when using asmall red-hot tungsten filament to 
give a light beam which was interrupted 800 times/sec. by a rotating slotted disc, the 
resulting alternating current being amplified. Successive periods of 10 minutes 
heating in vacuum at about 350° c. (Se begins to vaporize appreciably in vacuum at 
about 200° c.) then caused progressive increase of the resistance, noise, signal and 
signal/noise ratio, as in figure 17, until a maximum was reached for all these pro- 
perties nearly simultaneously. By driving in excess selenium and then again 
heating 7m vacuo in successive stages, these properties again passed through a 
maximum and returned to a nearly constant low level again (figure 17). (Note that 
in figure 17 the periods of baking stated were at the temperatures indicated by the 
apexes of the triangles in the lower diagram, additional to the 4-5 minute heating- 
up.) These observations, together with the visible loss of selenium from the cell 
to form a condensed deposit in the cooler parts of the pumping stem, appear to 
indicate that the vacuum baking resulted in at least the upper layers of the PbSe 
deposit losing selenium and passing from a composition having excess Se, through 
the stoichiometric composition PbSe with no Pb or Se impurity centres which 
corresponds to the maximum of resistance etc., to one with a defect of Se (ie: 
excess Pb). In view of these very wide variations the constancy of the lattice 
dimensions is remarkable. Wide variations in electrical properties of PbS layers 
have been described by Hinterburger (1942) and of PbSe layers by Eckart and 
Raithel (1941). The progressive decrease of the successive maxima in figure 17 
may possibly be due to gradual oxidation resulting from the minute trace of residual 
air in the vacuum system or oxygen evolved from the selenium side tube, or it may 
be due to some other cause such as loss of selenium from the lower parts of the 
layer so as to make it tend towards stoichiometric composition (and later recede 
from it again due to Se defect ?), or perhaps a gradual aggregation to form larger 
more isolated PbSe crystals owing to the energetic diffusion of the selenium in the 
layer. Figure 18, plotted from these same measurements shows a progressive 
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increase in resistance of the deposit for a given signal/noise ratio. The initial line 
5-6-7-8 passes near the origin, corresponding toa relation of form (S/N) = Ae™, 
Similar conditions existed for vacuum heat-treatment of PbS cells. Heat 
treatment of PbS or PbSe in oxygen resulted in more complex variations which 
will not be described here. 


C. General summary 


To sum up all the above results, the preliminary electron-diffraction obser- 
vations have provided a very much more definite knowledge of the nature of the 
photosensitive PbS and PbSe layers than has hitherto been available. It has 
identified the composition and defined the form of the crystalline oxidation 
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product formed during sensitization of PbS layers in oxygen, and in conjunction 
with photoconductivity measurements it has indicated that the crystal orientation 
has little effect on the photosensitivity compared with other factors such as 
especially, very slight deviations from stoichiometric composition, the size of Bic 
crystals in the layer, and probably the distribution of the PbO. PbSO among the 
PbS crystals. The axial length of the PbS crystals was a=5-917 ie and that of 
Pb Se was 6-117kx., and no deviation greater than 0-1°% occurred whether the PbS 
contained excess S or excess Pb (as indicated by the photoconductive 


<a | roperties), 
or whether it was unoxidized or relatively heavily oxidized, Mek ' . 
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e Figure 18. 
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ABSTRACT. Three dimensional vector loci are used to express the properties of 
resonant circuits and of the tuned transformer. An indication is given of the application 
of the method to other problems, including a simple low frequency selector circuit. 


S41. INEROD UEC TION 

HE study of the steady state behaviour of systems possessing stiffness, 

inertia and resistance forms a large part of any treatment of alternating 

current phenomena, of acoustics and of vibrating systems in general. 
Tt is not proposed to add further detail to this type of study, but rather to discuss 
a method of representing the behaviour of such systems which has been found 
attractive to students, which presents the solutions in a form which appeals to 
the imagination, and which is thought to afford a convenient visual summary 
of the phenomena concerned. 


§2. SIMPLE RESONANT SYSTEMS 


Consider a system whose behaviour is represented by the harmonic equation 
(figure 1) 


d?g doa 
L—- —+=.g=V. 
fp Rae ee tO 


Solid diagrams illustrating resonance phenomena ee 


where gq is the charge on the condenser and V the applied electromotive force. 


(0) B— 


Figure 1. Figure 2. Vector locus of admittance (OP) of 
circuit of figure 1. | Syrnbols in brackets 
denote alternative quantitics. For the 
parallel circuit OP represents V/Vr. 
The steady state solution may be written in terms of the impedance vectorZ 
or of the admittance vector Y. Let these be divided into real and imaginary 
components; so that 


Z=R+jX 
and -  Y=G4jB , 
“(Rjxyimeexy fo (2) 


For variations of X, the locus of Z is a straight line in the complex plane and 
that of Y is a circle of radius 1/2R (Hague 1945) (figure 2). 

The latter result depends only on the relation (2) and not on the means whereby 
the variation of X is effected. If now a three-dimensional diagram be plotted, 
with the variable parameter responsible for changes in X along an axis at right 
angles to the G and B axes, a model of the type shown in figure 3 is obtained. 
It will be noted that the form of the curve is essentially a single circular turn 
on a straight line. 

In this model the third variable is w/w, where w is the angular frequency 
of V and w, the resonant value, 1/\/(LC). By suitable illumination the three 
projections are obtained, showing respectively the variation of B and of G with 
frequency and the circle diagram, the latter, of course, being the vector locus. 
of the admittance. 

When used to illustrate current resonance in a series circuit, it is noted that 
the current J is related to the current /,, at resonance, by 


IL, =(1 —jX/R)/(1 +X?/R*), 


so that the circle may be made of unit diameter. Here the chord from the point O 
(figure 2) to a given point on the circumference represents the vector value of the 
current, in terms of the current at resonance, for the frequency chosen. It is. 
immediately apparent from the model that (i) the variations of conductance and 
of susceptance through resonance are simply two aspects of one process, (i1) the 
total variation of the susceptance is equal to the total variation of the conductance, 
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(iii) the maxima and minima of the susceptance curve occur on the steep parts 
of the conductance curve, (iv) the phase change in the neighbourhood of resonance 
is very rapid, (v) the phase at resonance is zero and (vi) the current is a small 
quantity in quadrature with the voltage at points well removed from the resonance 
peak. 

$35 EE eieN ED TRANSFORMER 


Let the secondary circuit be of resistance R, and reactance X, and let it be 
coupled to the primary circuit by mutual inductance M. The increments of 
primary resistance AR and reactance AX, caused by the presence of the secondary 


are given by 


bes Ase 
AR = 0M? py cover (3) 
xX 
AX =o 2 ee (4) 
R24 X¢ 


showing an exact correspondence in form with the real and imaginary com- 
ponents of Y (equation (2)). 

The increase in impedance of the primary circuit can thus be represented 
by the same kind of solid diagram; in fact, the projections of the solid curve are 
perhaps most familiar as representations of the transferred impedance. Here 
the radius of the circle is w2M2/2R,, so that the representation is exact for tuning 
the circuit by varying C, or L,, but if the applied frequency causes the variation, 
the diagram is only accurate for sharply tuned circuits. 


§4. MORE GENERAL METHODS 


For purposes of calculation and for more general application (including 
parallel resonance when the tuning is reasonably sharp) the solution may be 
expressed in terms of Q and of the ratio of the impressed frequency w to the 
resonant frequency w, (where w, is defined as the frequency for zero phase) 
(Lee and French 1943). 

For the series circuit, since Q=w,L/R and w,= 1/./(LC), we have 
R=1/(Qw,C). Substituting these values in equation (1) leads to 

CV /q=(1 —o?/w,”) + (9/Q)(o/or), 
where V and q represent the vector values of the applied voltage and of the charge 
on the condenser respectively. q, the charge at the resonant frequency, is given 
by 
Gr/= w/o, =) 84 a w?/c,?) 
1 
and WH S'S NN Ge 
YN TQe,jo— aa) 

In the complex plane (5) represents a unit circle in which tan ¢ = O(w,/w —w/w,) 
and, since the angle at the centre is 2¢ the solid curve may be plotted by means 
of a piece of graph paper wrapped round a cylinder. This procedure, followed 
by shaping a piece of wire to fit the curve, was used in making the solid models 


illustrated. 
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The method has been found valuable for illustrating the properties of tuned 
transmission lines and of dielectrics containing charged particles, elastically 
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Figure 3. 
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Figure 1. Rubber sheet stretched by factors of 2°5 and 2°0 in two perpendicular directions. 
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Tigure 1 (a2). Form of unstrained sample. 


Solid diagrams illustrating resonance phenomena 135 


bound or aperiodic. In addition, a simple low frequency selector circuit has 
been developed. This, in effect, is an interstage coupling consisting of two 
transformers, one tuned slightly above and the other slightly below the required 
frequency, with the output voltages in opposition. It will be realized that, 
whereas the output of each transformer separately is represented by a line drawn 
from the origin to a point on the circumference of the circle (figures 2 and 3), 
the combined output is the chord joining the ends of two such lines. Thus the 
maximum output can be the same as that of either component separately, but it 
drops much more rapidly towards zero on either side of the peak. Satisfactory 
discrimination against the second harmonic is obtainable at 9 cycles per second. 
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ABSTRACT. A method is described by which a sheet of rubber may be subjected 
to the most general type of homogeneous strain, with simultaneous measurement of the 
principal stresses and the birefringence. As predicted by the molecular network theory, 
the birefringence is found to be proportional (a) to the difference of the squares of the 
principal extension ratios, and (b) to the difference of the principal stresses in the sheet. 

An examination of the stresses shows that these cannot be accurately represented by 
the equations derived from the network theory. A closer approximation is given by a 
more general theory of Mooney, the simplest form of which gives a satisfactory account 
of the behaviour of a swollen rubber. For a dry rubber, on the other hand, the Mooney 
formula, though better than the network theory, does not adequately represent the 


experimental data. 


SU NPR OD Cal ON 


N a previous paper (Treloar 1947 a) the author has dealt with the theory 
of the optical properties of strained rubber on the basis of a model in which the 
rubber was envisaged as a network of long-chain molecules of randomly-jointed 

links. These links were assumed to be optically anisotropic, having two different 
optical polarizabilities, in directions parallel and perpendicular to their length. 


10-2 
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The theory was applied to the case of the most general homogeneous deformation 
of rubber, and it was shown that, for each of the three principal directions of strain, 
there is an associated optical polarizability, by which the refractive index for light 
having the corresponding direction of electric vector is determined. If Aj, A, and 
A, are the principal extension ratios, and 7, my and n, the corresponding refractive 
indices, the double refraction for light propagated along the direction A3 may be 
represented by the equation 


+2) 2 
Doh = ee FEN (cy 9) 0738Q— AB), vee (1) 


in which 7 is the mean refractive index, N the number of molecular chains per c.c. 
and «, — #» the difference of polarizabilities for the link of the chain. The quantity 
v, refers to the volume fraction of rubber in the case when the rubber is swollen in 
an optically inert solvent. Comparison of this result with the corresponding 
equation for the difference of principal stresses (Treloar 1947 a), 1.e. 


tj = NRTOROE= M&M), eee (2) 


where ¢,, ts, t; are the principal stresses, & is Boltzmann’s constant and T' the 
absolute temperature, leads to a linear relation between the birefringence and the 
stress-difference, which has already been examined experimentally for the case of a 
simple elongation in an earlier paper (Treloar 1947 b). 


§2. EXPERIMENTAL METHOD 


In the present paper a study of the stresses and birefringence in the most'general 
type of homogeneous strain is reported. For this purpose a sheet of vulcanized 
rubber, compounded according to the formula given in the paper referred to 
(Treloar 1947 b), was stretched in two directions at right angles in such a way 
that the strain was homogeneous, whilst the two principal extension ratios A, and 
A, in the plane of the sheet could be varied independently by the application of two: 
sets of forces. The arrangement used permitted the simultaneous measurement 
of the principal stresses and the double refraction corresponding to any state of 
strain. 

It is not easy to devise a suitable method of applying two unequal strains to: 
rubber in two perpendicular directions. The method here adopted, though not 
ideal, particularly in that the maximum attainable strain was limited, seemed to 
provide the simplest means of achieving the desired result. It makes use of 
the principle of the guard plate to eliminate the effect of the non-uniformity of the 
strain in the vicinity of the edges of the sheet. ‘The shape of the test piece is shown 
infigure1. It was in the form of a square, with five projecting lugs on each side, 
to which strings were tied for the application of the loads. The surface of the sheet 
was marked out with a series of lines forming a square lattice, so that the state of 
strain at all points could be observed. With the sheet placed horizontally, the 
three middle lugs on one side were loaded by means of three equal weights attached 
to strings passing over pulleys, whilst the strings attached to the two outermost lugs. 
were secured to a rectangular frame, mounted horizontally. The five strings on | 
the opposite side were connected to the opposite side of the frame. Similarly, the 
three middle lugs on an adjacent side carried three equal weights, whilst the two, 
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outer lugs, as well as the five on the side opposite to them, were connected to the 
frame. The 14 strings connected to the frame were provided with tension 
adjusters, and these strings, as well as the pulleys, could be moved sideways 
independently, so that each set of strings could be kept parallel and the two sets 
maintained perpendicular to each other. ‘The appearance of the stretched sheet 
when the tension on the outer lugs was suitably adjusted is shown in figure 1; 
in all cases the strain within the rectangle ABCD was sensibly uniform. Under 
these conditions it was assumed that the stresses acting on the sides of the rectangle 
ABCD were defined by the loads applied to the three central lugs, and that the 
non-uniform edge region, though requiring the use of a slightly different tension 
on the outer lugs *, would be without effect on the stress in this inner, uniformly- 
strained region. 

For the optical measurement a beam of light from a polarizer set at 45° to the 
directions of principal strain was passed normally through the sheet, the birefrin- 
gence being measured by means of a Babinet compensator. 

The measurements were made at room temperature, which varied between 
17 and 22°c. 

93. EXPERIMENTAL RESUS 

The data obtained are givenintables1and2. ‘Table 1 refers to the dry rubber, 
and table 2 to the same rubber swollen in medicinal paraffin to v,=0-525. 
Readings were taken in the order given in the tables. Generally about 10 minutes 
elapsed between the application of a load and the taking of readings. ‘The 
forces f, and f, are the loads in grams applied to each of the three middle lugs; 
the areas of cross-section of the rubber in the unstrained state (Ay), over which these 
forces act are given in the table. The figures for optical retardation are given in 
terms of the number of sodium wavelengths phase difference between the two 
rays, this being the quantity actually measured. ‘To convert this to difference of 
refractive indices (7, — 7) it is only necessary to divide by the thickness of the sheet 
and multiply by the wavelength. ‘The extension ratio A;, and the corresponding 
thickness in the strained state, were determined from the measured A, and dg, 
assuming no change of volume on straining. 

In addition to the data represented in tables 1 and 2, figures were obtained for 
the case of simple elongation, i.e. with f,=0, using a parallel strip of rubber cut 
from the same sheet. ‘These additional data are included in the graphical present- 
ation of the results. 

The optical data are presented in figures 2 and 3, omitting points (near the 
origin) corresponding to the case when f,=/;. In the first case (figure 2) the 
birefringence is plotted against the difference of the squares of the corresponding 
extension ratios, A,2—A,2._ While the expected linear relationship (equation (1)) 
is approximately borne out, there are for the dry rubber, slight departures which 
are probably greater than the errors of measurement. For the swollen rubber, on 
the other hand, the agreement with the theoretical form is very close. In the 
second case (figure 3) the birefringence is plotted against ¢, — fg, the difference of the 
principal stresses. Again, the proportionality is more exact for the swollen than 


for the dry rubber. 


* It was found experimentally that the difference of tension on the outer strings was in fact 
very slight. 
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Table1. Dryrubber. Thickness=0-082cm. Ay =0-0648 cm? 


Retardation 
fi (gm.) fe (gm.) Ay rz (sodium 
wavelengths) 
100 100 1:07, 1-08; —0:02 
100 200 0:97, 1-297 —0°55 
200 200 1:16, 1-205 —(0-07 
300 200 1:46, 1:09, +-0°54 
300 100 1:58, 0-90; +1-06 
300 300 le Sing LA —0-05 
300 400 1:18, 1-88 —0-80 
200 400 0:93, 2-04 —1:50 
100 400 0:78 2-14 —2-00 
400 400 1-60 (LOW —0-13 
500 400 2:24 (lebiils, +0:64 
500 300 2:42 1-09, +1-39 
500 200 2°51 0-87; +2-05 
500 100 2:64 0-72, +2-80 
500 500 2:07 2-10 —0-05 
600 500 2:68 1:87; +0-54 
400 600 1e2 55 2-98 —1:55 
300 600 0-98, 3-06 —2:31 
200 600 0:79, 3-14 —3-05 
100 600 0:67, 3-21 — 3-76 
600 600 2:34 2:63 —0-18 


Table 2. Rubber Salba to v,=0-525 in paraffin. Thickness 
(swollen) =0-:099cm. A )=0-0985 cm? 


Retardation 
Fi (gm.) fz (gm.) Ay Az (sodium 
wavelengths) 
100 100 1-08, 1-08, —0-:02 
100 200 0:99, 1-309 —0-40 
200 200 1-21 ie), +0-:01 
300 200 1S 1:06, +0-51 
300 100 1-66 0-88, +0-93 
300 - 300 1-41 1-36, +0-03 
100 400 0-80 2-04 —1-46 
200 400 0:95, 955 —1-02 
300 400 1-24 1-83 —0-51 
400 400 1:65 L655 0-00 
400 500 1:50 2:12 —0-46 
500 500 1-97 1:94 +0-02 
500 300 Deel 1-10 +1-06 
500 200 2:43 0-85, +1-62 
500 100 2:48 On/Sie +2:07 


Figure 4 represents the difference of the two principal stresses, plotted against 
A,?—A,?, for both samples. As with the birefringence, the swollen rubber appears 
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Figure 2. Birefringence plotted against difference of Figure 3. Birefringence plotted against 
squares of principal extension ratios for dry and difference of principal stresses for 
swollen rubber. The crosses refer to simple dry and swollen rubber. The 
elongation. crosses refer to simple elongation. 


to agree with the theoretical form (2) within the experimental error, but the 
agreement in the case of the dry rubber is not exact. 


§4. FURTHER EXAMINATION OF STRESSES 

Whilst the evidence presented in figures 2, 3 and 4 suggests that the theoretical 
expectations regarding the dependence of the birefringence and of the principal 
stress difference on the two principal extensions in the plane of the sheet are 
approximately fulfilled, particularly in the case of the swollen rubber, an examin- 
ation of each of the principal stresses separately reveals consistent deviations from 
the theoretical form, as will now be shown. 

Since the stress f; normal to the sheet is zero, we should have, according to the 
theory previously outlined (‘Treloar 1947 a), 

= GG —=Ns)@ eee cee (2a) 

with a corresponding expression for ¢,. According to this equation, a plot of ¢, 
or ty), against A?—A3 or A3—AZ respectively, should yield a straight line. ‘The 
actual results, shown in figures 5 and 6, present a more complex appearance for 
both the dry and swollen rubber. 


30 T 


s (pli (kg/cm?) 8 


6 
N- 2% 0 2 4 non 6 8 10 
Dry rubber. Stress t, plotted against 
\y2—A;2. The black circles refer to simple 
elongation. f, and f, are expressed in grams. 


Figure 4. Difference of piincipal stresses plotted Figure 5. 
against difference of squares of principal ex- 
tension ratios. The crosses refer to simple 
elongation. 
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The method of plotting employed in these diagrams requires a little explan- 
ation. Firstly, it must be remembered that the stress (say ¢,) is determined by the 
force (f,) divided by the area of cross-section, which is itself a function of the strain, 
i.e. t,=A,f,/Ap. Hence, while one of the forces is held constant, the stress 
corresponding to this force varies with changes in the second force, since such 
changes alter all the extension ratios. In order to simplify the presentation, the 
stress plotted is referred to as ¢,, and the corresponding force as f;. Thus two 
points appear on the diagram for each state of strain, one for each of the two 
principal stresses. 

It is seen that the points fall on a set of lines, on any one of which f, is constant, 
while f, varies. This clearly means that ¢, is dependent on a more complicated 
function of the extension ratios A,, A, and; than the one represented by equation (2). 

This important departure from the theoretical form exhibited when the prin- 
cipal stresses are plotted separately is rather surprising in view of the fair agree- 
ment obtained when only their difference is considered. Before discussing its 
significance, it seems desirable to examine very carefully the nature of the evidence, 
in order to make sure that it is not due to some spurious effect introduced by the 
experimental technique. Three possible sources of trouble will be discussed, 
namely (a) relaxation, (6) anisotropy in the original sheet, and (c) non-uniform 
strain round the edges. 


(a) Relaxation 


It is true that the attainment of equilibrium values of the stresses in strained 
rubber is a matter of considerable difficulty, requiring the breaking-down of 
temporary intermolecular ‘‘cross-links”” by the use of higher temperatures or 
other means. No attempt was made to obtain genuine equilibrium values in the 
present experiments, and it is therefore necessary to consider whether the dis- 
continuous array of figures 5 and 6 might arise from this cause. 

There are two reasons for rejecting this hypothesis. Firstly, although it was 
observed that the values of the strain parameters obtained under the application of 
a given pair of forces depended to some extent on the order of application of the 
forces and on the previous history of the sample, nevertheless, the stresses, when 
plotted as in figure 5, invariably fall precisely on the appropriate discontinuous 
arrays. Secondly, it is clear that swelling should favour the approach to equili- 
brium, but the discontinuities are no less marked for a rubber swollen with 
nearly 100% of solvent than for the dry rubber. 


(b) Amisotropy in the original sheet 


It is conceivable that the departures from the theoretical form might arise from 
a difference of properties in different directions in the original sheet of rubber, 
introduced by the rolling process prior to vulcanization. Such an effect would 
arise from non-equilibrium or secondary linkages rather than from primary 
linkages, since the form of the stress-strain relations derived from the molecular 
theory is in no way dependent on the detailed structure of the molecular network 
(James and Guth 1943). If such anisotropy were present, it would be expected 
to reveal itself in a non-uniform swelling in different directions; but no such 
non-uniformity in swelling was observed. Furthermore, a latex rubber sheet, 
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prepared by a quite different process not involving rolling or pressing, gave 
a precisely similar type of diagram to that of figure 5. It is therefore concluded 
that an initial anisotropy cannot be responsible for the essential features of this 
diagram. 

(c) Non-unifornuty of strain 

It has already been suggested that the guard-plate principle used in these 
experiments should eliminate any effect of non-uniformity of strain in the neigh- 
bourhood of the edges of the sheet. 

In order to check this point, experiments were made by two independent 
methods. . ‘The first was by the use of parallel strips in simple elongation, using 
values of force per unit area of the unstrained section equal to the values of f, and f, 
in the tables. ‘The points derived from these experiments, separately marked in 
figures 5 and 6, will be seen to fall exactly on the previously obtained arrays. 
The second check was obtained from measurements on a circular sheet clamped 
round its circumference and inflated. Details of the method have been described 
elsewhere (Treloar 1944, 1944b). This corresponds to the special case where 
A,=A,. The experiment was made on the actual sheet of swollen rubber used to 
obtain the data of figure 6. The results cannot conveniently be shown on this 
diagram, because the force could not simply be chosen to correspond to the values 
of f, and f, previously used. A plot of ¢, against A{—A?, however, agreed to 
within 3% with the data in figure 6 for the case f, =/fy. 

These independent experiments relating to the two extreme cases of simple 
elongation and 2-dimensional extension with A, =A, prove conclusively that the 
unexpected form of the results is not due to some defect in the technique of 
measurement. All the evidence suggests that the observed phenomenon repre- 
sents a fundamental property of the rubber. 


So LHP ORE PIiCAL INDE REPRE LAT TON OF DATA 
Mooney (1940) has discussed the form of the general stress-strain relations for 
a rubber on the basis of an assumed stress-strain relation in simple shear. For the 
particular case of a linear stress-strain relation in shear, he obtains, for the work of 
deformation 


ety a al 
W,=C(+B+9—3) +0255 + 55 + 8-3); een (3) 


where C, and C, are physical constants of the material. ‘The more general case 
of a non-linear shear relation may be resolved by adding further terms to equation 
(3). ‘Thus the next term in the series would be of the type (Mooney, loc. cit. 
equation (40)) 


if 1 1 
We A+ +8 —3)4B 5+ at 3); oat (4) 


and the work of deformation would then be 
W=W,+W,. 


From a rather more general standpoint, Rivlin (unpublished work) has argued 
that the work of deformation, or stored-energy function, must be a function of 
certain ‘‘strain invariants”, which are themselves symmetrical functions of the 
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even powers of the extension ratios Aj, A, and As. The simplest functions which 
may be chosen for the stored energy are 


W,=GOA2+A2 403-3), ates (5) 
and Wr=K(i +53 +g -3): Res 2 (6) 
Next in order of complexity are 
Wy=L(AgtAgt+rg—3), ne ees (7) 
gl gil: 3! 
and Wy=M (p+ + 3): socen (8) 


Rivlin’s approach is perhaps simpler, and avoids the assumption that the form 
of the stress-strain relation for shear in an isotropic plane is independent of a 
stretch in the direction normal to that plane. However, the resultant stress- 
strain relations are the same whichever theory is adopted. 

The stored-energy function derived from the molecular theory (Treloar 
1947 a), from which equations (2) and (2 a) are obtained, is 


WaGOLLMN EN 3), 6 en eee (9) 


and is therefore equivalent to Mooney’s equation (3) with the constant C, put equal 
to zero. 

Let us now examine the stress-strain relations resulting from the first equation 
of Mooney (equation (3)). Corresponding to this stored-energy function, the 
principal stresses are given by equation (25) in his paper, which may be written in 
the form * 


i= =(GL RX) ee (10) 


which, as will be seen, is a rather more complicated expression than equation (2) 
derived from the molecular theory, using the stored-energy function (9). In 
this expression G and K are physical constants expressible in terms of C, and 
C, in equation (3). In the case when the forces per unit area of the unstrained 
rubber, Ff’, and F, are given, and F,=0, we have the following equations relating 
the strains with the forces (or stresses) 


c= AVF, ae (G is KAS) F A3); 
p= Ay F=(EL RCC >) eee (10a) 


These simultaneous equations enable the extension ratios to be determined when 
the applied forces F, and F, are given. Their solution has been obtained graphi- 
cally, using the particular values G=1-0, K=0-1, for six values of /, combined 
in all possible ways with the same values of F,. This solution is represented 
in figure 7 in terms of the stress ¢, plotted against Aj—A3. Detailed comparison 


with figure 6 shows a very close quantitative correspondence between the 
theoretical and experimental points. 


Another method of comparing the experimental data with the theoretical 


equation (10) is by plotting the experimental stress ¢, against (1 + KA3/G)(Aj —A3), 


* 'The constant G is not the same as Mooney’s G. 


Stresses and birefringence in rubber 143 


choosing a suitable value of K/G. This type of plot is shown in figure 8, in which 
curve (b) relates to the swollen rubber. With K/G=0-1, the points fall very 
nearly on a straight line. 

Thus Mooney’s equation (3) gives a completely satisfactory representation 
of the behaviour of this particular rubber in the swollen state. 

It is interesting to examine some further properties of Mooney’s equation. 
Firstly, for a pure shear, we have one dimension (say A,) unchanged, and equation 
(10) reduces to 

t, =(G+ K)(Az—A}). 
In this case, therefore, the behaviour is indistinguishable from that given by 
the simple equation (2a). Secondly, for the difference of the two principal 
stresses, in the general homogeneous strain, we have 


PGR ee eee (105) 
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Figure 6. Swollen rubber. Stress t, plotted Figure 7. Stress calculated from Mooney’s: 
against A\?—A,;?._ The black circles refer equation, using G=1:0, K=0-1. 


to simple elongation. 


If the forces are applied in the directions A, and A, the extension ratio corre- 
sponding to the thickness direction is always less than unity, and A} is generally 
quite small. If, therefore, K is also small compared with G, the term KA3 
may be neglected, except at very small strains, i.e. equation (105) reduces to the 
simpler form (2). ‘This accounts for the observed linear relation between ¢, —t, 
and A? — 23. 

Turning now to the dry rubber, we find rather a more complicated picture. 
By a suitable choice of K/G it is possible to represent the stress ¢, as a continuous. 
function of (1+ KA3/G)(Aj—A3), but this function is not linear (figure 8 (a)). 
From this observation it is concluded that Mooney’s equation represents a part, 
but not all of the deviation from the molecular theory (equation 2a). ‘To account 
for the whole discrepancy it would probably be necessary to add further terms 
to the stored-energy function, such as those represented by equations (7) or (8). 
Since Mooney’s equation (3) is based on a linear stress-strain relation in simple 
shear, it follows that the dry rubber is non-linear in shear. This fact may also 
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be deduced by direct interpolation from the experimental data. In other rubbers 
examined by the author, a non-linear shear relation has been obtained (Treloar 
1944), but Mooney (1940) has reported a linear relation up to 200°% shear. 


§6. CONCLUSION 

The conclusion to be drawn from the preceding analysis of the stress-strain 
relations may be summarized in the following way. As a first approximation 
the equations derived from the statistical treatment of a molecular model provide 
a basis for the interpretation of the elastic properties of rubber. A closer approxi- 
mation is obtained by including an additional term in the stored-energy function, 
i.e. by the use of Mooney’s equation (3). This appears to give an accurate 
representation of the properties of a swollen rubber, but is still inadequate when 
applied to a dry rubber, for which a third approximation, including at least 
one more general term, is likely to be required. 


30 
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Figure 8. Data from figures 5 and 6 plotted with (1+ KA,?/G)(A,?— A3) as abscissae. 
(a) Dry rubber, K/G=0-05. (6) Swollen rubber, K/G=0-1. 

It is necessary to emphasize that the use of the second or higher approximations 
does not in itself throw any light on the physical mechanism responsible for the 
observed behaviour. It is, in fact, the 3-dimensional analogue of simple curve- 
fitting. It must also be borne in mind that any physical model that might be 
postulated must necessarily lead to results which may be represented by some 
combination of the stored-energy functions discussed in § 5. 
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ABSTRACT. Measurements by a number of different methods of the dielectric constant 
and loss angle of water and heavy water at three widely separated wavelengths in the region 
of anomalous dispersion, A=10cm., 3 cm., 1:25 cm., are described. Results show that 
the Debye equations are exactly obeyed, there being a single time of relaxation varying with 
the temperature ; the ratios of the relaxation times of water and heavy water are found to be 
in quantitative agreement with the ratios of the viscosity over a temperature range, suggesting 
that the re-orientation mechanism in these liquids is the same as that of viscosity. The 
value of the activation energy of the re-orientation process is discussed. The value of 
55 found for the optical dielectric constant gives a reasonable result for the dielectric 
constant of water on the Onsager theory. This theory is discussed and shown to give a 
relaxation equation of the Debye type for alternating fields. 


Soe NER ODI C moO) 

INCE 1900 the anomalous dielectric dispersion of water at very short wave- 

lengths has been widely investigated, and has been summarized by Dorsey 

(1940). ‘The main qualitative features of the phenomenon are no longer in 
doubt, but quantitatively the results of different workers show wide divergencies. 
The significance of the. problem lies in its theoretical bearing on the structure of 
water, which is not only of fundamental importance but also of technical and 
biological interest. In addition, a knowledge of the dielectric properties of water 
is needed in calculations on the transmission of microwaves. 

It therefore seemed important to apply the recently developed high-frequency 
techniques to obtain an accurate measurement of the refractive index n and 
absorption coefficient « of water. It was also felt that parallel measurements on 
heavy water would provide comparative data of theoretical value. In this paper 
measurements are described at three widely separated wavelengths, namely 
10-0cm., 3cm., and 1:25 cm., in the temperature range 0°-75°c. 

The theory of an ideal polar dielectric in an alternating field was first given by 
Debye in terms of a single relaxationtime. It can be shown that the relation below 
holds generally for linear dielectrics, as is discussed later in this paper: 
ea, (1) 
1+jur’ 
in which 7 is a time of relaxation, and « is the complex dielectric constant at the 
angular frequency w; €, is that part of the dielectric constant which is due to the 
atomic and electronic polarization, and is assumed to be real and independent of 


w. The static dielectric constant is «,. This expression may be expressed in 
terms of the optical constants as follows: 


(3S Ay 


Es ie €9 4h 6e= Es - i) at € 
depo a LAVA 
eee (<, <a €y)@T ms (€, = €9)A,/A 
lta 1+ AP’ 
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2nk =e'' 


146 C. H. Collie, #. B. Hasted and D. M. Ritson 


where is the real part of the refractive index, « the absorption coefficient, «’ and 
<’’ the real and imaginary parts of the complex «, A the wavelength corresponding 
to a frequency w, and A, the wavelength which corresponds to frequency 1/7 
(the so-called ‘‘ Sprungwellenlange ’’). 

These equations represent the fall of the dielectric constant from its static 
value «, to its optical value «o, the fall being accompanied by a single broad 
absorption band in the neighbourhood of the characteristic wavelength A,. The 
present measurements were made with the object of finding out whether the 
simple expression with a single relaxation time does adequately represent the 
change in optical constants with frequency at a fixed temperature. As will be 
seen, this is in fact the case, so that one can also investigate : 


(1) the variation of A, with temperature. ‘T’his can then be compared with 
other physical properties, such as the viscosity, whose temperature 
dependence is known. 

(2) the ratio of A, for water and heavy water. 


§2. METHODS OF MEASUREMENT 

The optical constants vary so rapidly with wavelength and temperature that 

no universal method is available. 

The following methods, some of which have been described in previous 

publications, have been used :— 

(1) The change of the resonant frequency and damping of an Eg) resonator, 
when an axial capillary tube filled with water is inserted, is measured 
(Willis Jackson 1946, Collie, Ritson and Hasted 1946). 

(2) A similar measurement using an Ho, resonator (Penrose 1946, Collie 
Hasted and Ritson 1948). 

(3) Measurements of the absorption coefficient are made by direct trans- 
mission in a waveguide filled with the liquid under examination (figure 1, B). 
If two measurements are made of the absorption coefficient, one, «,, in 
a rectangular guide of width a,, and a second, «sy, in a guide of width a, 
near cut-off, then it is easily shown (loc. cit.) that 


21a?) (R2«? — R2«3 
Ries (C/o M¢ ae gD) (x2 +43), 
Si Sy 


é’ = 2k, fe’ — R32 c2/w? +K7]}, 


in which k, = 7/a, and ky =77/a, for the Hy) wave in a rectangular waveguide. 

This method affords a satisfactory way of measuring bothnandx. As 
can be seen from the form of the equation the method is a practicable 
adaptation of the ideal method of determining « by determining the atten- 
uation in free space, and m by measuring the critical wavelengths of 
propagation down a guide. ‘The quantities which have to be measured 
accurately are the width of the guide and the wavelength. These are the 
two quantities which can be measured with the greatest precision, and the 
latitude in the choice of a, and a, is sufficient to allow the two measurements 
of the attenuation in the guide to be approximately independent measure- 
ments of n and « the optical constants in free space. This method is not 
limited to rectangular waveguides, and some of the later measurements 
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have with advantage been made using an H,, wave in cylindrical waveguides 
of radii 7, and 7,, in which case k, = 1-841/r,, and k, = (S41 i7,; 


(4) Direct measurements of the absorption in a sufficiently wide waveguide 
yield « directly (figure 1 A). 


Power out 
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mevement 


Figure 1. 
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The methods which have been used are shown in table 1. 


Methods of determination of m and x 


Principle of absorption measurement at 1 cm. and 10 cm. 


D,0 


Absorption in 
two rectangular 
guides. 


Absorption in 
two circular 
guides. 


Absorption in 
two rectangular 
guides. 


Absorption in 
two circular 
guides. 


Table 1. 
H,O 
| Wavelength n kK 
(cm.) 
10 Eoi9 resonator. 
Absorption in | Absorption in 
two rectangular | two rectangular 
guides. guides. 
Absorption in | Absorption in 
two circular | two circular 
guides. guides. 

Bey Absorption in | Absorption in 
two circular wide rectangu- 
guides, lar guide. 

Absorption in 
two circular 
guides. 

125 Ho, resonator. Absorption in 


§ 3. 


wide rectangu- 
lar guide. 


Hp, resonator. 


Absorption in 
wide rectangu- 
lar guide. 


MEASUREMENT OF 2 AND xk AT 10 CM. 


These measurements are relatively easy to make using radar receiving equip- 
ment whose main features were standardized during the war. 
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Power from a well screened reflex klystron (type CV 35) is led through a 
coaxial lead to a coaxial socket of the type first used by the General Electric 
Company in their research laboratories at Wembley. Commercial pyrotenax 
cable (type 147/1) is very suitable for the coaxial transmission line since the con- 
tinuous outer coating facilitates the very high degree of screening needed if the 
fullest use is to be made of the high sensitivity obtainable (10~" watt) with a 
superheterodyne receiver. 

From the coaxial plug the power reaches the receiver, either through the 
absorption cell or through an H,, piston attenuator. 

The construction of these absorption cells is shown diagrammatically in 
figure 1 B and the absorption cell using circular guides is shown in detail in figure 2. 


Figure 2. Figure 3. Figure 4. 
10-cm. absorption cell. 1-cm. absorption cell. 3-cm. absorption cell. 


Power from the coaxial line E is fed into the rectangular wave guide K by means of a 
probe; a match is obtained with two stubs 3-75cm. apart at H. A small fixed 
loop in this section enables the power level to be monitored on a crystal detector. 
From the top of the guide the power passes through a }mm. mica window into a 
circular guide filled with water. The window is held in position by a plate 
secured by three countersunk screws. The joint is rendered watertight with 
i eae silicone stop-cock grease which provides an effective seal even at 

ac 

The interchangeable circular guide is made from copper tube straightened ona 
silver steel mandrel. ‘The two diameters used were 1-261cm. and 0:-576cm. As 
far as could be judged by measurement of the open ends and the fit of a steel plug 
the guides were uniform to 1 part in 200, and this is confirmed by the absence of 
systematic variations of the measured absorption coefficient. 

After transmission through the water the power is again picked up by a loop at 
the end of a movable pyrotenax coaxial cable C. The loop can be withdrawn and 
its position read by means of a 12cm. screw movement A. 
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The loop is held central in the guide by means of a loosely fitting plunger (see 
inset). The loop extends about 14mm. into the liquid and is nearly the same 
diameter asthe tube C. Spacers can be inserted into the movement to allow differ- 
ent portions of the wave guide to be explored by the loop. 

A subsidiary electrode L protruding from the cable C enables the conductivity 
of the water to be measured during the course of measurement. The absolute 
conductivity is obtained by calibrating the apparatus with N/50 KCl. The true 
value of ”’ is obtained from its apparent value <,’ by the relationship €;’ =" +470/w 
in which o is the D.c. conductivity in £.s.u. If conductivity water is used and 
the joints are made so that as little solder as possible is in contact with the water the 
conductivity correction can be kept below 2% of <’’ at 80°c. The magnesium 
oxide filling of the pyrotenax must be prevented from coming into contact with the 
water by a well baked layer of shellac varnish. 

For measurements at 3-2 and 1-25 cm. this correction is sufficiently small to 
justify the omission of a continuous measurement of the conductivity. ‘This was, 
however, always measured at the end of a run to verify that excessive corrosion had 
not occurred. 

At all frequencies freshly boiled water must be used to prevent the formation of 
bubbles. For the same reason the temperature should be progressively lowered 
during a series of measurements. The heavy water used throughout was of 
99-5% purity. 

All the relevant parts of the apparatus were water jacketed, and water from a 
thermostat was circulated through the jackets through lagged pipes. 

The temperature was measured in the water stream by a thermometer which 
had been compared with a set of thermometers by Anschiitz, which had themselves 
been compared with thermometers carrying an N.P.L. certificate. 

The measurement of the absorption coefficients is made in two parts. First 
the receiver responses k,, R,....R, are measured for a series of positions of the 
loop in the absorption cell. ‘The piston attenuator is then substituted for the cell 
and a measurement made of the position of the attenuator loop at which the receiver 
responses are the same as one of the previously noted values R,, R....R,. To 
verify that the power output of the transmitter and the sensitivity of the receiver 
have not changed during 


the measurement, the re- | | Wide circular 

sponse R, is checked after © hg 
every reading of R,. In ae = 
this way it is possible to a : : 
plot a curve whose ordinate = a : 
is the displacement of the co ie 
loop in the absorption cell 2 : 
and whose abscissa is the €&, Mees 
displacement of the loopin 2 rectangular a g 
the H,, attenuator which §& , 1-25 cm.59°86. i 
produces the same change = 

in receiver response. This ~ ef | 1 | - 
curve should be a straight siteniet or ant a ae 

line whether the receiver 

response is linear or not. Figure 5, ‘T'ypical absorption plots. 
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This offers a valuable check and enables errors in the experimental arrange- 
ments to be detected and eliminated. The two most obvious sources of 
error are pick-up and the production of unwanted modes in the absorption cell. 
Since the former is most troublesome at low signal strengths and the latter when the 
path in the water is small and the signal consequently strong, it is very improbable 
that these sources of error will compensate each other. ‘Typical straight lines are 
shown in figure 5. The slope of these lines was evaluated by the method of least 
squares and the absorption coefficient calculated from the theoretical attenuation 
of the attenuator. 


S45 PHB RE Gy Vara 
An Admiralty receiver type P 50 was used in these measurements: The band - 
width of 3 Mc/s. was sufficient to embrace the small frequency wanderings of the 
transmitter and local oscillator klystrons. Although more complicated than a 
simple crystal rectifier and sensitive galvanometer it has so many advantages that 
its use is almost essentual if an overall accuracy cf 1% is aimed at. 
The main advantages are :— 


(2) A constancy of calibration lasting for several hours, and in favourable 
circumstances over a period of days. 


(b) A factor of at least 10* in the usable sensitivity. 
The increased sensitivity is of great importance in the type of work being described 
because: 
(i) It enables all coaxial leads to be of sufficient length (about 2 m. of pyrotenax) 
to damp out the standing waves due to inevitable mismatches at plugs and 
junctions. e 


(11) It enables power levels to be measured with a piston attenuator; otherwise 
the use of these admirable instruments is precluded by their large insertion 
loss. 

The use of an attenuator allows the receiver to be used as a monitor so that it is 
sufficient to read the second detector current on a 0-100 microammeter. A 
suitable resistance net work allows a 40 db. range to be covered without changing 
the gain of the 1.F. stages. 


§5. MEASUREMENTS AT 1-274 CM. 


Accurate measurements at this wavelength offer some difficulty since the 
degree of frequency stability required is relatively much greater if an 1.F. amplifier 
of 3 Mc/s. bandwidth is to be used, while the valves are inherently less frequency 
stable than klystrons operating on a longer wavelength, owing to the very small 
clearances and high power dissipation involved. Satisfactory measurements were 
made possible by modulating the frequency of the output klystron by means of a 
10-40 volt saw-tooth voltage of recurrence 500 c.p.s. applied to the reflector. The 
output of the receiver is then to a first approximation independent of the absolute 
frequency, since the 1.F. frequency now sweeps throught the whole effective 
spectrum of the LF. receiver. Satisfactory results were obtained provided 
outside conditions were kept as constant as possible. This was done by enclosing 
both the output and local oscillator klystrons in copper screening boxes cooled by a 
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constant stream of compressed air; in addition the voltage control was improved 
by replacing the neon stabilizers in the power pack by dry batteries. Some 
attempts were made to use the “two waveguide” method which had proved so 
convenient at 10cm. It was not very satisfactory owing to the difficulty of match- 
ing into and out of a waveguide only 0-9 mm. in diameter, and the method was 
abandoned. The refractive index n was measured by the resonator method 
previously described (Collie, Hasted and Ritson 1948), and « was determined 
by direct measurement of the attenuation in the absorption cell, shown in 
figure 3, and diagrammatically in figure 1B. The water is contained in 
the upper part of the vertical rectangular guide G of internal cross section 
10-75 mm. x 4-5 mm., supported ona thin (0-15 mm.) mica window firmly clamped 
and greased in position. After traversing the water the wave is picked up by a 
movable guide (internal cross section 8-75 mm. x 3-5 mm.) fitted with a distrene or 
rubber window about 2mm. thick. This type of pick-up is very satisfactory and 
is much easier to match into the water than the loop and coaxial line used at longer 
wavelengths. The wave is led to the mixer by the right-angle bend at D and a 
length of flexible waveguide. A pair of stubs provides sufficient matching. 
In view of the rather open design of the cell it was thought worth while to measure 
the difference in temperature between the water in the cell and the water circulating 
in the water jackets. Measurement with a thermocouple showed that after ten 
minutes at any fixed temperature the difference in temperature between the water 
in the cell and the thermometer chamber was less than 0:1°c. The attenuation 
produced by a known thickness of water was measured as before by direct com- 
parison with a piston attenuator. 

The greatest source of error in absorption measurements of this sort lies in the 
possibility that new modes are created when the wave passes through the window 
into the water-filled guide. Each mode will be differently attenuated so that little 
precise information can be obtained from the measured attenuation curves, which 
are not exponential. In these experiments the production of unwanted modes was 
minimized by using a thin window whose flatness could be investigated by reflected 
light after it had been mounted inthe cell. Distortion due to differential expansion 
was avoided by greasing, but not sticking, the window to the copper flange. A 
watertight joint up to 80° was ensured by using non-melting silicone grease. The 
freedom from unwanted modes, and overall accuracy of measurement, is shown by 
the achievement of straight-line plots for cell movement against piston attenuator 
movement. ‘l’ypical plots at all wavelengths are shown in figure 5. 

At this wavelength good straight lines through more than 3 mm. of water were 
obtained with this design of cell. There is, however, disagreement between these 
results and those we reported in an earlier paper, which is due to two facts: — 


(i) The design of the earlier cell did not permit the propagation of pure modes, 
so that curved logarithmic absorption plots were obtained; this curvature 
was erroneously described as geometrical in origin. 


(ii) In order to eliminate this effect measurements were taken with a large dead 
space of liquid, but through small thicknesses, with very high receiver 
sensitivity. Repetition showed that these measurements were vitiated by 


pick-up. 


Lite} 
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In the revised method both these sources of error have been eliminated beyond 
reasonable doubt, and the new results are taken as accurate to +1%. 


§6. MEASUREMENTS AT BpoWaey (IMC 


Measurement of « at this wavelength was carried out with a rectangular wave 
guide cell similar to that described for use at 1:274cm. It is shown in figure 4. 

A type CV 87 reflex klystron was used as the source of power; this valve has a 
coaxial line output and is connected to a waveguide system by a probe. After 
passing through a unidirectional feed and a flap attenuator the power is matched 
by stubs through a thin mica window G into the water-filled cell. ‘The small 
amount of power (about 1°%,) led off by the one way feed is coupled into a high-Q 
resonator, the output of which is fed to a crystal rectifier and a high sensitivity 
galvanometer. In this way both the power level and the frequency stability of 
the wave entering the cell are continuously monitored throughout a series of 
readings. 

The method of measurement was similar to that used at 1:274cm. The cell 
was so designed that a number of fine cylindrical waveguides could be substituted 
for the rectangular guide, and a concentric line and pick-up loop for the output 
guide. These parts are shown in figure 4, and when they are in use the apparatus 
is similar to that used at 10cm. for the measurement of both m and x. The 
diameters of tube used were 2:51 mm.and3-77mm. As can be seen from figure 7, 
excellent agreement was obtained between values of x obtained by the two methods, 
and the fine cylindrical wave guide has been used throughout for « measurements 
on heavy water, since less liquid is required. he cylindrical tubes were used 
throughout for measurement of 7. 


§7. MEASUREMENT OF WAVELENGTH AND ATTENUATION 


The wavelengths used were 1:27,cm., 3-21,;cm., and 10-00. cm. respectively. 
The two former were measured with a simple absorption wave-meter using 
an H,, mode. The absolute accuracy depends upon the diameter of the reamed 
hole. They are uncorrected for the dielectric constant of air, and probably 
accurate to better than }°/,. The 10cm. wavelength was measured with a high 
precision resonant cavity wavemeter (Admiralty type G93), which could he 
calibrated against a crystal-controlled oscillator. 

Three inductive piston attenuators using an H,, mode were used in these 
measurements. Their attenuation constants are given below :— 


Attenuator, 
Way cooet diameter A 
(cm.) (mm.) 
1-274 4-00. rae 
3-213 9-525 ae 


‘The attenuation constant « is calculated from the relationship 


__ displacement of the attenuator 
displacement through the water 
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The attenuation constant A is calculated from the relationship 
— 417? A? /d? + (1-841)?/7? = w?/c?. 


An account of the theory and construction of these instruments has been given 
inter alia by Clayton, Houldin, Lamont and Willshaw (1946). 


Wo, WS SWUML ANS 

Values of m and « have been obtained at a sufficiently large number of temper- 
atures to enable smooth curves to be drawn through the points and values at ten 
degree intervals from 0° to 75°c. for water and from 5° to 60°c. for heavy water to 
be given at each frequency. ‘These values, which are estimated as accurate to 
+1%, (10:1 probability) and + 0-4°c. in temperature, are given in tables 2 and 3, 
together with the values calculated from a single value of A, which is also tabulated at 
each temperature. Some measured values at 3 cm. and the smooth curves through 
them are shown in figure 7. It will be seen that except for a very few values there is 


Table 2. Dielectric properties of water 


éy=9°9 
IKKECe €s Ns X=1-27 cm. A=3-21 em. A=10-:0 cm. 
Px n ix K n K n K 
0 4-94 2:86 28 2°86 9-03 WeS7/ Gay. 


88-2 8°34: Api 2°81 esl 2°83 9-04 | 1-37 calc. 


DED) oil deus 2°43 8:89 | 0-982 exp. 


a 84-2 2:39 5-60 | 2-92 | 7-88 | 2:39 9-00 | 0-989 calc. 

6-25 | 2-86 | 8-08 | 1-97 | 8-83 0-739 exp. 

| i 80°36 | 1-80 6-20 | 2-83 | 8-46 | 1-96 | 8-87 | 0-737 calc. 

| 6-75 | 2-67 | 841 | 1-57 8:78 | 0-558 exp. 

o 76-7 1:39 6-71 | 2-64 | 8-25 | 1-58 8-70 | 0-559 calc. 

7-0 2:40 | 8-20 | 1-29 8-53 | 0-442 exp. 

| * 73-1 1-12 702 | 239 | g-21 | 1-28 8°51 | 0-440 cale. 

| | 7.26 | 211 | gas | 1-05 | 8-28 | 0-350 exp. 

= 69-8 (O:O02 1007-23: | 2:10) 8-12 ||) 1-04 8-33 | 0-348 cale. 

iS. 7:30 | 1:87 | 7-99 | 0-865 | 8-09 | 0-281 exp. 

sa 66-6 0-760 | 7:34 | 1-83 | 8-01 | 0-856 | 8-14 | 0-283 calc. 

3 anlents5 200 (a7o8 a8) MOO 2 era 7e7S 0-212 exp. 

| 62:1 0-608 | 7:34 | 1:50 | 7:80 | 0-663 | 7-87 | 0-218 calc. 
eee ees eee ee 
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Table 3. Dielectric properties of heavy water 


én) 
T(t) € ds N= NAT Cron A=3:°213 cm. A=10:-0 Se 
=f n K n K n K 
4-75 2-66 6°82 2:84 — 12525 vexp: 
: 85-8 3°84 4-54 2:66 6°84 2:89 — 1252) icale: 
ul 2-76 7:24 2:66 — 1-2 aexp 
- 83°8 Sel, 4-91 2:78 Tuy 2:69 — 1-26 Neale: 
a 5-65 2:84 HoTks} 2:26 — 0:94 exp 
ag 80:1 eo 5-56 2°83 LD LDS) — 0-93 cale. 
6:25 2 70 7-98 1°85 — Os 7iaexp 
a 76°5 1-76 6-09 MET eo 1-88 a 0:70. calc. 
_ On 2:54 7-98 | 1-47 — 0:54 exp. 
a 39 1:36 6-60 2:56 8-04 | 1-49 = Oz5Saneales 
a= — 7:98 28 | — 0-43 Pes 
50 
69-8 1-11 — — 8-02 1-24 | — 0-425 calc. 
¥ 2 = = 8-00 | 1-035 | x, 0-34 exp. 
60 : = 
66:7 0-923 = = 7-94 | 1-025 | — 0-343 calc. 


a SS ee ee eee 


agreement to + 1% with the values calculated assuming a single time of relaxation. 
The results were analysed using the value of « at 10 cm. and 1-274 cm. as a basis for 
finding 4,, the best fit being subsequently obtained by trial and error. At low 
temperatures it is useful to plot 1/[1 +A,?/A?] against ¢’ for different frequencies, to 
obtain straight lines passing through e, and ¢, when 1/[(1+AZ2/A7]=1 and 0. In 
this way ¢, is found to be 5-5 + 1, and the single time of relaxation is conveniently 
displayed; figure 6 shows straight lines obtained at 0° and 20°c. for water. The 
value of 5-5 may seem at first to be high, but is in agreement with the recent results 
of Saxton and Lane (1945), and with measurements of the refractive index in the 


far infra-red which show a value of n =2, rising as the frequency decreases (Rubens 
_ 1915, Cartwright 1936}. 


§9. PREVIOUS. MEASUREMENTS 


There is a wide field of literature on the dielectric properties of water at high 
frequencies. Although the measurements of individual workers are in many cases 
self-consistent there are wide discrepancies between the results of different 
authors. The early results with damped sources being generally regarded as 


The dielectric properties of water and heavy water 155 


Single point measured 
in cylindrical waveguide | 


Figure 6. 


Relaxation-time analysis. 


0 20 


Figure 7. 


40 0 
C in degrees centigrade > 


80 


Absorption coefficient H of 


water at 3:213 cm., measured direct- 
ly in rectangular waveguide, plotted 


against temperature. 


Syuares in- 


dicate accuracy claimed, + 1% and 


£0°4° c. 


unreliable, only those obtained with continuous wave sources are shown in 


table 4. 


Date > 


1937 
1937 
1937 


1939 
1939 
1940 
1943 
1947 


1946 
1947 


Table 4. Values of A, for water at 20°c. 


Authors 


A. Esau and G. Baz 
H. W. Knerr 
W. Hackel and M. Wien 


Slevogt 

W. Kebbel 

Divilkowski and 
Masch 

C. P. Connor and W. P. 
Smyth 

eeeA ver oaxtonmancdan): 
Lane 

P. Abadie and J. Girard 

Present authors 


Method 


Free wave 

Wire wave 

Thermometric at 
metric wavelengths 

Wire wave 

Free wave 

Thermometric method 


Wire wave 
Free wave 
Wire wave 


Waveguides and 
resonators 


Approx. value 
of A, 
(cm.) 
1-85 

1-6-1°8 
1:60 


For bibliographical references see end of paper. 


These results show that measurements in this field are peculiarly susceptible 


to concealed systematic errors. 
(i) The assumption of a square law voltage current law for crystals. 


The following may be cited as typical :— 
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(ii) Physically attractive but mathematically invalid approximations. ‘The 
so called ‘‘ wire wave” methods are open to such approximations as can be 
seen from the rigorous treatment which has been given by Slevogt. 


(iii) he relationship 
m+n? +1—2n 
4241420 


between the reflection coefficient R and the optical constants, shows R to be 
insensitive to changes in 7 and methods based on the measurement of R are 
correspondingly uncertain. 


(iv) Owing to lack of power, experimenters using the optical or free wave 
method have sometimes had to make measurements before the region of 
steady exponential attenuation had been reached. 


The authors hope that they have avoided or overcome most of these difficulties. 
As has been already mentioned all power measurements have been made with a 
piston attenuator, so that the receiver acts only as a monitor. Where possible 
independent methods have been used to cross-check some of the results. 

The most extensive series of independent measurements was made atA= 10cm. 
at which wavelength the E,,) resonator method was compared with both versions 
of the two-waveguide method over a large temperature range. The results given 
in table 5 show satisfactory agreement. AtA=3-2cm. « was measured using both 
cylindrical and rectangular guides. Reference to figure 7 shows that the two 
methods agree. Up to the present no method of cross-checking the 1:25 cm. 
results has been possible apart from doubtful tests of self-consistency and con- 
sistency with the results at otherswavelengths. 


Table 5. 10cm. results showing agreement between various methods 


n n K K 
P(e Gs) n Rectangular Cylindrical | Rectangular Cylindrical 
Resonator waveguide waveguide waveguide waveguide 


0 8:99+0:09 9-03+0-05 9-03+0-05 | 1:3740-01 1:37+0-01 
10 — 8:89 8:89 0-98 0:982 
20 8-86 8-83 8-83 0-74 0-739 
30 8-68 8:78 8-78 0:56 0-558 
40 8:50 8°53 8°53 0-435 0-442 
50 8:28 — 8:28 0-345 0-350 
60 8-09 — 8-09 0-28 0-281 
70 7°87 — 7°78 0-23 — 
80 eH — — 0-20 — 
90 Heys — — On? — 

100 7°38 — — 0-15 — 


910. THEOREDICAT 


It is significant, though not unexpected, that water and heavy water, in common 
with many other chemically pure liquids, obey the Debye equations exactly anda 
have a single time of relaxation. From these data the values of dielectric constant 
and loss angle can be calculated for any frequency. For completion of the work, 
measurements at wavelengths of 0-5cm. and 1mm. would be desirable, but 
technique is not yet advanced sufficiently for these wavelengths to be available. 
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By such measurements, greater accuracy than +1 in the value of €, would be 
obtained. 

It is also significant that the single time of relaxation should hold at all temper- 
atures, indicating that the change between the various structures suggested by 
Bernal and Fowler (1933) is continuous. It would be difficult to imagine dipoles in 
two different crystalline structures simultaneously present at any given temperature 
taking the same time of relaxation. It is apparent that the absorption lines shown 
by water vapour in this region have no effect on the dielectric constant of the liquid. 

The expression €—€ )=(e,—€,)/(1+jw7) which reproduces the experimental 
results was originally obtained by considering a dielectric in which the internal 
field was a linear function of the polarization. It is well known that this assump- 
tion over-estimates the interaction between the dipoles so that the dielectric model 
becomes spontaneously polarized except at very high temperatures. This. 
difficulty is considerably enhanced by increasing the value of e9. Moreover, there 
is a discrepancy between Debye’s 7 and the experimental 7 which arises from the 
assumption of the Lorentz internal field (see Polar Molecules, p.92). tis therefore 
important to consider how the dielectric constant would depend upon frequency 
on theories such as those of Onsager or Bottcher (1942), which avoid these diff- 
culties by making the interaction between the dipoles depend on the dielectric 
constant. The following treatment of Onsager’s theory follows closely that of 
Fowler and Guggenheim (Statistical Thermodynamics, 1939) except that an 
alternating field and a relaxation time are introduced. 

Consider first a molecule radius a, with a polarizability y, reduced to a refractive 
index 1 by the Lorentz-Lorenz formula as follows :— 


V= Orie). 9 | eae (1) 
If the electric moment of the molecule is , and it is immersed in a continuous. 
medium of dielectric constant ¢,, then by Onsager’s theory there will be a 
‘“‘reaction””’ field R given by 
R= pw 2legel)i(Ze, Nae.) 0) | hae ates (2) 
Where by definition w=,u+yR and u is the unit vector in the direction of the 
moment. Therefore 


=u by 2, 1)/(Zeseija ee (3) 
Combining (3) with (1) and re-arranging, 
= p(n? +2)(2e,4+ 1)/3(2e, Fm"), © canteen (4) 
and combining (4) with (2) 
R202 (gus Zee a. (5) 


Now consider the application of an alternating field E of frequency w, and let the 

dielectric constant corresponding to this frequency bee. Let the additional field 

this causes at the cavity be E’. 

Then by the Onsager theory 

Be ee Eye) 

~ 2e+1 (2e + 1)a* 

and by a similar argument as was used for equation (4) 
ee Ee 2yei(Qenn*). aes (7) 


E’ 
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The total electric force F acting on the molecule is given by 


FEE +R) ~ = eee (8) 
En’ +2je , Ute ee (9) 
so that F= Det ne 3(2e, + na Ko 


One may now follow the conventional treatment as given, for example, in Fowler 
and Guggenheim. In evaluating the mean polarization, the effect of the relaxation 
time can be introduced as in Debye’s original treatment by writing /3R M1 +jor7) 
instead of the expression ji9/3R7 appropriate to the static case. With this modifi- 
cation one obtains after some reduction the result that 


in which A is written for the expression 
(n?4+2)?.(2e,+1) 4 
(2e,+7).6RT "°° 
Since] /2¢is small comparedwith unity no great errorwill be involved by substituting 
e=e,/(1+jwr). 
Thus” «— n= A/(1+jwr) — A/2«,. 


As w-> 0,€,—(n?—A/2e,)> A. Ife, is defined so that A =e, —¢,, thenn?—4 = £9: 
We have, therefore, «—ey=(e,—€,)/(1+jwr), with éo+tt=n". The effective 
deviation is of the order of 1/4<2 and will therefore be vanishingly small. 

We can therefore say that it is impossible to differentiate experimentally 
between a linear internal field and an Onsager field. 

The above calculation also shows that the dielectric constant differs from the 


value calculated from the relaxation time curve by } as w— co on the Onsager 
‘model. 


It is just possible this might be found experimentally, but small deviations of 
this order could be attributed to many factors. 

It can definitely be stated that the experimental value of <q for water is higher 
than the currently accepted value of 2:0. It is interesting to note that if we sub- 
stitute €)=5-5 in Onsager’s equation instead of €9 = 2-0, the calculated value of the 
dipole moment is found by to be 1:5 Debye units. 

The accepted value of jy is 1-87 in the vapour phase, and it is therefore apparent 
that if the high value for «, is accepted, Onsager’s formula gives, unexpectedly, a 
value of ¢,= 130 which is greater than that found experimentally. 

However, this large value for €9 is not without its own difficulties; the value of 
<o for ice is 3, so that one might expect a considerable decrease in dielectric constant 
when water freezes. Also, if one substitutes the high value of €y in equations (3) 
and (9), one finds a ratio p/u)=2-6. Since Mo 1s already 1-87 Debye units this 
would mean so great an electronic polarization that y could no longer be regarded as 
constant. 

The only possibility is to leave pu throughout the equations, which would finally 
give 
cant te ae ee 

V Ze, 1 BRT +jwr) 2e+1 
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It is possible that the anomalous dielectric properties of “associated ” liquids 
are due to the fact that the atomic and electronic polarizations are very much higher 
than is generally assumed. It can be seen qualitatively that association might 
produce high polarizabilities. 

It is of interest to test the Debye relation 7 =82na?/2kT experimentally. 
We find for water and heavy water a fair approximation to the proportionality of 7 
and 7/7, which is shown in figure 8. It appears that a power law slightly different 


Sif T T } 


| D,0 
| | H,0 
| 


1/T «10-° poises per degree 


Figure 8. Relation between relaxation time and viscosity. 


from unity is obeyed; logarithmic plotting shows it to be a 0-97 law, from which 
extrapolation will give the real relaxation time at any temperature for which 7 is 
shown (e.g. for super-cooled water). Calculation of 7 on the accepted basis of 
a=1-38a. givesa value of r=0-85 x 10-11, in good agreement with the experimental 
value of 0-9 x 10-". 

Further, the ratio of viscosities and of relaxation wavelengths for water and 
heavy water show fair agreement over a temperature range, as can be seen from 
table 6 (viscosities from Lewis 1933). 


Table 6 
Smocthed values of : 
\sp,0/ 81,0 Ip,o' 14,0 
TC.) +2% 

10 1:30 1:29 
20 7 1-25 
30 1:24 2? 
40 1-21 1:19% 


* Extrapolated 


The close numerical agreement between the observed and calculated values of 
is surprising, and probably fortuitous, but the close dependence of 7 upon viscosity, 
both in water and heavy water, seems to us significant. It is evident that viscous 
flow and dipole orientation are influenced by the same factors, and can thus, in 
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some ways, be said to favour a hole theory such as that of Eyring (Glasstone, 
Laidler and Eyring 1941). If one accepts Eyring’s result 


eRe re 
— =, edSih  ¢ SHiRe 
is 


connecting the thermodynamic properties of the activation process involved in 
re-orientation, one can calculate AH from the observed temperature coefficient of r. 
The value obtained is about 5 kcal/gm. mol. at 0°c., which is roughly the same as 
the surface energy. The observed value of AH also falls with rising temperature 
as the surface energy is known to do. 
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ABSTRACT. An attempt is made to calculate the heat of solution of gold in silver on 
the basis of the quantum theory of metals. These two metals are chosen because they have 
the same atomic volume, and therefore are the simplest case. 

The steps in the argument are as follows : Suppose that a gold atom replaces a silver 
atom in the lattice. Then, to a certain approximation, one can represent the substitution 
of the silver ion by a gold ion, in its effect on the electrons, by a “ potential hole ”’ of depth 
AE and radius 79. This potential hole will alter the energy of the conduction electrons. 
To a first approximation the change in energy is just 4E, which would give zero heat of 
mixing. A second-order term of order (4E)?/Ey always gives a positive heat of mixing ; 
Ey is here the Fermienergy. This term is calculated exactly by wave-mechanical methods ; 
it gives 0-69 ev. per atom. ‘The same calculation shows, however, that there is a concen- 
tration of charge in the gold atom in excess of that in the surrounding silver atoms ; this 
alters the potential in which the electrons move, so that a self-consistent calculation is 
required to obtain the true energy. For this the labour required would be almost prohibitive ; 
therefore we use instead the Thomas-Fermi method and obtain 0:45 ev. We thus find 
0-15 ev. per atom for the heat of solution, which compares well with the observed value 
‘0-13 ev. 

With the help of the potential obtained with the Thomas-Fermi method the residual 
resistance of gold in silver is found to be 0:16 micro ohm cm. for 1% solution. The 
considerable discrepancy as compared with the experimental value 0:38 seems closely 
connected with very similar discrepancies found in other theoretical work on temperature 
resistance of the noble metals. 


Si INTRODUCTION 


N many applications of the electron theory of metals to alloys, it has been found 

possible to treat the different constituent atoms merely as sources supplying 

different numbers of metallic electrons which belong to the metal as a whole. 
The irregularities created by the foreign atoms have been entirely ignored and the 
electronic band structure has been treated in the same way as in pure metals 
An exception is the work on residual resistance of alloys, for the residual resistance 
is a direct result of the local irregularities in the potential field caused by the 
difference in the nature of the constituent atoms. However, most of this work is 
essentially formal and succeeds only in correlating certain facts, and leaving the 
question of the nature of the irregularities untouched. ‘The only attempts at 
quantitative estimation of absolute magnitudes of residual resistance seem to be 
those made by Mott (1936). An attempt to obtain some idea of cohesion of alloys 
with the simple model employed by Mott in his work on resistance has naturally 
led us to a closer examination of the irregularities in the potential surrounding the 
foreign atoms. It is the purpose of the present paper to investigate their nature 
more deeply and give a semi-quantitative discussion of the cohesive energy and 
resistance. For our actual discussion, the alloy Ag—Au will be used, because the 
constituent atoms have practically the same radii, so that the effects due to distortion 
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of the lattice will not be present. Similar effects to those discussed here will occur 
also in Cu-Ag, Cu—Au, but will be complicated by effects of distortion both in the 
consideration of resistance and energy. 

The sum of the sublimation energy and the ionization energy of the atoms in a 
metal will be called the binding energy of the metal; this is the energy necessary 
to separate the individual ionic cores and valence electrons into free ions and 
electrons. Then the binding energy of a monovalent pure metal with “ almost- 
can be put rather crudely as 


22 
w=-3(E+ 52), ee 


free electrons 


< 2m 


E° is aconstant for the metal concerned, and k, is the wave number vector of the 7th 
electron, so that the second term alone is just the Fermi energy of the electrons. 
Formula (1) can be derived theoretically (e.g. Seitz 1940, p. 345) when, owing to 
correlation, only one electron is in each atomic sphere. E°+ H?k?/2m is thus the 
eigenvalue of the zth state in the Hartree sense. (By adding these eigenvalues 
the mutual electronic repulsive energy is counted twice, and this serves very 
conveniently to take care of the repulsion between the ions). ‘Thus E£° appears 
as the constant potential energy of the Sommerfeld model. Mott (1936) has 
actually shown that to the same approximation this simple interpretation can be 
extended to alloys. ‘Thus the difference in the ionic field, say, of a gold atom ina 
silver lattice can be regarded as a change in the potential energy from E%, to E,, 
in the atomic sphere surrounding the gold ion. This idea is of importance in 
making possible at least a semi-quantitative treatment of effects involving such 
irregularities in the potential field. However, simple consideration will show that 
it is necessary to elaborate the square-hole potential, before it can be used as the 
effective field in which the electrons move. 

Consider, for instance, an Au atom in an Ag lattice. As the atomic radii are 
the same, the Fermi energies are the same. From (1) it follows that 

ER, = ye = W an a W se = 2:60 CVs) ee ee (2) 
where the values of the binding energy are obtained from the sublimation energy 
and the atomic ionization energy (Landolt Bornstein 1936). If the square-hole 
potential 
U(r) = EX, 7 >1)(=1-59-x 10cm!) 
= EX, — 2-60 Tat 

is taken as the effective field in which the individual electrons move, where the Au 
atom is supposed to be situated at the origin and 7» is the atomic radius, the eventual 
distribution of electrons can be obtained very simply by the Thomas-Fermi 


method with the particles treated as non-interacting. The number of electrons in 
any volume element d7 is given by rite: 


pdr =(31/3h*)[2m(E,—Ur)P2dr, a (3) 


where E,, is the energy of the highest level filled. Let p, be the density of electrons 
in pure Ag. ‘Then if the metal is taken large enough, obviously 


Po =(87/3h3)[2m(E,, — £3,)]?2, 


that is, the constant density outside the potential hole is identical with that of the 
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pure Ag metal and £,,— EF, is equal to its maximum Fermi energy:<. 5:5:6vs 
From (3) it follows that the ratio of the density for 7 <7q in excess of py tO py is 
given by 

(5:5 + 2-6)? — (5-5)3? 

(5-5)82 

showing a considerable accumulation of charge at the Auatom. The electrostatic 
potential produced by the excess of charge can be calculated directly. Acurve(C) 
giving the potential energy of an electron caused by this potential together with a 
curve (A) showing (p—pp)/pp is-plotted in figure 1. The energy is seen to rise 
above 10-0 ev. at the centre of the Au atom. The use of the simple square-hole 
potential alone as the field in which the electrons move, therefore, does not lead toa 
self-consistent result. 


= 0-78, 
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Figure 1. 


A= Ap/py by Thomas-Fermi method. 

B= Ap/py) by wave mechanics. 

C= Potential energy due to 4p by Thomas-Fermi method. 
D=Potential energy due to 4p by wave mechanics. 


Qualitatively it is clear that when the electronic interaction is properly con- 
sidered, the potential-hole will be negatively charged due to the preferential 
accumulation of electrons, and its immediate neighbourhood will be positively 
charged to screen off the field produced by the negative centre. The rather 
intricate nature of this charge arrangement suggests the method of Thomas-Fermi 
(with interacting particles) as the suitable method to treat the problem, at least as a 
first approximation. ‘T’o avoid excessive complication the exchange and correla- 
tional effects will not be considered, so that the Thomas-Fermi method can be used 
in its original form. ‘This procedure amounts to ignoring the change in the 
exchange and correlational effects brought about by the change in electron density. 

Since it is difficult to predict the accuracy of the Thomas-Fermi method on 
a priori grounds, the next section will be devoted to finding exact wave functions 
for electrons in the potential hole, and the results obtained will be compared, 
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Further developments in §3 concerning energy changes of these cases will be 
important for later discussions of the energy of cohesion in §5 
§ 2. 


EXACT WAVE FUNCTIONS AND CHARGE DISTRIBUTION FOR 
NON-INTERACTING ELECTRONS 

For the present purpose it is convenient to suppose the metal to be in the form 

of a large sphere of silver of radius R with a single gold atom situated at its centre. 

‘The normalized wave functions of the electrons have the general form 


G,(r 
Nua : ) ©,,,(66), 


where ®,, (04) is the normalized spherical harmonic for a state with orbital quantum 
number / and azimuthal quantum number m. 


N,, is the normalization factor for 
G,, which satisfies the radial wave equation (e.g. Mott and Massey 1933, p. 22) 


d?G;, i+1 
Cn +(#t= at a) G,;=0 TT 
dr 
“G; Ui+1) 
7a +(# = us) Gye eis opodor (5) 
hk’ h?k? 
where oes =HE-E,, oor =E-E, 


denote respectively the kinetic energies in the two regions 
E,, E, will be written in place of E2,, and EX, 


Here, for simplicity, 
The solution of (5) that is regular at the origin is 


Gyr) =AV Jp 522B'7) r<1o 


= BV Fs) kr) + CV 1 _asapy(kr) TAT 9; essa (6) 
Jj 44/2) J_a41/2) are the Bessel functions of half odd-integer orders 


. Since G,, and 
dG,,/dr are continuous at 7=7, the ratio of the constants A, B, C, is fixed so that 
only a multiplying constant is left arbitrary 
asymptotic form at infinity 


To define this, let G,,(7) have the 


De 
G40)~(5) sin (kr — 412 + my). 


corte (7) 
Then if the radius of the sphere R is taken sufficiently large 
DRE 
(Ni) 7 = al o (kr —4la + y)dr 
mess 
atk’ 
which determines N,, 


For the state represented by (4) to be a stationary state, the condition 


G(R) =0 
has to be fulfilled. When & is large (7) may be used. Thus 


Heat of solution and residual resistance of gold in silver 165 


where » is an integer. From (9) it follows that the number of states with definite 
values of /, m, and with k in the interval k to k+dk, is given by 


2R 
alo ea ake 2 ee (10) 


The factor 2 allows for the two states of spin. With (8) and (10) the electron 
density can be written down directly in terms of G7). Since the solution 
corresponding to G;,, where the gold atom is replaced by a silver atom is just. 
V/1J).1)(kr), which satisfies (7), it is seen that the change in density is given by 


m Kin Ohi irr, f2 - 
Xp=28{ 3 a6 %4),4(04)} | "(GO — hall) pay 
t W=—m YO 
ae he hm (Gat) — Thiol) 
= 5 E(2+1)] (Ce ee 0 ee (11) 


k,, is the wave number corresponding to the highest filled state. Owing to the 
complicated natures of the coefficients A, B, C, involved in G,,, the integral can 
only be evaluated numerically. The result is shown in figure 1, curve B. The 
Coulomb energy of an electron due to this excess charge density has been worked 
out from this curve by numerical methods and the result is represented by curve D: 
in the same figure. Owing to the discontinuous nature of the density function as 
given by the Thomas-Fermi method, the general agreement between A and B is 
as good as can be expected. ‘The potential is not discontinuous, and the agree- 
ment between C and D is remarkably close. Since the case treated actually 
involves a discontinuous drop in the potential, and since nevertheless the simple 
Thomas-Fermi method still yields reasonably good results, it seems reasonable to. 
believe that the method will in general lead to fairly reliable conclusions. 


S27) BS NER GY] CHANGES IN° THE CASE, OF HYPOTHETICAL 
NON-INTERACTING ELECTRONS 

It is well known that the Thomas-Fermi method fails generally to take account 
of the kinetic energy associated with rapid variation in electron density and gives 
too low a value for the energy when rapid changes in potential are involved. As. 
the square hole potential will still be used in the discussion on cohesive energy, 
it is important to be able to estimate the inaccuracy caused by its discontinuous 
change in potential, and to correct for it. ‘he most direct, and probably the best, 
course is to calculate the energy change caused by the presence of the Au atom in 
the hypothetical case that has been treated by both methods and use the difference 
to correct for the inaccuracy in later discussions. 

In the Thomas-Fermi method the average energy of an electron in any 
particular volume element is taken to be the sum of the potential energy at this 
point and 3/5 the maximum kinetic energy. Asp, corresponds to just one electron 
per atom, it follows from §1 that the number of electrons in the gold sphere is. 
given by (8-1/5-5)??, and the total energy of these electrons is 


(Fs) ( oo +2). or (12) 


Of these, for a pure metal, one electron is normally in the same sphere and the 
others are to be regarded as drawn from the top of the Fermi distribution of the 
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whole metal, so that the other electrons might be regarded as not perturbed by the 
presence of the Au atom. Therefore the energy of the “same” electrons, before 
the Au atom ion replaces the Ag ton, is 


3 8-1 3/2 8:1 3/2 
3 55+ (F3] -1|x55+(35) <ho  eeee (13) 


Subtracting (13) from (12), and using the value E,— Fy = 2-60 ev. as given by (2), 
the following value for the energy change is obtained 
—3-59 ev. 
A more significant quantity is given by 
= 3-59 -2(, — E,) ==0-99 ev.) eee (14) 
which is the energy change accompanying the transfer of one gold atom from 
bulk gold to pure silver in this hypothetical case. - 
The energy change corresponding to the rigorous treatment in §2 can be 

expressed in terms of the phase shifts 7, introduced in (7). We recall that the 
solution corresponding to G,,(7), before the Au atom replaces the Ag at the origin 
is \/rJ 141 0(kr), so that the phase shifts are originally zero. It then follows from (9) 
that if R is taken to be very large, the value of k corresponding to the same quantum 
number 1 is altered by the infinitesimal amount 
dk= —1,,/R, 

when the Au atom is introduced. Accordingly, the change in energy is given by 

(8) 1h nw, 
2m mR 

which is independent of the azimuthal quantum number, so that on summing 
over all the filled states with the help of the density of states given by (10), the total 
energy change is obtained in the form of a sum of integrals 


he km 
4 (=) 2 2(21+ 1) i, Riqydhi-. es oeke eee (15) 


In connection with the evaluation of the expression (15), a simple and rather 
interesting identity may be deduced, which makes it possible to obtain a result of 
the desired degree of accuracy without numerical work involving terms with />2, 
so that the necessary values of 7,,; can be taken directly from results obtained in 
working out the last section. 

It is well known from the theory of collisions that an approximation to Nr 1S 
given by (Mott and Massey 1933, p. 28) | 


amt (~ 
ne Fal, UV Oy getkr 27 ree (16) 
where U(r) is the perturbing potential which vanishes at infinity, so that 
Un= —AE= —(E,-E,) T<Np sass oe le (17) 
=() Tf, 


in the present case. Substituting (16), (17) into (15), the following approximate — 
value for the change of energy is obtained: . 


km To 
—AEE2(21+1)[ "kdk | (J,saio(kr))?r dr. 
aaa 0 
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On the other hand, on using the expansion (Watson 1944) 


a \l2 © 
a (a) & it(21+ 1) P0086) FesslRr), 
it is easily established after integration over 6 that 


3 
aah 
‘o) 


=| ew thr x ek dy 
) atomic sphere 


77 iio 
=F PAA+Y| rraylhePar, | 
On substituting (19) in (18), the following identity is obtained: 
r hin ‘al ar 
25.455 UAB | dk) (Iuyolkr))°r dr | 
J0 J0 | x 
ns f teadies (20) 
=~ ae (5 )(kyroy = —AF, J 
where use is made of the fact that the distribution corresponds to a density of one 
electron per atomic sphere. (20) states the fact that when the approximation (16) 
for 7.18 used, the energy change obtained is just — AZ, the value that would have 
been obtained if the first-order perturbation theory had been used. 
We use this result in the following way. We desire to find the change in 
energy using exact values of 749, 74, only. With the help of (20), the change of 
energy can then be written as 


kin To Km Pry 
Be Nee NB 2. Rdk| (Jyjo(kr))2r dr +6 kdk|  (Jg(kr))2r dr 
a) 0 0 a 


0 


f2 km phem 
E 52 hypoth +6 | kuadh | Oe tS ote) (21) 
AVR oy 


which corresponds to (15), but with values of y,, given by (16) for J>0. The 
integrals in the second term can be expressed in terms of the Si functions (Jahnke 
and Emde 1945) and the integrals in the third term are evaluated by numerical 
integrations. ‘The result is 
—AE+AE(0-957) — 3-18 = — 3-29 ev. 

The first two terms together give the total contribution of all states with />2. 
It has the very small value 0-llev.; the assumption that deviations from. (16) 
could be ignored for states with />2 is therefore justified. The more significant 
quantity corresponding to (14) of the Thomas-Fermi method (the heat of mixing) 


is in this case 
—3:29+ABF=—0-69eV. 4 waceee (22) 
Comparison of (22) with (14) shows that the energy value given by the Thomas- 
Fermi method is too low by 0-30ev. Although the discrepancy is not very great, 
considering the abrupt change in the potential, it is, however, clear that the 
conclusions regarding cohesive energy developed in later sections cannot be 
expected to be more than semi-quantitative. For the ultimate quantity of 
interest (heat of mixing) corresponding in the real case to (22) and (14) in the above 
hypothetical case should be of the order 0-1—0-2ev. So the correction of 
0:30 ev. is by comparison very considerable. 
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$4. APPLICATION OF THE THOMAS-FERMI METHOD 
TO THE ACTUAL CASE 

‘The considerations of § 1 show that it is necessary to take electronic interaction 
properly into account to obtain a self-consistent result. The square-hole potential 
describes the difference of the ionic field within the atomic sphere of the gold atom 
from that of the rest; so that if the electron cloud were uniform, it would be the 
effective potential in which the electrons move. ‘Therefore, in order to treat the 
actual case, a natural procedure would be to take the potential as the superposition 
of the square-hole potential and that produced by a uniform space distribution of 
positive charge that would neutralize the electronic charge completely if the latter 
were uniformly distributed. Let (Ur) be the square hole potential as already given 
by (17) and V(r), —e be respectively the electrostatic potential and the electronic 
charge. The density of electrons at any point is then (87/3h*)[2m(E,, — U + Ve)]?*, 
which leads to the following Poisson’s equation for V: “| 


3 3e 


: 2ne aoe 
V?V = —47p= ay 8 [2m(E,, -U+eV)]* air 39. Sea (23) 


m 


for the positive charge density is e(477,?/3)-'._ Equation (23) can be simplified by 
introducing the dimensionless quantities 


og =(E,, — U(r) + eV(r))/Ex, 
5 ( ) ( ) | ( ))/ | is EIA Pte (24) 
E=7)/(roEy/3e?)"?7o, J 
Omr\2/3 ~h2 i Hl 
where Ey= (FZ ) 2mree | 


is the maximum kinetic energy of a free electron gas of the same average density. 
In terms of (24), (23) becomes 


Ind fe des , 
B AG =) =p 9 eee (25) 
Introducing f defined by 
FQ=o0), = ee (26) 


a Z [(Z y “2 1] OS See (27) 


Since the change in electron density is small compared with the average density, — | 
it is advantageous to introduce A(é) defined by | 


KO=E 4 AC) se (28) |i 

In all cases treated, A/é can be treated as small and it is sufficient to solve (27) to | 

the second approximation. We thus substitute (28) in (27) and expand /3? in — 
powers of A/é to the third term. ‘The result is 
HN ae cS 

de 2 se Oe ave (29) 

Before solving (29) it is useful first to examine the boundary conditions to be_ |} 

satisfied by A. As A must tend to zero as € tends to infinity, it is clear from (29)% 

that A must have the asymptotic form 
Arde 324) Sot Se eee (30) 


where A is an arbitrary constant. 


(25) reduces to 
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At the origin, as the charge density must be finite, it follows directly from 
Gauss’s law that 


dV 
chad 0 as r—0, 
dg 
and hence that —2. —>(), 
dr 


Therefore it follows that 


dX 
A> €— and ONS =. 


dé dé 
Since g—1 must remain finite, it is necessary that 
A(é)+ 0 as Ce Oe ou eatets (31) 
Therefore in the solution at the origin only one arbitrary constant remains, 
defined by 
cee) 
Ge) 25 


which, together with the arbitrary constant A in the asymptotic solution at infinity, 
is to be determined by the boundary conditions on the surface of the atomic 
sphere €=&). On this surface the electrostatic potential and field must be 
continuous. It follows that 


Aa(Ex) ~ Bal) = El gaa) ~gu(€)) =f EPH) — 6, Ee, 
(e), ~(Gz), = eles) —axtéo) + fo( (2) — (2), ) 
= U,(&) = U3(6) = i _ aes (32) 


where the suffixes 1 and 2 indicate the limiting value for €<& and €>€, 
respectively. (32) will determine A and B, which are left arbitrary in (30) and (31). 
Let 4(€) be a solution of the homogeneous part of equation (29), so that 


d*h 3 
Te i ne (33) 
Construct the integral 
Ss d2A d*d 
— —A—5 ) dé. 
[, (6a “aa 


A simple application of Green’s theorem leads, with the help of (29) and (3), to 
Es 3 5 A? 
ie —— 1c aan) OA 3 
| tp 4% (34) 


the following relation: 
d dA ay 4 = 
dé dé}|:, 8 


which makes it possible to give explicitly the value of A and dA/dé at any value of , 
both for >, and €<&,, in terms of the constants A, B in (30) and (31) and the 
integral involving the small quantity of second order A®/g*._ For the region € >), 
it is only necessary to take €,= 00, &;=€ and put ¢ equal to eV** and e-V 3/28 
successively. The limit on the left-hand side at €, is determined with the help of 
(30) and two linear equations are thus obtained with A and dA/dé as unknown 
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and may be readily solved. For €<€, the procedure is similar, the only difference 
being that €, is taken to be at the origin and (31) is used accordingly. A first 
approximation to A involving only the constant A or B may be obtained by leaving 
out the integral involving A2/é. When these first approximations are used for A 
in the integral, A and dA/dé are obtained in an explicit form to the second approxi- 
mation. ‘The following are these solutions: ; 


“ 9 


(tee ee: 
E>&, NE Jy) 20Ae s+ 


[eE1( — 3a€) —e-* Ex( — aé)]} | 
ne (35) 


aX*-1 fe Aza bass 4 
fe hee ap) —a§ 7 pag Say at Hy ( —o * 
de 1 2% Ae 3 [e%* Eu( —3x€) +e Ex 7E)I 6 i 
1 : Bel, pam 2 eee “at fe ox Sinh? x 
255 ee {2B sinh a€ + ak ie et : dee | e 5 av} 
gx Bay ei easing aes , sinh? x | |} 
See a 2 ve fit aa 2 p—aé x 6 
dé 5| 2B eoshse + aL hi if oes 
oeuvre (36) 
where for brevity « is used for V 3/2, and 
DD ee 
—Ei(—x)=| — dx. 
Ei( — x) er x 


(35) and (36) are not particularly convenient to give the general form of A, but they 
make it possible to determine directly A, B from the boundary conditions (32). 
Once they are determined, a single numerical integration of equation (27) provides 
the complete solution in numerical form and serves furthermore to check the 
accuracy of (35) and (36). The constants determined are 


A= 5:50, B=0-1122, 


ii 


The error in A, A’ is less than 1°% and is accordingly smaller in fandg. The total 
effective electronic potential energy as function of € calculated from the numerical 


solution is represented in figure 2. In the same figure is also shown the corre- 


sponding distribution of charge. It is seen that the qualitative feature of the 
charge distribution is just as expected, a negatively charged core surrounded by an 
equal ainount of positive charge in its immediate neighbourhood. The effect on 
the potential is to fill up the potential hole to a very considerable extent and should 


tend to reduce the electrical resistance; this will be discussed quantitatively in §6. — 


§5- HEAS OF SOLUTION OF A DILUTE SOMUDLON Ofc nNEAG 


The heat of solution is a convenient measure for the cohesion of analloy. F or | 


a dilute solution of Au in Ag, it is the energy change accompanying the transfer 
of one gold atom from bulk gold into bulk silver. First consider the energy 
change in the system, when an Ag ion in pure bulk Ag is replaced by an Au ion. 
The interaction of positive ions is not altered and may be left out of our consider- 
ation. If one adds up the energies of the electrons regarded as moving in the 
permanent field U— eV, one counts the Coulomb energy of the electron cloud twice 
over, so the energy of the system is obtained by subtracting from the sum of the 
electron energies the Coulomb energy of the electron cloud. Because of the 


a |) 
ty 
: 
: 
' 
Wy 
: 


spherical symmetry, for purpose of calculating change of energy it is only necessary | | 
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to consider the electrons in a finite sphere of radius R if R is taken large enough. 
Therefore the energy change is given by 


ig (3 g 2 
J =(2 Se, R* 74 
Bee elt G eV) 5 |v (5 =) | 


To 


lim 


R>x¥ 


eae —e Se (Re 7 fe 
all pods Br + Ey 5) | ac =) | pa ee (37) 


The first integral corresponds to the solution after the Au replaces the Ag ion and 
the second integral corresponds to the pure Ag. In both integrals, the three 
terms give respectively the kinetic energy, potential energy and half of the Coulomb 
energy due to the electron cloud alone. In (37) the Coulomb potential due to the 


oe : f\ : 7 
| 

: | 
{ 

1 
| 

u 


Potential Energy in cv. 


Deviation from Uniform Charge Distribution 


Figure 2. 
A=Deviation from uniform charge= Ap/pp. 
B= Potential energy curve U—eY. 
C= Original potential hole. 


electron cloud alone has been obtained by subtracting from the total Coulomb 
potential V, the potential produced by the uniform positive charge, namely 


3e(R_ 
ee NE 


Expressing (37) in terms of the dimensionless quantities g and €, and making 
use of the fact that the solution satisfying (30) corresponds to a neutral distribution 


R 
(mle tar =O, <1 ou, eevee (38) 
R>on- 0 
(37) becomes 


ee * 2 ge(& . U-a 3)- an 5) 
3Ex(=) lim [eaele (5+ OR, +5 Tot aahh? chee (39) 


R->m 
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which leads to 


3 “Rf é2 2 . g Ur & 
36o(E) im [tine + om” + (CGE) io ee +a) 


R->@ - 


on account of (25). The integral can be simplified by integration by parts. 
It is to be noted that both g and U are discontinuous at &,, so the integration has to 
be carried out in two sections. Then after a certain amount of manipulation with 
the help of the boundary conditions satisfied by g, most particularly those at ) as 
given by (32), the following expression is obtained : 


(E, — E2) 365 (BE, -£,)& Ga. oF al. = (4) dg — : , Soe 1) dé} : 
ne (40 


‘This expression gives the energy change when a gold atom is substituted for a 
silver atom. 

The most important feature of (40) is the fact that (£,—£,) is naturally 
separated out from the rest. It is the energy change that would be obtained if the 
electronic wave functions were not distorted by the potential hole; and at the 
same time, it is, according to (1), the energy difference between a metallic Au atom 
and a metallic Ag atom (i.e. average energy per atom in bulk material), so it leaves 
the rest of (40), namely 


_3J2 » (dg ip ed Ge Sie 
set Bve(E), + 70, OCH) 31, Se Dap,  e ae 
as the heat of formation per atom of Au in Ag. This is a very small quantity 
compared with (£,—£,), so the natural separation of \E, — E,) is of considerable 
significance in the actual evaluation of the heat of formation. The last integral 
in (41) contains an integrand which takes on both negative and positive values and 
nearly cancels out when integrated over the whole range. It is, however, possible to 
make use of (38) and obtain an approximate but more convenient expression for 
it of sufficient accuracy. If we write (38) as 


0= | ede(g2—1) 

~ [7 eae[ (143) Jf" [384 3(8) Jee 
Or 
[itera [(—) ede 3" (2) a= 3)" eas 


We introduce an error less than 1% in (41). Thus (41) becomes 


2 ay hte Le ef dee om 
=n en8(), + ao], PCG) +3), Cle Death 68) 
The integrand is now always positive and is convenient for evaluation, and further- 
more, it shows quite generally that the heat of solution must always be positive. 
The use of the numerical solution obtained in the last section gives the value 
0-45 ev. Applying the correction to take account of the kinetic energy associated 
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with the rapid change of potential of the potential hole obtained in §3, the value 
0-15 ev. is obtained. |The value obtained by extrapolation to infinite dilution 
of Wachter’s data (1932) from measurements made on electrochemical cells is 
O0-l13ev. The agreement is close. In view of the considerable correction that 
has to be made, however, the value 0-15 ev. should probably not be taken as more 
than an indication of the order of magnitude. 


§6. RESIDUAL RESISTANCE 


It is fairly straightforward to work out the residual resistance due to the dis- 
solved Au atoms from the numerical solution of the field obtained in § 4 and repre- 
sented in figure 2, The residual resistance is determined by the scattering of the 
electrons on the surface of the Fermi distribution (e.g. Mott and Jones 1936, 
p. 258). If the wave functions of such electrons are of the form given by (4), the 
equation satisfied by them is (Mott and Massey 1933, p. 22) 

a7G;, U(r) i+1) 

der)? 3c (1 Ey (rp ) Gr=0i) We eee (44) 
where U(r) is the potential energy of an electron moving in the effective field, the 
zero being chosen so that U(r) 0 away from the Auion. As U(r) is only known 
in numerical form in the present case, it is very much simpler to treat this collision 
problem by the method of phase shifts (Mott and Massey 1933). If the phase 
shifts 7, are defined in terms of the asymptotic expression for G, by (7), the scatter- 
ing cross-section into a solid angle dw is given by 


I(0)dw=|f(0)\2dw, aa (45) 


where /(@) is given by the general expansion 


1 0 
FO\==5— 2% (14127 = 1)Plcos8), anes (46) 
21R py, 1=0 


‘The resistivity due to Au ions in a dilute solution in which the ions may be regarded 
as scattering electrons independently of one another, is given in terms of J(@) by 
the formula (Mott 1936) 


hR, 
pa—eA ae (47) 
with A=\| (—cos6) 1(6)2nsin@ 26, 
Jo 


where x is the atomic percent. ofthe Auatoms. Substituting (45), (46) and (47) 
and integrating over the angle @ with the help of the well known properties of the 
Legendre polynomials, it is found that 
Ce at [(27+ 1) sin? n, —2/ sin sin 7, cos(m_y—7)]. «se (48) 

m l= 
Owing to the somewhat larger size of the scattering centre (figure 2) as compared 
with the potential hole treated in §2 and $3, contributions from /=2 term are also 
considered. 1, and 7, are obtained with the first approximation given by (16) and 
nm has been obtained with better approximate methods so that the error is of the 
order of 1°4,. The values obtained are 


19207189, —-7,=0-005,  y.—=—0-025. sss. (49) 
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In this case 7, is smaller than 7, for U is partly negative and partly positive and the 
maximum of G, lies somewhere in between, whereas G, has its maximum entirely 
in the region where U is negative. Using (49) in (48) and converting the result 
into practical units, the following resistivity for 1% of Au atoms is obtained: 


0-16 microohmcm. 


The agreement with the experimental value 0-38 micro ohm cm. is poor. It is 
particularly to be noticed that the close agreement of the potentials of the hypo- 
thetical case treated in §1 and §2 calculated with the Thomas-Fermi method and 
the rigorous method indicate clearly that the application of the statistical method | 
‘o obtain the field should be almost quantitatively correct. It is, however, believed 
that this discrepancy does not reflect any particular inadequacy of the model 
used for the dissolved ion. In the first place the resistance depends on t1e square 
of the depth of the potential hole, so it is a very sensitive quantity. Secondly, it 
only depends on a minute fraction of the bulk of the electrons in the metal which are 
most likely to be affected by the zone structure. Furthermore, it should be 
remarked that the only serious attempts at calculating theoretically the temperature 
dependent part of resistivity of the noble metals by Bardeen (1940) and Peterson 
and Nordheim (1937) show very similar discrepancies with experimental results 
(Bardeen’s result shows almost exactly the same proportional error as the result 
just obtained, and Peterson and Nordheim’s result is too small by a factor 2). It 
appears therefore reasonable to suppose that all these discrepancies are due to 
treating the electrons on the surface of the Fermi distribution as free. 


$7, CONCLUSION 

It is clear from the preceding considerations that when a foreign atom is 
dissolved in a metal, an appreciable polarization of the electron cloud will in general 
be caused in its immediate neighbourhood. Either a negatively charged centre is 
created, surrounded in its immediate neighbourhood by an equal amount of 
positive charge, or the reverse takes place. Any theory of cohesion must take this 
factinto account. Owing to the limitation of the method employed, the reasonat le 
agreement of the result obtained in this paper with the observed value could not be 
looked upon as more than semi-quantitative. It is to be appreciated that the small 
value of the heat of solution as compared with the primary perturbation in con- 
junction with the delicate question of polarization makes the problem difficult to 
treat with great accuracy. Finally, the polarization of the charge tends to fill up 
irregularities in the potential that would otherwise be there. This reduces the 
residual resistances caused by the presence of foreign atoms very considerably. 
In the case treated in the paper, namely, Au in Ag, this leads to a result rather at 
variance with experiments but at the same time brings it into line with other 
theoretical works on temperature resistance, indicating therefore a general 
inaccuracy in treatment rather than particular inadequacy connected with the 
model employed in this paper to represent a foreign atom. 
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ABSTRACT. Theoretical expressions are obtained for the momentum distribution 
of the electrons in atoms of different atomic number. ‘This enables density curves to be 
drawn for the momentum p. The mean momentum, which is also calculated, shows a 
steady increase with atomic number. Finally, the new calculations of momentum distribu- 
tion enable earlier work on the shape of the Compton profile to be extended from sodium 
to potassium, in monatomic form. 


Si, INTRODUCTION 


HE spatial distribution of electrons in an atom has been fully studied in 

recent years, and the use of variational methods and self-consistent-field 

modifications for solving the complete Schrodinger wave equation for 
an atom have enabled electron density curves to be drawn, in ample and satis- 
factory agreement with experiments using X-ray scattering. But on the other 
hand, the velocity distribution is relatively unknown. Yet it is, as Dirac (1947) 
has shown, equally fundamental. In addition, as Jauncey and DuMond 
(DuMond 1933) have explained, it is the velocity distribution which determines 
the shape of the Compton line in x-ray scattering; and, as Hughes has demon- 
strated both experimentally and theoretically (e.g. Hughes and Mann, 1938) 
this velocity distribution also describes the inelastic scattering of electrons. 
It has seemed desirable, therefore, to make explicit calculations of the velocity 
distribution in simple atoms (Z7<20); the results of such calculations are pre- 
sented in this present paper. . 


§2. METHOD OF CALCULATION 


Following the methods developed by the writers in a series of papers since 
1941 (see Duncanson 1943 for a complete list of references), we calculate the 
momentum p of an electron rather than its velocity. As we use atomic units 
throughout, the units of momentum and velocity for an electron have the same 
numerical value, viz.c/137. The distribution of p for a single electron is governed 
by a momentum wave function y(p); there is one y(p) for every type of atomic 
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orbital, and y(p) is best found from the normalized space wave function (r) 
by using a fundamental formula of the Dirac transformation theory: 


A BNE (2n)-92| exp (<iprid(rydry 2 ee (1) 


In this form y(p) is already normalized. Associated with y(p) is the mean 
radial distribution function /(p). This function is defined (Duncanson and 
Coulson 1945, referred to hereafter as D.C.) in such a way that I(p) dp is the 
probability that the electron has a momentum whose magnitude, independent 
of direction, lies between p and p+dp. Evidently 


W(p)=|x(P)x*(P)P"deoy, nee (2) 


where the integration is taken over all values of the solid angle w,. In this way 
I(p) is normalized so that 


pO 


| ; i(p)dp=1t . (3) 


‘The mean momentum # for this orbital is given by 


nS i ig ee ee ey (4) 


Our present calculations start with a presumed knowledge of #(r) for each orbit 
in the atom: (1)-(3) then allow us to determine y(p), I(p) and p for the separate 
orbits. Details of the analysis are given in § 3. 

Now, if there are m electrons in an atom with coordinates r,...r,, and momenta 
P.--P,, the total wave function y(r,...r,) which describes the complete electron 
configuration is related to the total momentum wave function X(p,.-.p,,) by 
a formula similar to (1) for a single electron: 


X(Py 2.8 °P,) = (2ny-3"?| 28 t(D : Ti ol teats +P, : r, Fr, oe .)dr”. 


But ‘F(r,...r,,) is the sum of one or more determinants (Slater 1929), and. so; 
on expansion, it may be regarded as the sum of a large number of distinct terms, 
each of which is the product of m atomic orbitals. Since the coordinates of each 
electron appear once and once only in each of these terms, the integration (5) 
for any one term is immediate; the result is simply that which would be obtained 
by replacing each space atomic wave function by the corresponding momentum 
wave function. By addition, therefore, it follows that the complete wave function 
X(p,..-P,) is formally precisely the same as V(r, .. .r,,) except that each space 
orbital is replaced by the corresponding momentum orbital calculated as in §3. 

Knowing X(p,...p,,) the calculation of the generalized mean radial distribution 
function I(p,...p,) follows at once from an integration completely analogous 
to (2). All the terms are known, so that I is determined for all atoms. 


§3. WAVE FUNCTIONS AND FORMULAE 


For hydrogen, with only one electron, an exact space wave function ¢4(r) is 
known. In atomic units it is simply (z)-”2 e. For atoms in the second row 
of the periodic table we have used the analytical wave functions originally proposed 
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by Morse, Young and Haurwitz (1935) and corrected by the present writers 
(1944). The resulting mean radial distribution functions have been given in 
D.C. (p. 193) and need not be reproduced here. For the remaining atoms 
(11<Z<20) we have used wave functions suggested by Slater (1930); that is 
to say, the complete wave function for an atom is made up in the usual way 


(Slater, 1929) as sums of determinants in each of which the electrons are defined 
by atomic orbitals of the form 


Brim (1) = NONG NG te Pe (cos Ayer, 


The N’s are separate normalizing factors for the radial and the two angular terms 
in %, n* is the effective principal quantum number whose value is given by Slater 
in terms of the true quantum number, cis a quantity related to the nuclear charge Z 
and the screening constant s by the formula 


In this way all the coefficients in (6) are supposed to be known. 

We have also used Slater functions for some of the lighter atoms (Z7<10). 
The results, which we do not need to describe in detail, since they are presumably 
less accurate than the Morse functions, show that for these atoms there is very 
little difference between the momentum distributions predicted by the two types 
of wave function. It seems highly probable, therefore, that the same conclusion 
will hold for the heavier atoms (11 <Z<20) for which the more detailed wave 
functions are not known. 

The transformation from the complete space wave function for all the electrons 
in an atom to the complete momentum wave function requires, first of all, that 
we should be able to transform the individual atomic orbitals v(r). This trans- 
formation from ¢b(r) to y(p) proceeds as follows. Using (1) and (6): 


x(P) = (2n)-82N,N,Ng|r”*1e-"P(cos Bete dr... (8) 
Now Bauer’s expansion of e~'P" gives 
e~'Pr = 5(2j + 1)(—2)'P, (cos 9) f,(pr), 


j=0 


where # is the angle between r and p, and the functions f, (pr) are defined by 


f(pr)=sin pripr, E(pr)=V/(m/2PNVpaplPr)- vases (9) 
Substitution of these values in (8) allows the integration over the solid angle 
dw [dr=dwdr] to be made at once. And if we write p,6,,4, for the polar co- 
ordinates of p, we have 
x(P) =(—2)'4/(2/7)N,NoNgP% (cos 6, )e#"*?R(p), ....0- (10) 
where 


Rp) = [ eh pdr a, (11) 
eto) 
The fact that x(p) depends on angle in precisely the same way as 7/(r) has already 


been noted, for hydrogen-like wave functions, by Pauling and Podolsky (1929) : 
our argument shows it to be true for wave functions in any central field. 
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Using (10) and the definition (2) of I(p) it follows that 


(p)=(2/m)N2PPR(PY, nee (12) 
where, by a simple integration it may be shown that 
N2= (26)? "+/P(Qn* 41). eed (13) 


The calculation of y(p) and I(p) for any one electron now depends only on a cal- 
culation of Rp). This must be made separately for s, p, d,... electrons for which 
1=0,1,2,... If n* is integral, as it is for the first three rows of the periodic 
table, the integration (11) yields algebraic functions. These functions are 


32c>p? 32c>p?( p? — 3c”)? 
L,(p) = preys ) I,,(p) = an c2)8 : ’ 
512c7p4 1024c°p*( pp? =<)? 
I,,(P) = er art I,.(D) a en ee z 
1024c7p4(5c? — p*)? 4096c%p® 
I;,(P) = ener ? I3(P) = Eee +) lug “Spa rgeces (14) 


‘Care must be exercised in this and later formulae to remember, from (7), that c 
has a different numerical value for each orbit in each atom, and the appropriate 
value must always be used in any particular case. Different values of c in any 
single formula, however, merely correspond to changes in scale of p. 

For electrons of higher quantum number than those in (14), 2* is not integral, 
and the evaluation of R,(p) is more clumsy. It may be systematized most simply 
in terms of two new sets of functions C(7) and S(z), defined by 


CG) — ii Tea COS Pi Ute aes S(n) = ig Pee Sin pr dre knees (15) 
‘C(n) and S(n) are functions of c and p, as well as of 7, and if we put 
p=ctana, (Cae a/2) 7 eee (16) 
we may write C(n) =T(n+1) cos (n+ 1)x cos”+z, 
9(2) = D(a +1) sin (2 1)a,. cos” 40 eee (eu) 


In terms of these functions it may be shown that 


PRo(P) = S(n*),- 
PPRy(p) = S(n* — 1) —pC(n*), 
PPRa(p)= 35S(n*—2)— 3pC(n* — 1) —p?*S(n*), 
p'R,(p) =158(@* — 3) —13pC(n* —2) — 0 Se — 1) Cir) ee eee (18) i] 
From these relations, with (10) and (12), both y(p) and I(p) may be calculated | 


for all electron orbits of s, p, d and f type. The simple results in (14) are, of 
course, merely particular cases of (18) and (12) in which n* is integral. 


$4. THE CAL CULATIONSSAN DS RESWilains 
As determined above I(p,...,) is a function of the magnitudes Pinte 
The simplest way of expressing it pictorially is to integrate over all values of 
p2.-+p, from 0 to o, obtaining a function I(p,), which we may call the mean 
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radial momentum density for electron 1. Now, by the very nature of the 
determinantal wave functions used (which results from the fact that all electrons 
are identical), I(p,)=I(p,)=.... I(p), say. So we may regard I(p), which 
1s normalized to unity, as measuring the velocity distribution in the atom, just 
as the more familiar density p(v) measures the space distribution. It is this 
quantity I(p) which we have calculated. 

If all the original atomic orbitals (6) are orthogonal, it follows that I(p) is 
simply the sum of contributions such as (12) from each electron present. In our 
calculations for the atoms of the first two rows of the periodic table, this condition 
is accurately fulfilled. But for the others it is not, and although orbitals with 
different / values are rigorously orthogonal on account of the surface harmonic 
terms in (6), there is a small non-orthogonality between 1s, 2s, 3s, 4s, and between 
2p, 3p, 4p etc. Calculations 
aC yp: 193) for ka and'C 
show that the result of 
neglecting this non-orthogon- 
ality (which is rather clumsy 
to include) is to change the 
width of the Compton profile 
(§5) by less than 1°/,. Now 
it would be presumptuous to 
claim an accuracy as good as 
this for the initial wave 
functions: accordingly, for 
the heavier elements for 
which 11<Z<20, we have 
made our calculations of I(p) 
on the assumption that all 
the atomic orbitals may be 
regarded as orthogonal. The 
curves of I(p), therefore, 
which we show in figures 1 
and 2 are simply sums of 
curves corresponding to (14) 
and (18), suitably weighted 
according to the numbers of 
each type of electron present. 
In every case, of course, 


r T T 


Fr Li 


Figure 1. 


[. I(p) dp=1. 


Certain deductions from these curves suggest themselves at once :— 

(a) First of all, there is considerable resemblance between the I(p) curves 
in figures 1 and 2 and corresponding electron density curves. For example, 
just as in the case of space wave functions, there are momentum shells within an 
atom. But the behaviour of these shells is exactly opposite to that with space 
functions. For as we proceed along any row of the periodic table, figures 1 and 2 
show that the shells expand (instead of contracting), and when we start a new shell, 
it appears at the inside (instead of the outside) of the ones already there. The 
writers have never seen this expressly statéd befote, but it is imiplicit in the 
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Fourier relationship of the two types of wave function; for it is a known result 
that a function and its transform cannot both be made dense in any given region: 
and this, of course, is merely the Heisenberg Uncertainty Principle in its. 
fundamental form. 

An alternative physically significant interpretation of the behaviour of each 
new shell is that the electrons in the new shell are at first loosely bound, and 
therefore have, on the average, 
small momenta. As the shell 
is completed, the binding in- 
creases, so that the mean 
momentum increases also. 

The steady drift to larger 
values of momentum with 
increasing nuclear charge is 
well illustrated by the be- 
haviour of the 1s electrons. 
According to (14) the peak of 
the 1s curve lies at p=c/4+/3. 
But c varies from 1 for hydro- 
gen to 18-7 for potassium, so 
that the peak moves steadily 
from 0-58 a.u. to 10-8a.u. At 
the same time the height of 
the peak falls, for the height 
per electron is 27/87c. When 
the contributions are weighted 
according to the number of 
electrons present, this means 
that the maximum contribu- 
tion from the 1s electrons falls 
from 1-05 for hydrogen to 0:01 
for potassium. After this stage 
we may say that the 1s elec- 
trons are scarcely relevant to 
the density function except for 
very large values of p. This is 
quite different from the be- 
haviour of the space wave 
functions, where (1s) crowds 
steadily in towards the nu- Fiareno 
cleus, and may eventually be 
regarded as completely screening the nucleus. 

(b) We also notice that s-electrons give uniformly narrower, or more compact, 
I(p) curves than the corresponding p-electrons. ‘This means that the alkali 
atoms for which the outer shell is simply s, and, even more so, the atoms of the 
next column of the periodic table, for which the outer shell is s?, show a pro-. 
nounced sharp peak in the I(p) curve. But progressive addition of p-electrons 
smoothes this out so that, for example, whereas in Li there is a very clear 
separation between the 2s and the 1s contributions, by the time that we have reached. 
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carbon the distinction is obliterated and one single wider I(p) curve results. 
There is, in fact, a gradual transition as we proceed along the second row of the 
periodic table from Li to Ne, and also along the third row from Nato A. In each 
case, at the beginning of the next rows (Na, K) a new inner peak of great sharpness 
appears. 

(c) In the particular case of potassium, shown at the bottom of figure 2, we 
have drawn both the total I(p) curve and the components, with correct weight, 
that comprise it. In the total curve, one can clearly distinguish the inner peak 
due to the 4s contribution, followed by a second peak arising from 3s and 3p 
electrons; next comes a still wider peak from the 2s and 2p electrons and there 
is a very low flat peak at large distances arising from the Is electrons. Similar 
discussions could be given for the other atoms, but this sufficiently illustrates 
the way in which the different electrons contribute to I(p). 

A quantity of some importance is the mean momentum p. According to 
what we have said about I(p) and the formula (4), this is the sum of suitably 
weighted contributions p,,, P2,,... Explicit calculation shows that 


Pis=8c/37, Poe=8c/Sm, Pg=128c/1057, Py) =128c/457, yy =9472c/47250. 


Values of p for other orbitals may be obtained most easily by numerical integration. 
Figure 3 shows p plotted against the atomic number Z. Along each row of 
the periodic table the 
variation of p is almost 5-0 
linear, but the slope of 

each successive shell 

is less than the one va 
preceding, and there 

is only asmall change '!& 
in P on converting a 
rare gas to an alkali 
metal atom. All this 
is rather as one would 
have expected because 
p must be related to 
the total electronic 
energy E insome such 
way as p~,/(2E/Z), 
all in atomic units. 
Along each row of the 
periodic table there is 
a steady increase in 
E/Z, but on account of the poor binding of the one valence electron in the alkali 
atom, there is only a small change in £ on starting a new shell. 


Se 
(=) 


2-0 


Mean Momentum 


6 20 24 


8 12 | 
Atomic Number Z% 
Figure 3. 


oe COMP LON SPROE TUES 


It was shown in D.C. §2 that a knowledge of I(p) enables us to calculate the 
shape of the Compton profile. All Compton profiles for a given atom (8, all 
incident wavelengths and all angles of scatter) are included, by suitable change 
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of scale, in asingle formula, where the intensity J at “‘ reduced ” distance ¢ from the 
centre of the displaced band, is given by 


yo=1{ “OS. ee (20) 


Full particulars are given in D.C. and do not need to be reproduced here. In 
that paper J(q) curves were plotted for all atoms up to Ne. Our present calcula- 
tions allow us to extend this work up to K. Five representative examples, each 
giving one-half of the complete profile, are shown in figure 4. From these 
the general trend of 
the curves is apparent. 
These curves are nor- 
malized in each case so 
that J(0) =0. This, while 
making the areas under 
the curves different, 
greatly facilitates com- 
parison of different 
profiles. In general, 
the behaviour of these 
curves agrees closely 
with that found for the 
lighter atoms in D.C., 
that is, s-electrons give 7 

narrow profiles and p- Figure 4. 

electrons wider ones. s 

Thus the sharpest curve in figure 4 is for Mg, whose outer structure is (3s)?, and 
the widest is for A, where the outer structure is (3s)?(3p)®. In the cases of Na and 
Mg there is an appreciable “lag ”’ at the bottom of the curve, similar to that found 
before for Li and Be; this might easily cause trouble when separating the con- 
tributions to the full experimental profile that arise from satellite lines (e.g. the 
Ka, and Ka, doublet) in the incident radiation. 

Finally, we give below, 
in tabular form, the theo- 
retical widths Ag at half- 
maximum value, together 
with the resulting A/ values, 
in X-ray units, for back- 
scattering (9=180°) of Mo 
Kea radiation. This table, 
which may be compared 
exactly with table 1 in 
D.C., p. 195, shows again 
a steady increase, except for 
Na to Mg, along this, the 
third row of the periodic 


table, and the smaller value i 3 = ~ m 
with which each new row Atomic Number Z 
begins. Figure 5. 
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Atom Ne Na Mg 1p A K 
Aq Seo 0-80 0-68 1-48 2:24 1-06 
Al(x.u.) for MoKa A 3 = 
—180° if 4 83 Ho 15 23 11 


Thus in the case of the Compton profile, the changes to be expected between 
one row of the periodic table and the next row are less than between the elements 
of either row alone. his is shown very clearly in figure 5, which collects together 
the Ag values for the atoms both of this paper and the preceding one. 


REFERENCES 


Dirac, 1947, The Principles of Quantum Mechanics, Chapters 3 and 4 (Oxford : University 
Press. Third edition). 

DuMonp, 1933, Rev. Mod. Phys., 5, 1. 

DUNCANSON, 1943, Proc. Camb. Phil. Soc., 39, 180. 

DuNcANSON and CouLson, 1944, Proc. Roy. Soc. Edinb., A, 62, 37. 

DuNCANSON and CouLson, 1945, Proc. Phys. Soc., 57, 190. 

Hucues and Mann, 1938, Phys. Rev., 53, 50. 

Morse, YOuNG, and HaurwiTZ, 1935, Phys. Rev., 48, 948. 

PAULING and PopoLsky, 1929, Phys. Rev., 34, 109. 

SLATER, 1929, Phys. Rev., 34, 1293. 

SLATER, 1930, Phys. Rev., 36, 57. 


Some Experiments with Adjustable Geiger-Muller Counters 


By MOHAMMED CHAUDHRI, Ph.D.(Cantab), M.Sc.(Alig.)* 
ANDY A. Gy FENTON, B.Sc. + 
Physics Department, University, Birmingham 


* Of the Muslim University, Aligarh, India. + On leave from the 
University of Tasmania, Australia. 


MS. received 12 June 1947 


ABSTRACT. Special counters are described in which it is possible to alter the effective 
length as well as the material and diameter of the anode without opening them. ‘This is 
a distinct advantage in that the effects of these factors on the counter characteristics can 
be studied without complication due to changes in the gas filling and the nature of the 
cathode surface brought about by opening the counters or by making separate counters 
with anodes of different diameters. 

It was observed that, while in a given counter the anode material and the nature of its 
surface have little or nothing to do with the characteristics, there is an optimum anode 
diameter for which the best plateaux are obtained, a factor of obvious practical importance 
in counter design. 

The photosensitivity to visible light and longer wavelengths caused by high counting 
rates or by “‘ continuous’ discharges was shown to be due to some surface phenomenon 


on the cathode. 
It has been shown in a simple manner that the discharge is not localized, but spreads 


along the length of the counter. 

It was observed that the plateau length and slope improve at first, but then deteriorate 
when the effective length of the counter is reduced by shielding the anode wire at both ends, 
with glass capillary tubes. With metal capillaries, however, the plateau characteristics 
are best when the full length of the cathode is used. 
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$14. INTRODUCTION 
ucn work has been done on Geiger counters of the self quenching type 
(Korff 1946) but very few systematic investigations have been carried 
out which yield reliable data as to the relative importance of the various 
components of a counter during its operation. Moreover, there are still a number 
of factors influencing the plateau characteristics of counters which have not been 
given due consideration by previous workers. 
It is the purpose of the present paper to study the effects of the following 
factors on the plateau characteristics :— 
(1) The material of the anode wire. 
(2) The diameter of the anode wire. 
(3) The shielding of the anode wire. 
(4) The cathode surface. 
It has also been found during the course of these experiments that the cathode 
surface can be made photosensitive to wavelengths up to about 1 micron, and that 
the degree of the photosensitivity and its duration can be controlled within limits. 


§2, APPARATUS 

The experiments were performed with four counters having nickel cathodes 
and one counter with an aluminium cathode. The nickel cathodes were made 
from sheet nickel 0-1 mm. thick formed into cylinders 2:2cm. in diameter and 
10cm. long, spot-welded along the seam. ‘The aluminium cathode was a tube 
0-1mm. thick, 10cm. in length and 2-0cm. diameter. The cathodes were 
cleaned with nitric acid and were rinsed thoroughly with water and dried before 
assembly into the counters. ~ 

One of the counters with a nickel cathode has an anode made from a number of 
wires of different materials and of different gauges spot-welded together end to end. 
In the discussion which follows we shall refer to this as counter number 1 and we 
give below a detailed description of it and of the other counters. 

The envelope of the counter is a Pyrex glass tube 60 cm. long and 2:5 cm. in 
diameter. The nickel cathode C, figure 1, is closely fitted into the glass envelope. 
- Electrical connection of C with the external circuit is made by means of a nickel ~ 
strip spot-welded to the cathode at one end and to the Pyrex-tungsten seal S at the 
other. 

The anode A consists of wires of gold, silver, platinum, copper, molybdenum, 
tantalum and tungsten, ranging in diameter from 0-065 mm. to 0-25mm. Any 
one of these wires may be set within the cathode by rotating the vacuum-tight 
ground-glass joints E which support the glass tubes F on which the anode is wound. 
The extreme ends of the anode are spot-welded to 1 mm. nickel wires which run 
through the glass tubes F and are held in position by the Pyrex-tungsten seals J 
so that electrical connection to the anode may be made from outside the che 
The anode passes through the 2mm. bore Pyrex capillary tubes D, each 15 cm. 
in length and 5mm. in external diameter. ‘The outer ends of the capillaries are 
tightly fitted into the soft iron cylinders B, each 3 cm. long and 2-3 cm. in diameter 
which slide smoothly within the outer glass envelope. The weight of the bias 
capillaries is supported by two glass tubes G with a bore just greater than the 
external diameter of the capillaries, internally sealed into the envelope. The 
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ground-glass joints K carry thick leads L through which current may be passed to 
heat the anode wire when desired. The capillaries may be set in any desired 
position inside the cathode by moving the soft-iron cylinders from outside with 
the help of a magnet. 

Counter 2 has a nickel cathode similar to that of number 1 but has a fixed anode 
of 0:2 mm. diameter tungsten wire sealed at the ends of the glass envelope, which is 
45 cm. long and 2:5cm. in diameter. The arrangement of the glass capillaries is 
the same as in counter 1. Counter 3 has the same construction as 2 but brass 
capillaries of the same dimensions as in 2 are used instead of glass. In counters 
4’and 5 a single adjustable Pyrex glass capillary (15 cm. long and 5 mm. in diameter 
and 2mm. bore), sliding over the 0:2 mm. tungsten anode wire is used, the motion 
of the capillary being controlled mechanically through a ground-glass joint. 
The cathodes are of the same dimensions but one is of nickel, the other of 
aluminium. 


Figure 1. 


The counters were evacuated with a mercury diffusion pump and were out- 
gassed by warming and by letting the pumps run continuously for several hours. 
The counters were tested for vacuum tightness by allowing them to stand under 
vacuum for several days. 

The gas filling used was an argon-alcohol mixture which was allowed to stand 
ina flask for over 24 hours to mix thoroughly before being admitted to the counters. 

Counters 1, 2 and 3 were filled with the same argon-alcohol mixture containing 
23%, alcohol, to a pressure of 10cm. Hg. Counter 2 had been used previously 
with fillings of different gases but 1 and 3 were new. Counters 4 and 5 were filled 
with an argon-alcohol mixture containing 10% alcohol, to 10cm. total pressure. 
Once filled to the desired pressure the counters were not opened or refilled during 
the experiments described below. 

No wax joint was used in any part of the counters or apparatus and the grease 
used for the taps and ground-glass joints was apiezon M. 

A'scale-of-100 recorder was used to record the counts, and the potentials for the 
counters were obtained from a stabilized H.T. set giving up to 1650 volts. 
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§3, RESULTS 


Effective length and plateau characteristics 


Preliminary experiments with pa tee Re 
‘ : . Glass capillaries (source at 60cm, 
counters 4 and 5 indicated that the 5 Brase) Ube Macemtene een) i 
distance between the ends of the aa 


capillary tubes within the cathode 
determines the effective length of the 
counter. Moreover, it was observed 
that the plateau slope improved when 
the effective anode length was reduced 200; 
by sliding the capillary into the 
cathode. Detailed measurements were 


400 


Counts per minute 


0 2 4 (3 8 10 
Anode length (cm) 
Figure 2. 
8 
Volts Volts 
16 250 ue 
re » 
2 Se 4 
om = 
S12 200 = & -12+200 
ra 22) 2 
8 = S 10) = 
5 ea S 
~ 08; 150 — gr 087150 
g ae. S 
ie = 3 | = A. Pl | 
~ . Plateau slope 
2 -04 A. Plateau slope 100 ae 
s eee = B. Plateau length 
a 


0 2 4 6 8 10 ft) 2 4 6 
Anode length (cm) Anode length (cm) 


Figure 3. Figure 4. 


400 
300 


200 


Counts per minute 


100 


0 
1200 1300 1400 1500 1600 
Anode potential (volts) 


Figure 5. 


Experiments with adjustable Geiger-Miiller counters 187 


made with counters 2 and 3, in which Volts 
both capillaries are adjustable. Figure ie 
2 shows the relation connecting the 
length of the anode wire exposed to 
the cathode between the ends of the 
capillaries and the counting rate at a 
working potential of 100 volts above 
the threshold (i.e. about the middle 
of the plateau) when a gamma-ray 
source is kept at a fixed distance from 
each counter. The relation between 
the effective anode length and the slope 
of the plateau for counter 2 is shown in Figure 6. 

figure 3, curve A, while curve B shows 

the relation between the effective anode length and the length of the plateau. 
The curves indicate that as the anode length is reduced from the maximum of 10cm. 
the plateau slope and length improve at first but later deteriorate. Similar curves 
in figure 4 show that for counter 3 the plateau is best at 10cm. anode length and 
begins to deteriorate as the metal capillaries are slipped into the cathode. Figure 5 
shows the characteristic curves of counter 2 for different lengths of the anode. 
The variation of the threshold voltage with the effective anode length is shown in 
curves A and B, figure 6, for counters 2 and 3 respectively. As the effective length 
is reduced below 2 cm., the threshold voltage increases rapidly until at about 1 cm. 
there are no counts even at 1600 volts. 


A. Glass capillary shielding 
B. Metal =» ” 


1400 + 


Threshold potential, 


i300 + 


1200 n ma: er re 4 


Anode length (cm) 


§4. INFLUENCE OF THE GAUGE AND MATERIAL OF THE 
ANODE ON PLATEAU CHARACTERISTICS 


The plateau characteristics of counter 1 were obtained using anodes of different 
materials and of different diameters. As the anode could be changed by rotating 
the ground-glass joints, other conditions could be kept constant during the 
measurements. A 4cm. anode length at the centre was used throughout these 
tests, as earlier work had shown that good plateaux are obtained with this length. 


Figure 7, curve A, shows the NU nateau sions 
relation between the plateau nie ee PSs pe os 
slope and the wire diameter, 16 {1400 : 
while curve B shows how the 
length of the plateau varies with 
the wire diameter. It was found 
that anode wires of a given 
diameter give the same plateau 
characteristics although made 
from different materials. ibe 

Figure 7, curve C, shows 
that the threshold voltage de- 1000 , _ eee | 
creases when the diameter cf ‘ ae es ieee oa 
the anode is reduced. Fionre 7 
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§5. PHOTOSENSITIVITY OF THE COUNTERS 


Early experiments with counter 4 showed that after passing a “‘ continuous” 
discharge * through the counter (by setting the voltage a few hundred volts above 
the plateau region) the background counting rate immediately afterwards was of the 
order of 10000 counts per minute, but after 24 hours it had fallen to 200 per 
minute, the normal background being about 50 counts per minute and the plateau 
characteristics were slightly better than before the discharge. When counter 5, 
which has an aluminium cathode, was new, it had plateau characteristics as good as 
counter 4, but after a continuous discharge it required a considerably longer time 
to recover and was thereafter unreliable for counting rates of over a few hundred 
counts per minute. In subsequent measurements, which were carried out in 
greater detail with other counters, it was found that these effects were primarily due 
to changes in the behaviour of the cathode surfaces, produced by the discharge. 

Counter 2, with which a few million counts had already been recorded over a 
period of some weeks without making it photosensitive, was run at a high counting 
rate at an operating voltage about 100 volts above the threshold. In a particular 
case, an anode length of 3 cm. was used and gamma rays were counted a a rate of 
about 6000 counts per minute. Background counts taken with the room lights on, 
following doses of varying duration, showed that a detectable photosensitivity was 
built up after counting at this rate for 10 minutes. After a few minutes, the 
background counting rate fell to within probable error of the normal value. ‘The 
degree of photosensitivity and the period for which it lasts depend upon the dose 
given. 

The photosensitivity caused by passing a continuous discharge was studied in 
some detail with counter 1 in order to find out what part of the counter becomes 
photosensitive. ‘The procedure adopted was as follows: 

The capillaries in counter 1 were set so that a 3cm. length of the anode, 
starting 1 cm. from one end of the cathode, was used during the discharge. Since 
the cathode is 10cm. in length, a fresh portion of it is used when the capillaries 
are adjusted for a similar 3cm. length at the other end of the counter. After 
passing a discharge for 5 minutes, the background counting rate with the room 
lights on was about 20 times the normal value, but the rate on moving the capillaries 
to the other end of the counter was normal. Thus it appears that this portion of the 
counter had not been affected at all by the discharge, and that the gas filling retained 
its quenching characteristics. Next, the portion of the anode which had been 
used during the discharge, was wound along until it became the anode for the other 
end of the counter. Again, a normal counting rate was found, while the photo- 
sensitivity persisted at the end in which the discharge had been passed. This 
means that the portion of the cathode which was used during the discharge had 
become photosensitive. The sensitivity is greater and persists for a longer time 
when the duration of the discharge is increased. 

After a discharge had been passed in this counter for 1 hour, the photosensitivity 
decayed as shown in figure 8, curve A. The sensitivity decreased rapidly at first 
and then more slowly. In one counter appreciable photosensitivity remained after 
it had been standing unused for a period of several months. The data of curve A, 
figure 8, are shown on a double logarithmic plot in curve B. 


* When the potential is set above the Geiger region, the discharge nee passes is not really 
continuous, but consists of discrete pulses in rapid succession. 
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Two days after a continuous dis- re 
‘charge had been passed for 1 hour in 
‘counter 1, the spectral sensitivity was 
tested qualitatively with a set of Wratten 
filters and a tungsten-filament lamp, as 
after this time the decay of the photo- 
sensitivity was too slow to alter appreci- 
ably during the counting intervals. It 
was found that, under these conditions, 
the maximum sensitivity in the visible 
region was at the short wavelength end 
of the spectrum and there was no de- Time (min) ° 
tectable sensitivity above about 5500 a. ‘0 a a 5 é 
Immediately after a continuous discharge Log time (minutes) 
however, the counter was always sensitive eases 
to longer wavelengths than this, but measurement was difficult because 
of the rapid decay rate. 

§6. EFFECT OF A CONTINUOUS DISCHARGE ON THE 
PLATEAU CHARACTERISTICS 

The plateau characteristics were obtained for two 3cm. effective lengths of 
counter 2, each set at a distance of 1 cm, from the ends of the cathode. A discharge 
was then passed for 1 hour in one of these sections, and immediately afterwards it 
was observed that the slope and length of the plateau as well as the background in 
the other section with room lights on, were, within probable error, the same as 
originally. On the other hand, the background for that section of the counter in 
which the discharge had been passed was too high for the recording apparatus to 
follow. In darkness, however, the background for this end of the counter was the 
same as before the discharge and equal to that of the other end in light. ‘Though 
the background of this part of the counter was too high in light for several days for 
the recorder to follow, its plateau characteristics taken with room light off were as 
good as those obtained in light before the discharge. 
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§7. EFFECT OF A CONTINUOUS DISCHARGE ON THE 
BACKGROUND IN DARKNESS 


During the experiments with 35718 x10? 
counter 1 it was observed that 
the background counting rate 
with the tube in darkness was 
higher than normal immediately 
after the discharge had been 
passed in it for an hour, but 
returned to the normal value 
within 48 hours. This phenom- 
enon was studied in detail with 
counter 4, in which the effect was 
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darkness during the 20 minutes interval immediately after a continuous discharge 
of 30 minutes. It was found that this effect disappeared altogether after a 
number of discharges of short duration had been passed in the counter, although 
it could still be made highly photosensitive. This was found also to be the case 
with other counters. 


§8. DISCUSSION OF RESULTS 


During the experiments the pulses from counter 2 were examined with a. 
cathode-ray oscillograph and it was observed that for a given operating potential 
the pulse amplitude increased almost linearly with the effective anode length, 
which is consistent with the assumption that the discharge is not localized but 
spreads along the length of the wire. This demonstrates in a simple way the 
conclusion arrived at previously by other workers (Wilkening and Kanne 1942, 
Stever 1942) and studied recently by Hill and Dunworth (1946), Craggs and Jatte 
(1947), and Wantuch (1947). 

The linear relation (figure 2) between the separation of the ends of the 
capillaries and the counting rate shows that the efficiency of the counter is un- 
altered by reducing its length and also that for practical purposes we are justified in 
referring to the separation of the capillaries as the effective anode length of the 
counter. At very short lengths, however, the distortion of the fields near the ends 
of the capillaries would be expected to have a noticeable effect and this appears to be 
the case, since the counter with glass capillaries was not sensitive to gamma rays 
at an anode length less than 0-8 cm. with a working potential of 1600 volts; while 
the metal capillary counter at the same operating voltage would not respond below 
1:25 cm. anode length. ‘ 

The slope of the plateau in self-quenching counters has been attributed by 
Korff (1944) to the occurrence of multiple pulses which, with a recording circuit 
of high resolution, give rise to spurious counts. On this view the more ions there 
are per pulse the greater is the probability of a multiple pulse occurring. ‘There- 
fore, when the pulse size is reduced by reducing the anode length of a counter, one 
would expect to obtain a longer and flatter plateau. This would account for the 
initial improvement in the plateau characteristics for the counter with glass 
capillaries (figure 3) and for the increase in the upper limit of the plateau as the anode 
length is reduced (figure 5). Figure 5 shows that as the anode length is reduced 
below the 4 cm., the threshold moves to higher voltages more rapidly than the 
upper limit of the plateau, resulting in a reduction of the plateau length. This 
increase of threshold potential with decreasing anode length (figure 6) may be- 
attributed at least partly to the sensitivity of the recording apparatus which does 
not respond to pulses of less than a certain minimum amplitude. Observation 
with a cathode-ray oscillograph during the experiments showed that this minimum 
was reached when the anode length was reduced to about 4cm. with a working 
potential of 1250 volts. At shorter lengths than 4cm. the plateau did not begin 
until the pulse size was increased by using higher operating potentials. The 
increase in slope of the plateau for short anode lengths suggests that spurious: 
counts due to factors other than the number of positive ions per pulse were 
recorded. ‘These may be due to point discharges from irregularities on the 
surfaces of the electrodes or to distortion of the field at the ends of the capillaries, 
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both of which effects would be more noticeable at the higher potentials used. 
Furthermore, at higher potentials, the average energy of the positive ions being 
greater, there is a greater chance for some of them to reach the surface of the cathode 
and liberate secondary electrons, thus giving rise to spurious counts. In the 
counter with the brass capillaries, the best plateau was obtained at 10cm. anode 
length, suggesting that the expected improvement in the plateau during the initial 
stages on reducing the anode length was masked by effects caused by the presence 
of the metal capillaries. 

In order to explain the high background in darkness, observed after a 
continuous discharge (figure 9), one must look for processes occurring on the 
cathode surface in which electrons or negative ions are emitted. One possibility 
is that during the discharge active centres are formed on the cathode surface, and 
these undergo chemical change with the emission of electrons or ions. Also, if 
during this chemical change, positive ions are formed on an insulating layer, e.g. 
oxide, on the cathode surface, these may give rise to field emission electrons. On 
the other hand, the photosensitivity which is observed when the cathodes are in an 
activated condition after a high counting rate or a continuous discharge suggests 
that centres of low work-function are formed. With a sufficiently low work- 
function, and under the influence of the electric field at the cathode, an observable 
number of thermionic electrons would be able to escape at room temperature, thus 
accounting for the high background in darkness. For example, calculation shows 
that, if we regard the activated cathode as similar to an oxide coated thermionic 
emitter for which A in the emission equation 7= AT%e—*/F* is 0-03 amp. per sq. cm. 
and ¢ is 1-2 volts, the current at 300° K. is about 10 electrons per sq. cm. per second. 
This view is supported by the observation that the counter was sensitive to wave- 
lengths up to about 10000 a. when the high dark background was also present. 
Furthermore, the counter was sensitive to the long wavelengths only for a short 
period immediately after a discharge, but sensitivity at the low wavelength end 
of the visible spectrum remained for a long time, suggesting that the work function 
of the centres was low at first but increased with time. It is also worth noting that 
the counting rate in darkness increased when the counter was slightly warmed with 
hot air from a hair dryer, an effect quite distinct from the photosensitivity at long 
wavelengths. 

High background rates n darkness after high counting rates have been reported 
by several workers. Lauterjung and Neuert (1944) found the effect most marked 
for cathodes of low work-function such as Mg (photoelectric work function 2:4 
volts), whereas in the present work the effect was observed with nickel cathodes 
with a work function of about 5 volts. Christoph (1935) reported that in photon 
counters with Cd or K cathodes the background immediately after a glow dis- 
charge was much higher than normal. He also observed an increased photo- 
sensitivity after the glow discharge. Lauterjung and Neuert obtained linear 
decay curves using double logarithmic plots, whereas in our case the curve is 
linear at first and later changes slope (figure 9, B). 

The upper portion of the curve of figure 8, B, representing the first 30 minutes 
decay in light after a continuous discharge for 1 hour is linear and has a slope of 
unity which suggests that the natural decay curve follows a rectangular hyperbola 
law during that period but afterwards the decay becomes more complicated. 
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After 4 weeks the background in light was about twice the dark background and 
the decay curve became linear again, indicating that if it were to continue in this 
way, a measurable sensitivity should remain after 10® minutes (i.e. almost 2 years). 
Actually, one of the counters was found to be photosensitive 4 months after a 
continuous discharge of 3 hours duration had been passed in it. 

Spatz (1943) reported that counters with silver cathodes became photo- 
sensitive to visible light after long use, and showed how the plateau was destroyed 
by running the counters for some time at a high counting rate, but gradually 
recovered again over a period of several days. In view of the present work, we 
may suppose that the results obtained by Spatz were probably due to the effects 
observed by us and discussed in this paper. 

Figure 7 shows that for a particular counter there is an optimum anode 
diameter which results in a long, flat plateau, a point which is of great practical 
importance. 

It was invariably found during the work that the threshold occurred at a 
slightly lower potential and the plateau was longer when the counters were new or 
operated after a rest period of afew hours. Some minutes of operation, however, 
brought the threshold to its higher steady value and the plateau characteristics 
thereafter were reproducible so long as care was taken not to allow the counter to go 
into a continuous discharge. Observation with an oscillograph showed that 
after a rest period the pulses were larger, and this may explain the lower 
thresholds. 

It was also observed that the degree of oxidation of the anode wire does not 
influence to any appreciable extent the characteristics of a counter, because just as 
‘good plateaux were obtained using.a heavily oxidized copper anode as with a gold 
wire of the same diameter. Moreover, no difference was found when the anode 
wire was slackened so that it was about 2 to 3 mm. off the axis at the centre of the 
‘counter. 

It is generally believed that after prolonged use the plateau characteristics of a 
-counter deteriorate as a result of the dissociation of the quenching agent. Our 
results suggest, however, that this may be more closely connected with the changes 
on the cathode surface than with the breakdown of the filling. 

It is hoped to continue the investigation of some of the aspects discussed 
above. 
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ABSTRACT. Ascheme of measurements is discussed from which absolute determinations 
of sound field pressure can be made. The relation between the scheme and the reciprocity 
method is discussed. It is shown that the measurements also yield the axial pressure 
and “‘ projection efficiency’ of a transceiver, the “‘ projection efficiency” being defined 
as the ratio of the actual axial intensity to the axial intensity as it would be if the transceiver 
converted the whole energy absorbed into sound by vibrating 23 a simple piston. It 1s- 
suggested, without proof, that the upper limit of the projection efficiency is unity, and the 
relation between this quantity and the efficiency defined as the ratio of radiated acoustic 
power to consumed electrical power is shown to depend on a directivity coefficient and on 
the effective area of the radiator face. It is pointed out that neither of these quantities 
is readily determinable and that therefore the conversion efficiency of the transceiver, 
and hence the acoustic power radiated, cannot be obtained. 


Silk TWANAP ROD NUKE APIO] 

HE determination of the axial pressure, efficiency and directional 
properties of crystal transceivers used in submarine acoustic work 
presents many problems which, in the urgency of immediate need, 

have tended to remain without complete solution. ‘The present paper is a report 
of a theoretical investigation into the possibilities of making exact measurements. 
on the quantities specified for crystal transceivers.* 

An obvious initial requirement is a device for the absolute measurement of 
sound field pressure. Such devices are provided by the Rayleigh disc and the 
radiometer disc (which latter measures energy density in the wave) and present 
considerable technical problems of their own. On the other hand, it is known 
that by the use of the reciprocity theorem it is possible to measure sound pressures. 
without any ‘‘absolute” device. The method described in this paper is closely 
related to that based on the reciprocity theorem (Foldy and Primakoff 1945) 
but in its simple form it was developed independently of this result. 

* Security reasons still prevent the publication of actual numerical data, but the principles 


involved are of general interest in acoustical work and it has been thought worth while to prepare 
the theoretical work for publication. 
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The reciprocity principle presents the concept that two transceivers placed 
in a sound-transmitting medium are equivalent to an electric quadripole, and 
that the transfer impedance is therefore independent of the input direction. 
In other words, the ratio of open circuit volts on the receiver to driving current 
in the projector is unique. A concomitant property of the system is that the ratio 
of the receiving sensitivity of a transceiver, measured in open circuit volts per unit 
pressure, to its projection sensitivity, measured in pressure units at a given 
distance per unit driving current, is also unique. ‘This ratio has been called the 
““reciprocity parameter” and is denoted by J. 

The basic principle of the method described in this paper circumvents the 
reciprocity principle, in effect, by making the analysis rest on voltage measure- 
ments together with self impedance measurements on the transceivers. Never- 
theless, the full analysis of the method, as we shall see, cannot be carried out 
without invoking the reciprocity principle. 


$2, THE BASIC EXPE RENEE Ne 


Let us assume that we have a perfect transceiver, ‘[', for which axial pressure 
is required. We shall require an auxiliary hydrophone, H, and a projector, P, 
for producing a sound field. The hydrophone may be of any form whatever, 
provided that it is linear, and it may be supposed to be followed by an amplifier 
and an output meter. If the output meter records a voltage uw when the hydro- 
phone is placed in a field of which the undisturbed pressure is p, then we may 
write 

p=hu, sesmemicede 

where A is the calibration constant-of the hydrophone. 

The transceiver, I’, may be represented as 
in figure 1, in which J is the motional impedance 
of the device, e is the voltage produced in it by 
the sound pressure when it is acting as a receiver, 
and C is a capacitative load. ‘The first part of di 
the basic experiment is to establish a sound C 
field using the projector P, and to determine 


the voltage developed across C when T is placed e 
in this field. If p,, be the free field pressure, 
we have * 
2A 
e= + fos eee (oes) Figure 1. 


where A is the effective area of the face of T and k is a constant involving the 
piezo-electric constant of the crystal and the dimensions of the assembly. 
The output voltage, v, developed across C is 


tp 3 Oy Del al 
Lem PE has Anny ° 


where K is the conductance of C. 
Since the measured value of v gives no indication of phase, the modulus of 
the impedance factor 1+ KI must be used. 


* See (3.1) and (3.2), infra. 


Ree Ces 
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If, now, T is replaced by H, the output meter of H will record a voltage u,, 

given by 
ae i A eee (2.4) 

Let 'T’ now be used as a projector, a driving voltage V being applied to its 
terminals. ‘The power, W, radiated acoustically will be 

v2 
ne 
where S is the radiation resistance. If the dimensions of the radiating face of T 
are not small compared with the wavelength, A, the axial pressure, p, at a distance 
d( >A) is given by * 


W Se gi Oo) (2.5) 


k k VS 
p=5,VSW=x, bay te bh orcs (2.6) 
If this pressure be recorded on H, an output voltage wu will be produced, given by 
Dil Dod heli Se ES a eae (2.7) 
From (2.3) and (2.4) we have 
U 2Ah 
Ay = aaa SU B'S (2.8) 
From (2.6) and (2.7) we have 
u kS 
V = VAD Sw CC (2 9) 
From (2.8) and (2.9) we find 
2AS 
MyM, = d|T\}1+ KI)’ sheKeneieie (2.10) 
my 2Add (| 
ES h Timers 9 ee ee (PAA) 


For a transceiver with a radiating face of dimensions not small compared with i, 
we have (Vigoureux 1947) 
eS Ae DC ee we oP mee Ee Rereatate (2:12) 
where p is the density of the medium in which T is placed and c is the velocity 
of sound in the medium. 
If we can ensure that K/> 1, (2.10) and (2.11) become 


2AS 2k? iS 
eee eo ees at 
mila= TR Tp IK Adpg * [TTP canes 
2rd 
Ht |? = oer Kl Wee ee (2.14) 


The former relation depends only on the properties of ‘T’ but not of H, and the 
latter on the properties of H but not of T. In particular, (2.14) allows us to 
determine h, the hydrophone constant, directly in terms of voltage ratios, impe- 
dance and frequency. From (2.6) and (2.13) we have, further, 


bk 8 _& fmmdielR 
V~ Ad *|T| ~ Ad oe 
=f mmapelk ee (2.15) 
Ties: 


where f=c/A is the frequency of the sound. Thus the axial pressure is determined 
in terms of directly measurable quantities. 
* See (3.1) and (3.2), infra. 
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Finally, if we make the measurements at the resonant frequency of T, J 
becomes identical with S for our assumed perfect transceiver, so that (2.13) 
becomes 


2A 
MyM, = | K|adS’ CLOSED 
which suggests that by working at the resonant frequency of T we might determine 
its radiation resistance. 

We have shown then, in principle, that the system of measurements outlined 
above (to describe which the adjective “trimetric” is proposed) permits the 
absolute measurement of free field pressure and of certain constants of trans- 
ceivers without invoking the reciprocity principle explicitly, provided certain 
conditions are fulfilled. ‘The most important assumption made in the analysis 
has been that the transceiver has been vibrating as a perfect piston and has 
therefore possessed directional characteristics in accord with the geometry of 
the system. The implications of this will be discussed more fully in the 
following paragraph. 


$3. EP FE © TOF DIRT CAL@ ING Estey, 
The sound power, W, radiated by a perfect piston of effective area A, in an 
infinite baffle is 
W= Ae -- eee (3.1) 
where £2 is the mean square velocity of the radiating face (Rayleigh 1929). The 
axial pressure, po, at a distance d from this source is 
cA; 
tee = on re. (3.2) 
where A is the geometrical area of the radiating face (Vigoureux 1947). 
It is implicit in (3.1) and (3.2) that d>.1/A>A. Thus 
_ Ad2p2Ay 
= ~ pcA® ty  _ ng Fe See eee 
Consider now a radiator which is not a perfect piston and let its axial pressure 
at a distance d be p, when it is radiating a total power W. ‘Then the pressure 
at some point on the surface of a sphere of radius a centred on the piston can be 


written as ; 

ple. ) =p:f(0, $), Si ereeiehs (3.4) 
where 6 and ¢ are the angular coordinates of the point considered and f is some 
unknown function. ‘Then the rate of flow of energy radially through an element 
ds of the surface at this point will be 


{p(8, 4)}"ds 
pc 


W 


2 
= 7 {f(0, $)}*ds 


so that we have 


w= PE] (900, $)}%ds, 


where 


A=| (f(0F8)}*ds. 


Calibration of hydrophones and crystal transceivers 197 


For the perfect oscillator the relation between radiated power and axial — 
pressure is given by (3.3), so that from (3.3) and (3.5) we have 


; eA? a, 
Pig Po ane <a) ie ee (3.6) 
or 
DiS Oyler (3.7) 
where 


B=Ad/V A2A/A,. 


The factor B may be called the “directivity coefficient” for the radiator. 
(It differs from the “ directivity index” which has been used by other writers 
in that we have taken as a reference standard a perfect radiator of the same form 
as our imperfect one, instead of a spherical radiator.) 

It is worth while digressing for a moment to consider the form of A for the 
particular case of a circular radiator of radius a. This system possesses axial 
symmetry and 


2J,(bsin@ 
je ey 3.8) 
where J, is the Bessel function of order one and b=27a/A (Morse 1936). Then 
m{2 
= Snd?| [J,(6 sin @)/b sin 6]? sin 6d8, 
0 
d** ae) 
a, E eee) 8 ee (3.9) 


Since for a perfect piston we must have 6 = 1, it follows that 


2 2 J2)\2 
IN ea Re 


whence, since 
A=na’, 
J,(26) 1 
4j=4)1- eae (3.10) 


This is the result given by Rayleigh (1929) for a circular piston. In most of the 
transceivers used the value of 6 was such that {J,(2b)}/b was generally very small 
compared with unity, and we can take 4,= A. 

This, then, has given some idea of the part played by the directivity coefficient 
in the behaviour of the transceiver as a projector. We have now to decide what 
part 8 plays in the behaviour of a receiver. It is arguable physically that the 
generated voltage e (figure 1) for a “piston” of directivity coefficient B must 
be B times that for a perfect piston of the same materials and dimensions, since 
8 implicitly describes the mode of vibration of the radiating face however excited. 
It is, however, interesting to give a formal demonstration of this by invoking 
the reciprocity principle. . 

Referring to the crystal transceiver represented diagrammatically in figure 1, 
let the axial pressure w, produced by a driving current i be given by 


Sh eS ae eee (3.11) 
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and let the open circuit voltage v produced by a free field pressure wy be given by 
O=VOq. | | eer (3.12) 


Then the reciprocity theorem states that 
iad 4,2) ees (3.13) 


In figure 1 the actual current entering J is 


, dee 
ol Rhg 
and the power radiated acoustically is 
= Se (oa en AER AY 
W=S\1,|4= t+KN 2" : 
= pcA,&. 
The axial pressure for an ideal radiator is 
pcd ; 
Po= hd g, 
Ey el yee 
=" (eK Nowea a ee cde (3.14) 
The actual axial pressure p, is 
7, beppr oe ee SA 
Pi=BPo=. dy (iERTN oe eas a (3.15) 
sha eee 
Thus PS ae \TERTIN pe? (3.16) 
2dr 1 Nee 
d EIS oi) De ee eee 
an | v Ee 28 Tea rey We ere (334) 
Thus, in a free field p, 
V=vVp, - 
1 ie 
= 180 ier Ni aes le WE Saas (3.18) 
If e be the generated voltage 
1 
}1+AD)? 
SA? 
so that a ies ; ; 
e=2Bp roy Vannes pe (3.19) 


For an ideal piston the generated voltage is given by 


SA? : 
eo 2p, =¢/ Bs Ceceee (3.20) 


so that the generated voltage for an imperfect receiver is equal to its directivity 


coefficient multiplied by the generated voltage for an ideal receiver of the same 
materials and dimensions, att tt 
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§4. THE TRIMETRIC SCHEME: HYDROPHONE CALIBRATION 


We can now develop more fully the analysis given in § 2 for an ideal piston. 
We have, in fact, to rewrite the equations of that section in the light of the 
discussion given in § 3. 

We have instead of (2.3) 


pao aria 
a ear eee (4.1) 
and instead of (2.6) 
kVS A 
p Ad || ry ey ee NRL 6 sulsfeer a (4.2) 
Combining these with (2.4) and (2.7) we have 
5 RE le 
ee Wace ig PSCU AES: ote Ae: (4.3) 
i= sats BRS a 4.4 
Seen Pee ner (4) 
whence 
. PR ia bei) 
Hi WEE See uy | ee (4.5) 
m,  2A,dd e 
aes | wee a ee (4.6) 
Instead of (2.12) we must use the more precise form 
(sete ay Nore hae ak amar go (4.7) 
so that 
_ my 1+ ~). of 
mt Dn hee oe (4.8) 
which we may write as 
Stee ey ae ake 
a Fl ee es cv 9 (4.9) 


Thus the calibration of the hydrophone does not depend on the directivity 
coefficient of T nor on the relation between its effective and geometrical areas, 
nor does it depend on measurements being made at or near the resonant frequency 
of T. Ifthe load on T is a capacitance C, as in figure 1, having a power factor 5, 
then | Kj can be written as wC to an accuracy of 0-58, and since 6 is never likely 
to exceed 5°%, we can in general ignore this correction. 

The correction 1/KI will be small even if we work at the resonant frequency 
of T, for then J will approximate to the radiation resistance. This suggests that 
if the correction be not very large, use can be made of the ideal equation (2.16), 


which gives, 


va | 


Admym»,? 


|K|S= 


so that, bearing in mind that we are dealing with r.m.s, quantities throughout, 
we have Cop? a aie 
»_ my7Cp ( M1? M2r2 ) 
Ms aa 1+ maa gam): ee ene (4.10) 
We shall see, however, in general that [1/K/J]* is in fact negligible compared 
with unity. 


14-2 
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See Anelis, TRIMETRIC SCHEME: TRANSCEIVER CALIBRATION 
To determine the axial pressure of T we combine (4.2), (4.5) and (4.7), and 


obtain 
myMopf IK | 1 
pa vf mms wae (5.1) 


This again does not depend on 8, since the product mm, is experimentally 
determined, nor on the ratio of effective and geometric areas. 
Disregarding the corrections for 1/KI and 8, in the light of what has already 


been discussed, we have 
C 
Ds vts/ ae oe a ie (5.2) 


Both this result and (4.10) are independent of any resonance conditions, 
and therefore we can measure absolutely the axial pressure at any frequency 
and so plot the true acoustic projector response curve. 

We saw in §2 that the possibility of determining radiation resistance arises 
from the present series of measurements (equation (2.16)). Actually, it appears 
from the more complete expression given in (4.5) that the particular property 
of the transceiver which does exist and which is directly determinable from this. 
trimetric series of measurements is the quantity B?.A*,S/A,| Z|”. 

In order to find what significance can be attached to this quantity, suppose 
that measurements are made at the resonant frequency of '[—that is, at the 
frequency at which p/V, defined by (5.2), is a maximum. At this frequency 
|| becomes a pure resistance, the value of which may be denoted by S;. This. 
will in general be greater than the radiation resistance S, because it will include 
internal losses arising from mechanical and other causes. Any electrical method. 
which determines total power loss in the transceiver can be used to determine S, 
for, clearly, at a known r.m.s. driving voltage V the total power consumed will. 
be V?/S,. At resonance, therefore, we can determine a quantity R given by 

Ri Ay s ne = eae 
AGS .mmAdr{Cn 97 5) eae 

If the whole power consumed were radiated acoustically by an ideal piston 
of effective area A, identical with its geometrical area, the axial pressure py at: 
distance d would be given by 

esas 
Po rN2q2° so Tee. ioe oh Wooo 
The actual axial pressure, p, is, however, given by (5.2) 


TM, MspC 


Dyer 
p= pepe, 


“Thus 


pM ER ca eee een (5.6) 


5 Thus with a knowledge of R and S,, we can determine the ratio of the actual 
axial intensity to that which would be given by an ideal piston radiating the whole: 
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of the power as sound. ‘This ratio may be called the “ projection efficiency”’, 
and we can denote it by y).* Thus 


S,  mdCcSymym, 
ALS Aa ny, a 


It is obvious that the efficiency with which 'T converts electrical power into 
sound power is given by the ratio S/S,. This may be called the “ conversion 
efficiency ’’, 7, and we have 


aS ea 
The = iS ~ IP Baa 


Thus the actual conversion efficiency, or the actual sound energy radiated, 
cannot be determined from this scheme of measurements unless 6 and A,/A are 
known. ‘There is no accurate method of determining f: for axially symmetrical 
systems an estimate can be made from directional curves, but the method is 
unsatisfactory. Again, theoretical estimates of A,/A rest on the assumption 
that f is unity, so that it appears that 7, is a somewhat inaccessible quantity. 

This, however, need not concern us unduly, for the projection efficiency 
is an equally useful practical concept for asdic applications. While it is obvious 
that 7. 1, yet it is not immediately clear that y)>+1. ‘This has led to the enuncia- 
tion of the following proposition :— 

‘A plane radiator consuming, and radiating the whole of a given amount of 
power produces a maximum axial intensity at a given distance when it is vibrating 
as an ideal piston.” 

No formal proof has been evolved, but there is a strong supposition that the 
proposition is true. If it is true, then of course the upper limit of 7) would be 


unity. 
§6. CONCLUSION 


It has been shown that the scheme of measurements discussed in this paper 


can be made to yield the following results :— 
(i) The calibration constant for a hydrophone in “‘ absolute” units. 
(ii) The axial pressure, at any frequency, of a transceiver, in “‘absolute”’ 


units. 
These results depend only on directly measurable quantities such as voltage, 


capacitance, frequency and range. If measurements are made at the resonant 


frequency of a transceiver, then there is determined also: 


(iii) The projection efficiency. 
It is hoped that it may be possible later to publish practical details of the 
technique developed from this paper. The following results indicate the order 


* The choice of an ideal piston without edge-correction as a norm is arbitrary, and it is arguable 
that it would be better to take a simple piston with appropriate edge-correction as a reference 
If this is done, we find that p?|po2=(S1/R)(Ai/A), which is in general smailer than the 
A third possibility is to take a small spherical radiator as a 
po? would depend on A/d, and no convenient physical inter- 
The present suggestion for the definition of np is tentative, 


standard. 
quantity 7p as‘ defined in (Sy 

standard, but then the ratio of p” 
pretation can be placed on the result. 
and may require revision as experience is acquired. 
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of magnitude of the quantities involved in the calibration of a small tourmaline 
hydrophone :— 

Mm =9:33-.10" f=3-06 x10he/se 

fis— 1-087 0, p= LOA orice, 

C=4-01 x 10-® farad, Sb = 3-3 lec10Zohms, 

d=6:10< 107) em. 
We note first that 1/KI= —j/27x3-06x4-01 x13-31x101=—004 so 

that |1+1/KIJ| differs from unity by less than 1%. To this accuracy then we 
have, from (4.10) 


poe aCpf? 
Sine ed, a 
io ae ; 
or h=5'70f cae dynes/cm? per volt, 


where C is in pyr., fis in kc/sec., dis in cm. and zp is embodied in the constant. 


Sa 9-33 x 4:10 
== SCX 306, (Feceto x 10°), 
= 4-20 x 10° dynes/cm? per volt. 

This gives the hydrophone constant as already defined. Alternatively,. 
the sensitivity can be expressed as 1/h, i.e. as 2°38 ppv. per dyne/cm? 

It can be further shown that the quantity R (equation (5.3)) is, for the present 
cases, 6:1510* ohms, and S$?4/A,=S,?/R=18000 ohms. Approximately, 
A= A,, so that S8?=18000 ohms, which means that the actual radiation resistance 
is of this order, although nothing“is known of 8 except what might be inferred 
roughly from directional curves. Nevertheless, the example does serve to bring 
out the significance of the three ‘‘ resistances’”’ associated with a transceiver. The 
true radiation resistance, S, is that arising from the damping due to the medium 
(in our case, water), the dissipation resistance S, is the equivalent shunt resistance 
responsible for the absorption of electrical power, while the third quantity, which 
we have denoted by R, can be regarded as the radiation resistance of the equivalent 
perfect piston. ‘The magnitudes will be in the order R>S,>S. 
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ABSTRACT. ‘The method described is based on the Sommerfeld formula for diffraction 
at an edge, combined with the effect of reflection at the earth. A table is given which 
reduces the amount of numerical work involved. 


ST LNGRO DUC LLON 


N the winter of 1944-45 it was found necessary to use narrow-beam radar 
] equipments at wavelengths of the order of 10 centimetres for locating 

targets at small angles of elevation. - As the edge of the beam struck the ground 
near the equipment, there was a considerable amount of “‘ ground-clutter”’ which 
greatly increased the background of noise as seen on the display, thus swamping 
the echo from the target. With a view to reducing this clutter, calculations and 
experiments were carried out on the use of opaque vertical screens placed in front 
of the equipment, providing an artificial horizon slightly above the true horizon. 
The experiments were carried out by Army personnel in the campaign in Holland, 
and by Hey, Parsons and Jackson of Army Operational Research Group in this. 
country. 

By a well-known theorem, the polar diagram of an aerial system is the same for 
transmitting as for receiving; it is more convenient to consider the equipment as a 
receiver. ‘The treatment adopted is to sum with respect to phase and amplitude 
all the rays which reach the aerial, deducing the relative phases from the path 
differences, and the relative amplitudes from the Sommerfeld expression for the 
amplitude diffracted through a given angle. 


SA, WSN T ONG 


The arrangement of screen SS’ and receiver or transmitter ‘T considered is 
shown in figure 1. The calculations are simplified by considering T’ to be a 
receiver and by considering plane waves coming from the right at an angle 6. In 
all there are ten rays to consider, which can be put into four groups. Group 1 
consists of rays which have suffered no reflection; Group 2 consists of 
rays which have been reflected before reaching the screen; Group 3 consists of 
rays which have been reflected after reaching the screen ; and Group 4 consists of 
rays which have been reflected both before and after reaching the screen. 
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The component rays in each group are: 


Group 1. Direct ray TN, if @>A4; 
Diffracted rays TSN, and TS‘Ns. 


Group 2. Reflected ray TQ,N, if Des 
Diffracted rays TSQ,N, and TS’Q,Ns. 


Group 3. Reflected ray TQ,N, if 6>C; 
Diffracted rays TQ,SN, and TQ;S’Ns3. 
Group 4. Diffracted rays TQ,SQ;N, and TQ;S’Q.N;. 
All rays in a given group arrive at the equipment in phase except for the 


phase changes on diffraction.* If the resulting complex amplitude is A, for 
Group 1, A, for Group 2, etc., and the reflection coefficient is R, we have: 


—4rhy sin @ 
Serie 


The reflection coefficient R is in general a complex quantity which varies with 
the angle of incidence. For horizontal polarization R does not differ greatly 


Resultant field = A, + A,R?+ (A, + A3)Re 


Figure 1. 


from —1, and for vertical polarization and small values of 6 this is also approxi- 


mately true. In this problem the angles are always small, and so Ris taken to be 
—1. 


Resultant field is A, + 4,—(4y + A,)e ee 
If the free-space polar diagram of the equipment over the range of angles 
concerned is appreciably different from that of an isotropic source, the rays 
should be multiplied by a correcting factor F(@) for the direct ray, F(—6) for the 
reflected ray, and F(A), F(—B), F(—C) and F(—D) for the diffracted rays, 
where F(@) is the amplitude of the voltage polar diagram in the direction 0. 
* The phase is calculated on the basis of the path-length of the geometrical-optics ray (real or 


virtual) for the given group, and not that of the ray which travels via the diffracting edge ; this 


is not of course the actual phase of the diffracted ray, or of the resultant of the group, because of 
the phase change due to diffraction. 
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To calculate the diffracted rays, 
use 1s made of Sommerfeld’s solu- 
tion of the half-plane. In _ his 
solution there are two terms, one 
of which corresponds to simple 
diffraction and the other to diffrac- 
tion and reflection. ‘This second 
term is always small and is ignored. 
Figure 2 is a simplified part of ; 
figure 1, and the diffracted ray SUS 
TSN, is given by Sommerfeld as 


1 Oe 
~ E ue 0S) (Cc is) | 
where (C—2S) is the Fresnel integral function for 


a 2(22 sin i= a) 


The expression is taken positive when the diffracted ray is in the shadow region 
and negative when in the illuminated region, i.e. the sign of the function is the 
same as the sign of w. 

Similar expressions are found for the other rays. ‘The values of w for the 


rays are 
dd \"2 
Ray “/2(5) 


TSN, Sil ae 
Group 1 
TSN; Sy nes 


TSQ,N, sin —~— 
Group 2 
TS’O,.N, sin 


Group 3 
TQ sin 


We 
. C40 
TQ Q,SQ3N, sin OF% 
Group 4 D+6 


L 
( 
L 
be Q,5N, sin ae 
< 
| 
L 
( 
| 
] ros Q.N; sin — 
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$3. METHOD OR COME AVEO AN 
For simplicity of computation the Sommerfeld formula is split into real and 
imaginary parts: 
(1+C) 
Z 


(C—iS)=3[0 + C48) £16 =215)] 


LoS) 
a 


=M-+iN say. 


C and S can be found in tables as functions of u, and table 1 gives M and N as 
functions of wu. As the angles are small, we have approximately : 


7 2d 
ieee 


where F is the angle of diffraction and is given in degrees, or 


v=00247,/4 Aske 


h,—h hp—h 
also A=tan!—= 7 _ B20, 
h,+hp hy +h, 
ets pant te 
(Ce tan 7 D=stan 7 


Thus the method of computing consists of choosing a value of 6, finding the 
corresponding values of w, and then adding real and imaginary parts of the five rays 
in Group | and Group 4 and also the five rays in Groups 2 and 3. If @ is chosen in 
such a way that © 


4crhpsinO nz 


SS) SD 


A 2 


then the phase difference between the two sets of rays is a multiple of 90°. This 
simplifies the addition, as it only involves multiplying the resultant real and 
imaginary parts of Groups 2 and 3 by +1 or +7. The resultant intensity in the 
direction 6 is then the sum of the squares of the real and imaginary parts. This 
process is repeated for different values of 6, thus giving the polar diagram. In 
practice some of the above rays give a negligible contribution to the total amplitude. 
An example is given in the Appendix. 


§4. COMPARISON WITH EXPERIMENT 


Figure 3 gives a set of readings taken for a GL3 with a screen 9 ft. from top to: 
bottom, placed 35 yards from the set, the angle of elevation of the set being 2° 
and the top of the screen 15 ft. above ground. For the purpose of the calculations. 
A was taken as 10cm. and the height of the transmitter as 12 ft. above ground. 
The calculated curve is shown in figure 14. It will be seen that the calculated 
points agree reasonably well with the observed curve, and it therefore appears 
that the approximate method outlined above gives a sufficiently good approxi- 
mation for practical purposes in the field. The experimental curve was obtained 
by Mr. J. T. G. Milne, of the Army Operational Research Group, while tempor- 
arily attached to 21 Army Group. 
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It is important to note that the observations were made in Holland, where 


the ground is exceptionally flat over large areas, and that such agreement may not 
be found on more undulating country. 


Angle of elevation of target (degrees) 


2 3 4 
Observed signal strength 


Figure 3. Signal strength of target echo above a wire screen (observed). 


Angle of elevation of target. (degrees) 


° | 2 3 4 
Relative received signal (voltage) 


‘Figure 4. Signal strength of target echo above a wire screen (theoretical). 
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Wablege 


N 
0-000 
0-012 
0-025 
0-037 
0-050 
0-062 
0-074 
0-086 


0-098 
0-109 
0-120 
Osnsit 
0-142 
0-152 
0-162 
0-172 


0-181 
0-190 
0:198 © 
0-205 
0-212 
0-219 
0-225 
0-230 


0-236 
0-240 
0-242 
0-243 
0-244 
0-244 
0-243 
0-240 


0-236 
0-231 
0-225 
0-219 
0-211 
0-202 
O92 
0-181 


0-169 
0-157 
0-143 
0-128 
0-113 
0-097 
0-080 
0-063 


0-046 
0-029 
0-011 
—0-006 


WwW 
i) 
Un 


SS 
Nu 
St SS) 


400 


ee 
AnnakR RB 
MANoOnN Ub 
Sen Seger 


Bs 


600 
625 
650 
“675 
-700 
25 
*750 
rh) 


800 


_— RRR ee 


\o 0 CO CO 
ONUN 
Sy a ae 


pany 
Ke) 
No 
un 


1EO5 0 
O75 


2-000 
2°025 
2-050 
ALOGS 
2-100 
PENIS 
2-150 
ZA] S 


2-200 
PRINS 
BOY) 
2:275 
2-300 
MESS 
2:350 
2°375 


2-400 
2°425 
2:450 
2Ais 
2:500 
SIDS) 
2°550 
IPSS 
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Table 1 (contd.) 


N 
—0-081 
—0-068 
== (2055 
—0-043 
—0-030 
WO: 
—0-004 

0-021 


0-038 
0-056 
0-066 
0-073 
0-076 
0-076 
0-072 
0-063 


0-055 
0-039 
0-022 
0-006 
—0-010 
—Os023 
(D7 
—0-050 


—0-064 
—0-067 
—0-068 
—0-067 
W057 
—0-042 
—0-026 
== (Oil 


0-004 
0-018 
0-032 
0-045 
0-059 
0-062 
0-062 
0-060 


0-048 
0-032 
0-016 
—0-001 
= (O7/ 
——Or032 
—0-046 
—0-058 


(HY) 
—0-058 
—0-050 
—0-038 


u 
3-900 
3925 
3-950 
SED, 


4-000 
4-25 
4-050 
4-075 
4-100 
4-125 
4-150 
4-175 


4-200 
A225 
4-250 
APIS 
4-300 
A235 
4-350 
A379 


4-400 
4-425 
4-450 
4-475 
4-500 
4-525 
4-550 
4-575 


4-600 
4-625 
4-650 
4-675 
4-700 
4-725 
4-750 
4-775 


4-800 
4-825 
4-850 
4-875 
4-900 
4-925 
4-950 
4-975 


5-000 
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APPENDIX 
Computation of one point 
To illustrate the method of computation, the intensity for one angle of elevation 
is calculated in detail below. The following values of the relevant lengths are 
aised :— 


Transmitter height ae _ een ote 
Distance from screen to transmitter .. 100 ft. 
Height of bottom of screen... ot eet, 
Height of top of screen .. i; ek ores 
Wavelength ie he a Pe wise 
-The four fixed angles are, therefore, 
15 — 12 
— a if ’ 
A=tan™ 79 oh 
12—6 
= fant a ode 
B=tan i007 3°°45 
15 --12 
=f) Bs lien ; 
C=tan 100.7 apa 
12+6 
os af (anaes ee oF ; 
D=tan 1007 = a OD 
d 
and the parameter n =0.0247, /¢ F=0-428F, 
F being the angle of diffraction measured in degrees. For convenience we choose 
our steps of @ such that Fig | 20 
ry ete 
3 


i.e., 0 increases in steps of 0°:2. ‘This interval is sufficiently small for the resulting 
points to give a good indication of the curve. 

As an example, the calculation for 6=2°-8 is given. For this angle there is 
only one geometrical-optics ray, ‘TN, and so we have the following table. 


M N 
G R F 
FORD aed ; ts (from table 1) (from table 1) 
1 TN 0 0 1-000 0 
isn. —1-1 0-0471 —0-240 —0-204 
TSINs 6-25 2-674 0-066 —0:043 © 
2 TSQ:N, 455 1-925 0-114 0-026 
TSO.N: 0-65 0-278 0-355 0-132 
3 TQ,SN, 12:3 5-262* 
TO;GIN: 6°85 2-931 —0-007 0-075 
4 TO;SO,N, 17:9 7:658* 
TO;S OLN; 12-45 5-326* 
a ES ee eee 
1-288 —0-014 


* The table of M and N is only given for 0<u<5. For larger u the diffracted rav is negligible 
owing to the narrow polar curve of the equipment. ; 

The geometrical path of Group-2 rays is 3°54 greater than that of Group 1. As these rays 
undergo a phase reversal on reflection, all groups arrive in phase and the summation above is 
‘correct. Had we taken =2°-4, rays in Groups 2 and 3 would arrive“in opposite phase to those 
in Groups 1 and 4, and would, therefore, be multiplied by —1 before being added. 
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An Isotopic Abundance Rule and its Bearing on the 
Origin of the Nuclei 


The natural abundances of all stable isotopes of any element of higher atomic number 
than 43 are of similar magnitude, except when there exist natural isobars of lower atomic 
number ; those that have isobars of lower charge are comparatively rare. 

The limit at Z=43 is rather arbitrary ; indications of the same type of abundance 
pattern persist down to 3,Se but are not present among light nuclei. “‘ Of similar magnitude” 
means, commonly, within a factor of 2 ; at the worst, within a factor of 11 (88Ba : '%?Ba). 
““ Comparatively rare’ means, generally, between 10 times and 100 times rarer. There 
are only four actual exceptions, i.¢€. nuclei which, possessing lower charge isobars, are more 
abundant than one or more of their isotopes which possess none. ‘These are *49Cd, eShop 
QoNd, *soHg. Of those isotopes which do not possess lower charge isobars, the heaviest is 
occasionally rather less abundant than the rest ; in all such cases it possesses a higher 
charge isobar, which is less rare than usual (the four exceptions mentioned are such cases). 

The isotopic abundance pattern of ;,Xe (table 1) is particularly interesting, since, in 
obeying the above rule, it breaks the usual rule that every odd isotope of even atomic 
number has at least one more abundant adjacent even isotope. 

: The only other exceptions to the latter rule, beyond boron, are 43Sm, 14{Sm, 1%8Pt and 
22. 


Table 1. (From Mattauch, Kemphysikalische Tabellen (Berlin : Springer)) 


5a%e Mass No. 124 126 128 129 130 131 132 134 136 
Abundance % 0:094 0-088 1:91 26:23 4-06 21:18 26-98 10:55 8-95 
Lower charge isobar Yes Yes Yes No Yes No No No No 


In conjunction with Goldschmidt’s (1937) estimates of the cosmic abundances of the 
elements, the first rule may be generalized to the form : 

The total cosmic nuclear abundance at each mass number (above 60) is of the same order 
of magnitude ; where there is more than one stable nucleus of a given mass, the pre- 
dominantly abundant nucleus at this mass is that of lowest charge, and the others may be 
much rarer. 

The only grave exception to the rule in this form is the relative rarity of mass number 
185 (*52Re, possessing no isobar), by a factor between 100 and 1000. It is suggested that new 
determinations of the abundance of rhenium, particularly in meteorites, are called for. 

The rule in this latter form gives evidence in support of the hypothesis that in the origin 
of the nuclei, the heavier ones (as distinct from the lighter nuclei whose abundance is known 
to correspond approximately to a thermodynamic equilibrium (von Weizsacker 1938} 
Chandrasekhar and Henrich 1942)) have been built up by a sequence of neutron capture 
processes. Since the nuclei then assume their final charge chiefly by B-decay, the predomin- 
ance of the lowest charge isobars is a natural consequence. Detailed considerations of the 
néutron condensation process (to be published elsewhere) account also for the general form 
of the abundance curve for the heavier nuclei, and show that this process is the natural 
sequel to the thermodynamic equilibrium state indicated by the cosmic abundances of light 
nuclei. 

F, C, FRANK, 
H. H. Wills Physical Laboratory, 


University of Bristol. 


CHANDRASEKHAR and HENRICH, 1942, Astrophys. F., 95, 288. 
GOLDSCHMIDT, 1937, Skrifter Norske Videnskapsakad., iv. 
VON WEIZSACKER. 1938, Phys. Z., 38, 641. 
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REVIEWS OF BOOKS 


Achievements in Optics, by A. Bouwers. Pp. viii+125. (Brit. Emp. Distr. : 
Cleaver Hume Press, 1947.) 12s. net. : 


This book, one of a series of monographs dealing with progress of research in Holland, 
was prepared for the press before the liberation of the Low Countries, and deals almost 
exclusively with work done—often under conditions of secrecy—during the occupation. 

The first chapter deals with ‘‘ New Optical Systems”’. A clear and simple account 
is given of the advantages possessed by the concave spherical mirror, viz. small spherical 
aberration, complete absence of chromatic aberration, negative Petzval sum and absence 
of all orders of oblique aberrations when used with a stop at the centre of curvature. A 
series of systems is then described, comprising the Schmidt camera, the Maksutovy telescope, 
and derivatives of these two. Very notable results.seem to have been achieved. It is 
interesting to note that yet another claimant to priority for the use of weak negative 
meniscus lens for correcting the aberrations of a concave spherical mirror has appeared. 
This device has been generally attributed, so far, both in America and in this country, to 
Maksutov. 

Chapter II deals with applications of mirror systems to microscopes, telescopes and 
cameras, and contains little that will be unfamiliar to the informed reader. 

Chapter III deals with geometrical aberration theory. An interesting discussion of 
purely geometrical constructions for ray-paths and image points shows that even this 
somewhat dated approach can yet yield useful results. It is followed by a summary of 
the work of Korringa (Thesis, Delft, 1942) and Stephan, and Nijboer (Thesis, Groningen, 
1942), all of which starts from the analytical standpoint introduced by Hamilton and 
developed in England, in particular by T. Smith. The chapter concludes with a 
description of a method of dealing with the light distributions in images in the presence of 
larger amounts of spherical aberration. As with so many of the methods of applied 
mathematics, one can but say that the grave suspicions provoked by the analytical pro- 
cedure are allayed by the apparent success of the method. 

Chapter IV, “ Physical Optics’, summarizes the elegant treatment of aberrational 
diffraction theory given by Nijboer (zbid.) dealing with cases in which aberrational 
asphericities are small; and finally, Zernike’s remarkable analysis of the theory of the 
microscope is given, with particular reference to the phase-contrast method. 

The English used is but rarely foreign to a native ear. The only serious ambiguity 
arises when a full-circular aperture is described as ring-shaped. This would easily be 
taken to mean an annular aperture, and in the context the two have to be carefully dis- 
tinguished. Apart from this, and the use of “ feeded”’ and “ grinded” for the past 
participles of “to feed”? and “to grind” respectively, a very fluent and natural style is 
maintained. 

In sum, this book offers little really new information, but a wealth of scattered and 
useful material has been compressed into it. It can be wholeheartedly recommended to 
specialists and non-specialists alike. Finally, it is no mere courtesy to congratulate the 
Dutch workers who have contributed to optical research with such enviable enthusiasm 
and pronounced success under very unfavourable conditions. H. H. H. 


Microtecnic. International review for measuring and gauging technique, optics 
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